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We review state-of-the-art theory and experiment of the motion of cold and ultracold 
atoms coupled to the radiation field within a high-finesse optical resonator in the dis- 
persive regime of the atom-field interaction with small internal excitation. The optical 
dipole force on the atoms together with the back-action of atomic motion onto the light 
field gives rise to a complex nonlinear coupled dynamics. As the resonator constitutes 
an open driven and damped system, the dynamics is non-conservative and in general 
enables cooling and confining the motion of polarizable particles. In addition the emit- 
ted cavity field allows for real-time monitoring of the particle's position with minimal 
perturbation up to sub-wavelength accuracy. For many-body systems, the resonator 
field mediates controllable long-range atom-atom interactions, which set the stage for 
collective phenomena. Besides correlated motion of distant particles, one finds critical 
behavior and non-equilibrium phase transitions between states of different atomic or- 
der in conjunction with superradiant light scattering. Quantum degenerate gases inside 
optical resonators can be used to emulate opto-mechanics as well as novel quantum 
phases like supersolids and spin glasses. Non-equilibrium quantum phase transitions 
as predicted by e.g. the Dicke Hamiltonian can be controlled and explored in real-time 
via monitoring the cavity field. In combination with optical lattices, the cavity field 
can be utilized for non-destructive probing Hubbard physics and tailoring long-range 
interactions for ultracold quantum systems. 
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I. INTRODUCTION 

Laser light is a versatile tool to cool, prepare and ma- 
nipulate atoms. Laser cooling (Chu, 1998; Cohen Tan- 
noudji, 1998; Phillips, 1998) and optical pumping (Hap- 
per, 1972) relies on spontaneous emission, which is par- 
ticularly important if the laser frequency is tuned close 
to the energy of an atomic transition. It is suppressed if 
the laser frequency is tuned far from any internal excited 
atomic state. In this limit coherent scattering of photons 
dominates and the resulting light force, the dipole force, 
can be derived from an optical potential proportional to 
the laser intensity inducing a Stark-shift. This forms the 
basis for trapping and the manipulation of cold atoms 
(Grimm et ai, 2000), Bose-Einstein condensates (Cor- 
nell and Wieman, 2002; Ketterle, 2002), quantum gases 
(Bloch et ai, 2008; Giorgini et ai, 2008) and mesoscopic 
particles (Gordon and Ashkin, 1980), where spontaneous 
emission has to be avoided. In free space the back-action 
of the particles onto the trapping laser is negligible. In a 
microscopic picture, this means that the probability of a 
photon to be scattered by a particle is so small, that the 
chance for a second scattering event involving the same 
photon is negligibly small. Hence, the modifications of 
the field are not felt by the particles and the light forms 
a conservative optical potential. 

The situation changes drastically when the light field 
is confined in a high-quality optical resonator. Due to 
multiple round trips of intracavity photons not only the 
dipole force gets strongly enhanced, but also the back- 
action of the atoms on the light gets significant. Since 
atomic motion and cavity field dynamics influence each 



other, they have to be treated on equal footing. In most 
cases the dipole force then can no longer be derived from 
a conservative potential (Horak et ai, 1997) and the field 
dynamics get nonhnear (Vukics et ai, 2009). 

To get an intuitive picture, consider for example a mov- 
ing point-like atom, or an entire atomic cloud, forming 
a dielectric medium with refractive index inside a cavity. 
This induces a phase shift on the light field that depends 
on the position and shape of the medium relative to the 
resonator mode structure. Correspondingly, the cavity 
resonance frequency is dynamically shifted with respect 
to the empty cavity. If this shift is comparable to the cav- 
ity linewidth, the cavity field intensity, induced by an ex- 
ternal pump laser, can undergo a resonant enhancement 
and so can the back-action on the motion of the medium. 
For several atoms this coupled atom-field dynamics has 
the character of a long-range inter-particle interaction. It 
also generates a strong nonlinear field response, even if 
the particles are linearly polarizable, as e.g. atoms in the 
low saturation regime. Coupling to further light modes 
gives rise to interference effects, which are the origin of 
collective instabilities and self-organization phenomena. 
Photons leaking out of the cavity cause a damping of this 
coupled dynamics. This designable decay channel can be 
utilized to cool the motion of the medium independent 
of its specific characteristics. 

Historically, cavity quantum electrodynamics (QED) 
was born as a research field devoted to studying the ra- 
diation properties of atoms when boundaries are present 
(Berman, 1994; Haroche, 1992; Purcell, 1946). Advances 
in cavity technology over more than 30 years allowed 
to reach, both in the microwave (Raimond et ai, 2001; 
Walthcr, 2002) and optical (Kimble, 1998; Mabuchi and 
Doherty, 2002) frequency domains, the strong coupling 
regime where the coherent interaction between an atomic 
transition and a single radiation field mode dominates 
over all dissipation processes. As a next step, cold and 
slow atoms have been integrated successfully within op- 
tical cavity QED experiments, which led to significant 
coupling of the atomic motion to the cavity field. It 
got possible to generate sufficiently strong forces in order 
to trap an atom in the field of a single photon (Hood 
et at, 2000; Pinkse et ai, 2000). Several experiments 
achieved strong coupling even in the dispersive regime of 
cavity QED where the detuning between the light field 
and the internal atomic transitions is large. Although 
the resonant energy exchange between atom and field is 
suppressed in this regime, the position-dependent cavity 
frequency shift exceeds the cavity linewidth. Motion- 
induced changes of the effective resonator frequency and 
its back-action on mechanical motion is also the physical 
ground of cavity optomechanics (Kippcnberg and Vahala, 
2008) , which can be considered as an extension of disper- 
sive cavity QED towards macroscopic objects. 

In this review we survey the recent advancements of 
cavity QED systems in which coherent momentum ex- 
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change between particles and radiation field is the domi- 
nating effect of the light-matter interaction. The external 
degree of freedom of the material component ranges from 
the center-of-mass motion of a single atom, or a cloud of 
cold atoms, to the density distribution of a continuous 
medium such as the quantized matter-wave field of an 
ultracold gas. 

The review illuminates different generic features of the 
cavity-generated optical dipole force and is structured 
in three main sections. Briefly summarizing, in Sect. II 
we discuss the consequences of the retardation between 
atomic motion and the cavity field dynamics. This time 
delay leads to an irreversible dynamics that can be the 
basis of cooling schemes, as presented for single atoms 
in a cavity. Sect. Ill discusses how the field modification 
induced by an atom acts back on the motion of other 
atoms moving within the cavity. This cavity-mediated 
atom-atom interaction is a source of collective effects in 
atomic clouds. Finally, in Sect. IV we consider the col- 
lective dynamics of an ultracold gas induced by its strong 
coupling to the cavity field. Owing to the low temper- 
ature, the dynamics involves a reduced set of motional 
degrees of freedom, and the system becomes a realization 
of various paradigmatic models of quantum many-body 
physics and quantum optics (Lewenstcin et al, 2007). 

The most elementary situation that we will discuss 
is the dispersive atom-field dynamics of a single atom, 
or polarizable particle, inside a laser-driven high-finesse 
cavity (Domokos and Ritsch, 2003; Pinkse and Rcmpc, 
2002). The cavity field dynamically responds to the po- 
sition and velocity of the particle, thereby generating a 
time-dependent dipole force acting back onto the parti- 
cle motion. It is the finite response time of the cavity 
field which gives rise to the velocity dependent compo- 
nent of the force. It can have the character of a friction 
force shuffling kinetic energy from the particle to the 
cavity field and dissipating it via the cavity loss chan- 
nel. (Vuletic and Chu, 2000). This allows for cooling 
and self-trapping of particles in the cavity field (Maunz 
et ai, 2004; Nussmann et at, 2005). Sub-recoil cavity 
cooling of an ultracold atomic cloud has been achieved 
recently (Wolke et at, 2012), which paves the way to- 
wards reaching quantum degeneracy without relying on 
evaporative cooling techniques. Cavity cooling allows for 
slowing of any sufficiently polarizable particle with small 
absorption, without the need of a cycling transition. The 
possibility of extending the applicability of cavity cooling 
beyond atoms has been the subject of extensive research 
in the past years. The light field leaking out the cav- 
ity carries information on the trajectory of the particle 
(Hood et ai, 1998; Maunz et ai, 2003). Continous mon- 
itoring of the atomic motion, in turn, can be used for 
feedback control (Fischer et ai, 2002), which became by 
today the standard tool to capture single atoms inside a 
cavity for quantum manipulation (Kubanek et al, 2009). 

For cold atomic ensembles inside a laser-driven cavity 



(Elsasser et al., 2004; Kruse et ai, 2003b) the atom-field 
coupling strength increases and the dynamics becomes 
more complex. In many cases the effective coupling 
strength between particles and cavity field scales with the 
square root of the particle number (Raizen et al, 1989; 
Sauer et at, 2004; Tavis and Cummings, 1968; Tuchman 
et ai, 2006). As a consequence of this, the cooling of the 
center-of-mass motion is correspondingly more effective. 
Additional complexity arises from the relative motion of 
the particles, as the local intensity of the cavity field ex- 
perienced by one atom depends on the position of all 
other atoms (Horak and Ritsch, 2001c). This gives rise 
to an effective long-range (Miinstcrmann et ai, 2000), 
or global atom-atom interaction, described by an overall 
dispersive shift. The contribution of each particle to this 
shift depends on the local field intensity, which is pro- 
portional to the square of the cavity mode function at 
the position of the atom. In the low excitation regime 
this can be captured by a collective potential. Further- 
more, dissipative forces acting on the relative motion of 
the particles have been identified (Chan et ai, 2003) and 
interesting correlations between particles can build up 
(Asboth et ai, 2004). 

The cavity-meditated long-range interactions have a 
different character when the mode of the driving field 
is not identical with the cavity mode. In this case, 
the atoms can be considered as sources for the intra- 
cavity field and interference between these sources be- 
comes crucial. Correspondingly, the effective cavity driv- 
ing strength depends on the position of all atoms within 
the cavity mode profile and it is the field amplitude rather 
than its intensity, which mediates the long-range interac- 
tion. For the case of a transversally laser-driven atomic 
ensemble in a linear cavity, the long-range interaction 
causes a phase transition to a self-organized phase, in 
which the atoms arrange themselves in a checkerboard 
pattern, thereby maximizing scattering into the cavity 
mode (Black et ai, 2003; Domokos and Ritsch, 2002). 
In a unidirectionally driven ring cavity geometry, col- 
lective scattering between the two counter-propagating 
cavity modes results in a collective instability, referred 
to as collective atomic recoil lasing (Kruse et ai, 2003a). 
Various mean-field type theories can be used to describe 
the non-equilibrium dynamics and asymptotic behavior 
of large atomic ensembles, including the derivation of 
scaling laws characterizing the above described critical 
phenomena (Asboth et ai, 2005; Griefier et ai, 2010). 

Coupling ultracold atomic ensembles or Bose-Einstein 
condensates to the radiation field inside a high-finesse 
resonator, requires a quantized description of the atomic 
motion and reduces the number of relevant external de- 
grees of freedom (Brennecke et ai, 2007; Colombe et ai, 
2007; Gupta et ai, 2007; Slama et ai, 2007a). In the 
case of a laser-driven cavity, situations can be realized 
where the cavity field couples dominantly to a single col- 
lective motional mode of the atomic ensemble, provid- 
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ing a direct analogy to cavity optomechanics (Brcnncckc 
et al, 2008; Murch et al, 2008; Stamper-Kurn, 2012). 
Coupling a laser-driven Bose-Einstein condensate to the 
vacuum field of a cavity, leads to a quantum phase tran- 
sition between a superfluid and a self-organized phase 
(Baumann et al., 2010; Nagy et at, 2008). This provides 
an open-system realization of the Dicke Hamiltonian and 
its quantum phase transition (Dickc, 1954; Dimcr et al., 
2007; Hepp and Lieb, 1973; Nagy et al., 2010). The self- 
organized state can also be considered as a supersolid 
resulting from a broken Ising-type symmetry. More com- 
plex situations occur in highly degenerate multi-mode 
cavities (Gopalakrishnan et al., 2009, 2011b; Strack and 
Sachdev, 2011). 

Ultracold gases in optical lattices are one of the most 
intriguing systems in which the power of atomic and laser 
physics can be exploited to explore generic phenomena 
of solid-state physics (Bloch et at, 2008). The Hubbard 
model describing the dynamics of periodically arranged 
bosons or fermions can be de facto realized with ad- 
justable parameters and variable dimensionality. When 
the optical lattice potential is created by the field sus- 
tained by an optical high-finesse cavity, the correspond- 
ing cavity Hubbard-model predicts exotic new phases of 
matter (Larson et al., 2008a; Maschler and Ritsch, 2005). 
In many cases the cavity fields provide for a convenient, 
built-in real-time observation tool. Analyzing the emit- 
ted fields allows for dynamical monitoring of quantum 
phase transitions with minimum and well controlled mea- 
surement back-action (Mekhov et al., 2007c). 

II. SINGLE ATOMS IN A CAVITY 

A central objective of cavity quantum electrodynam- 
ics (QED) is the perfect control of light-matter inter- 
action at the single-atom and single-photon level in the 
regime of strong coupling where atom and cavity field 
form a single entity. A long lifetime of such an 'atom- 
photon molecule' requires slow and very cold atoms to 
ensure long interaction times and precise control of the 
atomic position. At sufhciently small kinetic energies, 
however, the light forces induced by even a few intra- 
cavity photons influence the atomic trajectory. The first 
cavity QED experiments with cold atoms have already 
manifested that the cavity light forces guide or deflect 
slowly moving atoms. In addition, extra diffusion takes 
place in cavity-sustained dipole traps which may remove 
the atom from the interaction volume. Clear signatures 
of such effects have been observed in transmission spec- 
troscopy experiments (Hood et al., 1998; Mabuchi et al., 
1996; Miinstermann et al., 1999). Time-resolved detec- 
tion of the transmitted light signal allowed for the re- 
construction of atomic trajectories (Hood et al., 2000; 
Pinksc et al., 2000). These experiments set the stage to 
include the atomic center-of-mass degrees of freedom and 



the optical forces in the cavity QED theory. In the fol- 
lowing decade, the theoretical and experimental efforts 
resulted in an extension of the interaction time from the 
transit-time range of microseconds to the range of min- 
utes (Figueroa et al., 2011; Kubanek et al., 2011). 

A. Mechanical effects of light on atoms in a cavity 

The theoretical description of the coupled atom-field 
dynamics has been presented in detail by Domokos and 
Ritsch (2003). Here, we recapitulate the notations and 
methods. Within the vast field of single-atom cavity 
QED, we restrict ourselves to the atomic motion in a 
cavity, in particular, to the important concept of cavity 
cooling. We review the recent experiments demonstrat- 
ing cavity cooling of single atoms. It is a manifestation 
of the time-delayed action of the electric dipole force on 
atoms within the cavity. The understanding at a single- 
atom level complements nicely another facet of cavity 
cooling which we will encounter in the case of many atom 
systems, where it appears in the form of the imaginary 
part of the collective excitation spectrum. 

1. A two-level atom in a cavity 

We consider a single two-level atom with transition 
frequency uja coupled to a single mode of the electro- 
magnetic field inside an optical resonator with resonance 
frequency loc- These frequencies will be referenced to 
the frequency uj of an external pump laser by defining 
the cavity detuning Ac = uj — loq and the atomic detun- 
ing Is. A — Lu — LJA- The two relevant atomic states are the 
ground state \g) and the excited state |e). We introduce 
the atomic raising and lowering operators, = |e) {g\ 
and a = \g) (e|. The cavity mode variables are the pho- 
ton creation and annihilation operators, and a, re- 
spectively. In the electric-dipole and the rotating-wave 
approximations and in a frame rotating at the angular 
frequency lj, the atom-field coupling is described by 

Hic/n = -Aca^a-A^(r)aV+^g (a^a/(r) - /*(r) aV) , 

(1) 

which is usually referred to as the Jaynes-Cummings 
Hamiltonian (Jaynes and Cummings, 1963) and, in the 
quantum optical context, has been reviewed by Shore 
and Knight (1993). The emphasis here is on that the 
position r of the atom is explicitly taken into account. 
The spatial dependence of the atomic detuning, A^(r) — 
A^ — A5(r), may account for a differential AC-Stark shift 
A5(r) which can be induced by auxiliary, far-detuned 
optical trapping fields. The coupling strength in Eq. (1) 
is spatially modulated according to the intracavity elec- 
tric field strength which is proportional to the cavity 
mode function /(r). For the effects reviewed in this pa- 
per, it is sufficient to consider modulations on the op- 
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tical wavelength scale, thus writing /(r) = cos{kx) for 
a standing-wave mode of a Fabry-Perot resonator, or 



for the running-wave modes sustained by 



a ring resonator {k = uj/c\s the optical wavenumber). 
The maximum coupling strength is given by the single- 
photon Rabi frequency, g — dyjhuic /'^^oV , where d is 
the atomic dipole moment along the cavity mode polar- 
ization and V = / (i'^r|/(r)p denotes the effective cavity 
mode volume (the maximum of |/(r)| is set to 1). The 
rotating wave approximation relies on that the character- 
istic frequencies of i/jc are much smaller than the optical 
frequency, i.e. , (|A^|, |Ac|,.g ^ t^). The atomic center- 
of-mass (CM) motion is a dynamical component of the 
system, which is described by the Hamiltonian 



i?moch = ^ h Kl(r) , 

2m 



(2) 



where m is the mass of the atom and the term Vd can 
represent an arbitrary external trapping potential. For 
the case of a far off-resonance optical dipole trap, this 
term, together with the differential AC-Stark shift A5'(r) 
in Eq. (1), fully describe the effect of the trapping laser. 
The characteristic frequency of the CM motion is given 
by the kinetic energy of an atom carrying one unit of 
photon momentum, |p| = fik. We will use throughout the 
paper the notion of recoil frequency (Cohen- Tannoudji, 
1992), with the notation 



2m 



The system can be excited with a coherent laser field 
at frequency lo, which either drives the cavity mode with 
driving amplitude 77, or directly the atomic internal de- 
gree of freedom at Rabi frequency 17, described by 



Hpum-p/fi = "iv — aj + iil h{r) (cr' — cr 



(3) 



For the case of pumping the atom with a standing-wave 
laser field from a transverse direction perpendicular to 
the cavity axis, the spatial mode function is given by 
h{r) = cos(fcz). -ffpump is effectively time-independent 
since we work in the frame rotating at the angular fre- 
quency uj of the monochromatic pump laser. 

Cavity QED systems in the optical domain are strongly 
influenced by dissipative coupling to the vacuum modes 
of the electromagnetic field environment (thermal pho- 
tons can be neglected at optical frequencies). Corre- 
spondingly, the dynamics of the system is described by a 
quantum master equation (Carmichacl, 2003) 



[H,p] +C 



(4) 



with H = Hjc + i/moch + ^^pump ^nd p denoting the 
density operator for the atomic (motional and internal) 
and cavity degrees of freedom. The dissipative processes 



are captured by the Liouville operators in Born-Markov 
approximation 

CcavP — {a^ ap + pa' CL — 2apa^^ , (5a) 

describing decay of the cavity field at rate k, and 

-2 j d^u iV(u)CTe-*'^-^">e*'^'^"V^) , (5b) 

describing spontaneous decay of the excited state |e) at 
rate 7 accompanied by the emission of a photon into 
the free-space modes of the electromagnetic field environ- 
ment. This process involves a recoil of kA = u>a/c « k op- 
posite to the direction u of the emitted photon, which is 
averaged over the directional distribution function iV(u) 
characterizing the given atomic transition. 

In general, the full quantum dynamics of the system 
defined by Eq. (4) including all degrees of freedom - the 
CM motion, the internal electronic dynamics and the cav- 
ity photon field - cannot be solved analytically even for 
a single atom. 

2. Dispersive limit 

For a broad class of cavity QED parameters, atomic 
saturation effects are negligible and the atoms can be 
considered as linearly polarizable particles. This holds 
true when the internal atomic variables cr, evolve on a 
much faster time scale as compared to the other variables 
due to a large atomic detuning IS. a or a large spontaneous 
decay rate 7. In either case, following the usual technique 
of adiabatic elimination, the atomic polarization operator 
a can be 'slaved' to the cavity mode and atomic position 
'master' variables. In the absence of direct atom driving, 
i.e. , = in Eq. (3), one obtains 



-i^A + 7 



(6) 



This approximation is valid if the population in the ex- 
cited atomic state is negligible (low saturation regime). 
By inserting the slaved variable a into iJjc and into the 
Liouville operator Eq. (5b), an effective master equation 
is obtained. Of particular interest is the large detuning 
limit in which the CM motion and the cavity mode are 
coupled dispersively by 



a' a . 



(7) 



It captures, on the one hand, the atom-induced disper- 
sive shift of the cavity mode resonance frequency which 
depends on the momentary position of the atom. On the 
other hand, the cavity field gives rise to an optical po- 
tential cx |/(r)p felt by the atom whose depth depends 
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on the dynamical photon number. Dissipation can be 
treated analogously and the effective Liouville operator 
was presented by Domokos et al. (2001). The dispersive 
and absorptive effects of the atom are expressed in terms 
of the parameters 



Uo = 



To = 



(8a) 
(8b) 



respectively. These relations reveal the connection be- 
tween the cavity QED parameters and the complex sus- 
ceptibility X = x' — ix" of a linearly polarizable object 
with electric polarization P = eoxE. With this con- 
nection at hand, the theory presented in this section 
can be used to describe a much broader class of parti- 
cles than only two-level atoms, and most of the findings 
can directly be applied to polarizable particles of sub- 
wavelength size. In Sect. II. C. 3, the linear polarizability 
picture is refined for the case of molecules. 

Using the dispersive interaction Hamiltonian, Eq. (7), 
the quantized one-dimensional motion of a single atom 
strongly coupled to single-mode cavity field has been nu- 
merically simulated (Vukics et at, 2005). The calculation 
confirmed the basic assumption of semiclassical theories 
(see Sec. II. A. 3), stating that the coherence length of 
the atomic wavefunction reduces well below the optical 
wavelength after a few irreversible scattering events. This 
happens although in the dispersive limit the coupling to 
the environment is provided by cavity photon loss rather 
than spontaneous photon scattering into free space. An 
efficient numerical code has been developed providing a 
general framework for Monte-Carlo wavefunction simula- 
tions of systems composed of the 'quantum optical tool- 
box' (Vukics, 2012; Vukics and Ritsch, 2007). 

If the atom is laser-driven from a transverse direction, 
i.e. , 7^ in Eq. (3), the adiabatic elimination of the 
internal degrees of freedom leads to 



gf{v)a + nh{r) 



-iA, 



■7 



(9) 



Consequently, additional terms appear in the effec- 
tive adiabatic Hamiltonian Eq. (7) and the Liouvillean 
Eq. (5b). In particular, coherent photon scattering be- 
tween the transverse laser field and the cavity mode gives 
rise to the effective cavity pump term 



H„/h = r^eff h{r) (/*(r) + f{r)a) 

with the effective cavity drive amplitude 

AAgn 



(10) 



Ai + 7^ 



The atomic recoil accompanied by photon scattering is 
accounted for by the spatial dependence of this term. 



3. Semiclassical description of atomic motion 

In many cavity QED experiments, cold atoms are re- 
leased from a magneto-optical trap into the resonator vol- 
ume. As the temperature T of the atoms is well above the 
recoil temperature, ksT ^ hujj^, where fcs is the Boltz- 
mann constant, one can assume that the reduced density 
matrix is almost diagonal both in position and momen- 
tum representation. This allows to treat the position r 
and momentum p of the atom as stochastic c-number 
variables. The regime of ultracold atoms ksT < tkun 
wiU be treated in Sec. IV. 



a. Langevin-type equation The separation of the quan- 
tized internal and the classical motional degrees of 
freedom was developed for the description of laser- 
cooling (Cohen-Tannoudji, 1992; Dalibard and Cohcn- 
Tannoudji, 1985; Gordon and Ashkin, 1980). This ap- 
proach has been adopted to the cavity QED scenario by 
extending the internal degrees of freedom to the com- 
bined space of the atomic polarization and cavity mode 
(Hcchcnblaikner et al., 1998; Horak et ai, 1997). By 
eliminating the internal degrees of freedom, the dynam- 
ics of the CM variables can be formulated in terms of a 
stochastic differential equation 



r = 



p = f + /3 



P , „ 



(11a) 
(lib) 



where f denotes the classical force and f3 a friction co- 
efficient in a non-conservative and velocity dependent 
force term. In general, (3 can be a tensor in the three- 
dimensional space as atomic motion along any direction 
gives rise to friction in all three spatial directions (Vukics 
et al., 2004). When the eigenvalues of the tensor /3 (or 
scalar in ID) are negative, one encounters cavity cooling. 
The noise term H induces the stochastic behavior. It 
has vanishing mean value and is defined via the diffusion 
matrix D according to 



{sit)om')) = 'DS{t-t'), 



(12) 



where o denotes the dyadic product. The exact noise 
correlation function has a width in the range of the dis- 
sipative parameters k and 7 of the internal dynamics. 
Therefore, it can be approximated by a Dirac-5 only 
on the much slower CM motion time scale set by the 
inverse of the recoil frequency, 0}]^^ 3> mm{K~^ ,"f~^}. 
The method of calculating the c-number parameters f , 
/3, and D of the Langevin-type equation from the master 
equation concerning the internal degrees of freedom was 
presented by Hechenblaikner et al. (1998) for the one- 
dimensional, and by Domokos and Ritsch (2003) for the 
three-dimensional case. This method accounts for the 
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quantum effects of the internal dynamics; hence the full 
approach is semiclassical. 

The practical use of this method is strongly limited: 
the nonlinear quantum master equation for the internal 
and cavity degrees of freedom has to be solved numeri- 
cally and for all atomic positions r. Moreover, the Hilbert 
space of the photon field has to be truncated at low pho- 
ton numbers. This approach was adopted by Doherty 
et al. (2000) and Fischer et al. (2001) to simulate the ex- 
periments conducted by Hood et al. (2000) and Pinksc 
et al. (2000). 

Analytical approximations can be obtained for the low 
atomic saturation regime (Murr, 2003) , where the atomic 
polarization can be replaced by a bosonic operator and 
hence the internal dynamics is described by linear equa- 
tions of motion. It is then possible to calculate the fric- 
tion coefficient /3 and, corresponding to its sign, the cool- 
ing versus heating regions can be mapped as a function 
of the detunings and A^, as shown for example in 
FiK. 1. 
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FIG. 1 (Color online) Cooling and heating regions as a func- 
tion of atom and cavity detunings. Shown are contour plots 
of the friction coefficient /3 (averaged over an optical wave- 
length) acting along the cavity axis on a laser-driven atom. 
Left: Bad-cavity regime, g = 7/2, k = IO7; Right: Good 
cavity regime, g = 87, «: = 7, where the dressed-state picture 
can be invoked for interpretation. Blue contour lines indicate 
cooling (C, P < 0), red ones heating (H, /3 > 0) regions. Note 
that the spatially averaged friction coefficient is shown here; 
tightly confined atoms locaHzed within a small fraction of the 
wavelength can follow completely different behavior depend- 
ing on their position. 



b. Semiclassical theory in the dispersive limit In the dis- 
persive limit of atom-cavity coupling presented in Sect. 
II. A. 2, there is an alternative semiclassical approach 
(Domokos et at, 2001). The Wigner quasi-probability 
distribution function can be defined in the joint phase 
space of the atomic CM motion and the cavity field am- 
plitude. The quantum master equation translates then 
into a partial differential equation for the Wigner func- 
tion. By dropping all terms containing higher than 
second-order derivatives, the resulting Fokker-Planck 



equation corresponds to the evolution of classical stochas- 
tic variables associated with the atomic motion and the 
cavity field. One can consider this approach as the con- 
struction of a semiclassical model which lies closest to the 
true quantum dynamics. As compared to Eq. (11), here 
large intracavity photon numbers are allowed, in fact, the 
validity of this approach requires photon numbers larger 
than 1. 

Consider the generic example of a single atom moving 
in one dimension along the axis of an externally driven 
linear cavity with the mode function f{x) = cos{kx), 
described by the equations (Domokos et al., 2001) 

P 



p=-hUo\a\'^f{x)+^,, 
a^ri~i (Uof^{x) - Ac) a 



{K + Tofix))a + ^^^ 



(13a) 
(13b) 

(13c) 



Apart from the noise terms and ^q, , these equations co- 
incide with the classical description in the initial cavity 
cooling paper by Horak et al. (1997). The force in (13b) 
acting on the atom is formally identical to the gradient of 
the optical dipole potential of the cavity mode. The am- 
plitude a, however, depends not only on the momentary 
position of the atom but has a memory effect because of 
the finite bandwidth k of its linear response. The actual 
force has therefore a velocity-dependent character which 
can give rise to a viscous friction force, that is, cavity 
cooling. Within this approach the friction cannot be de- 
termined by a single coefficient (3, on the other hand, the 
friction effect is correctly described for arbitrary velocity. 

The noise sources are taken into account in a consis- 
tent way as imposed by quantum mechanics. This results 
in non-trivial correlations (^pCa) 7^ 0. The result for the 
diffusion matrix as well as the generalization for several 
atoms was published by Asboth et al. (2005). The general 
model is used for numerically studying many-body sys- 
tems, c.f. Sect. Ill, well above the temperature of quan- 
tum degeneracy. 



c. Scattering model The semiclassical Langevin-type 
equation (11) can be constructed without mode decom- 
position of the radiation field. This approach is required 
when, instead of a simple Fabry-Perot-type cavity geom- 
etry, one considers an atom interacting with the radiation 
field of an interferometer which is composed of an ar- 
bitrary one-dimensional configuration of beam splitters. 
To deal with this situation a scattering model has been 
established by Xuereb et al. (2009) and solved for the 
force terms acting on a particle as in Eq. (11) by Xuereb 
et al. (2010). In the scattering model the atoms and 
beam splitters are treated on an equal footing as 'scat- 
terers' characterized by a single polarizability parame- 
ter. Thereby, a unified framework is created to describe 
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opto-mechanical systems in general, which has revealed 
the close relationship between cavity cooling of atoms 
and radiation-pressure cooling of mirrors (Arcizet et ai, 
2006; Gigan et al., 2006; Kleckner and Bouwmeester, 
2006; Metzger and Karrai, 2004; Schhesser et al., 2009). 



B. Cavity cooling 

One of the most promising results from the under- 
standing of the complex cavity QED dynamics involving 
atomic motion is the realization of cavity cooling, i.e. , 
the dissipation of kinetic energy through the cavity pho- 
ton loss channel in a controlled manner. 

Early ideas about using an optical resonator to en- 
hance the efficiency of laser-cooling relied on the modi- 
fication of the spectral mode density of the electromag- 
netic radiation field in the presence of spatial boundary 
conditions (Lcwcnstein and Roso, 1993; Mossbcrg et ai, 

1991) . In the most general form, the notion of cavity 
cooling in the perturbative regime has been expressed by 
Vulctic and Chu (2000). If an atom, placed inside an 
optical cavity, is laser driven at a frequency below the 
cavity resonance, Ac < 0, scattering favors the emis- 
sion of photons at frequencies higher than the pump fre- 
quency due to the increased mode density around the 
cavity resonance. The energy needed to upshift the pho- 
ton frequency is provided by the loss in kinetic energy 
in processes of inelastic scattering. With this very sim- 
ple picture, a robust three-dimensional cooling effect can 
be interpreted (Vulctic et ai, 2001). However, the pic- 
ture holds only true in the regime of weak atom-photon 
coupling, see the left panel of Fig. 1. When the reabsorp- 
tion of a photon starts to become non- negligible, which 
happens in a high-finesse cavity, the cooling mechanism 
substantially changes. This drastic change is illustrated 
in Fig. 1, where the right panel presents the friction co- 
efficient for a ratio g/n in the the single-atom strong- 
coupling regime. In the case of a standing-wave cavity, 
the dynamical cavity cooling effect can be interpreted in 
the frequency domain by means of a Sisyphus-type argu- 
ment (Horak et at, 1997), using the dressed-state picture 
of the strong-coupling regime of cavity QED (Harochc, 

1992) , see Fig. 2. For the case of a ring cavity, inter- 
estingly, the intuitive photon scattering picture can be 
pursued also in the strong-coupling regime and the full 
velocity dependence of the radiation pressure for arbi- 
trary coupling constant g can be obtained (Murr, 2006). 

In the following we will survey two regimes where cav- 
ity cooling has been demonstrated experimentally, and 
present the corresponding intuitive pictures of the cool- 
ing effect. Both regimes are in the dispersive limit of 
atom-photon interaction keeping the atomic saturation 
low. 





-1 1 
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FIG. 2 Sisyphus-type cooling mechanism underlying the 
hyperbolic-shaped cooling region in the right panel of Fig. 1. 
The atomic motion leads to a modulation of the internal 
dressed-state energy levels |±) which are linear combina- 
tions of the \g,l) and |e, 0) states with mixing determined 
by the mode function f{x) = cos{kx). This amounts to a 
state-dependent potential for the CM degree of freedom. For 

2 

A,4 < and Ac ~ —k + the transition from the ground 
state \g, 0) to the lower dressed state |— ) is resonantly excited 
at an antinode. Thus the excitation happens more likely at 
the minimum of the potential wells, whereas spontaneous or 
cavity decay transfers the atom-cavity system back to the 
ground state homogeneously in space. From Doniokos and 
Ritsch (2003). 



1. Cavity cooling with blue-detuned probe light 

Cavity cooling of single atoms has been first demon- 
strated by Maunz et al. (2004) via the observation of 
extended storage times and improved localization of sin- 
gle ^^Rb atoms in an intracavity dipole trap. The trap 
field was red-detuned with respect to the atom, however, 
the cooling was induced by a weak, blue-detuned probe 
field. The cooling rate has been estimated to exceed that 
achieved in free-space cooling methods by at least a fac- 
tor of five, for comparable excitation of the atom. Maunz 
et al. (2004) present a very intuitive interpretation of the 
cooling effect in terms of the classical notion of the refrac- 
tive index. Consider a standing-wave optical cavity res- 
onantly excited by a weak probe laser, Ac — 0, which is 
blue-detuned from the atomic resonance by Aa = 27r x 35 
MHz > 0, see Fig. 3. The resulting light shift parameter, 
Eq. (8a), far exceeds the cavity linewidth, Uq > 5k, so 
that even one atom can significantly infiuence the optical 
path length between the cavity mirrors. Due to Eq. (7), 
the atom placed at a node of the standing-wave mode 
profile does not couple to the cavity field and the intra- 
cavity intensity is maximum. In contrast, placed at an 
antinode, the atom shifts the cavity resonance towards 
higher frequency, i.e. , out of resonance with the probe 
laser, resulting in a reduced intracavity intensity. In a 
high-finesse cavity, however, the intensity cannot drop 
instantaneously when the atom moves away from a node. 
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The induced blue-shift of the cavity frequency at almost 
constant photon number leads to an increase of the en- 
ergy stored in the field, before the photons are able to 
leak out of the cavity. This occurs at the expense of 
kinetic energy of the atoms. The reverse, accelerating 
effect occurring when the atom moves from an antinode 
towards a node is much weaker, because the cavity is 
initially out of resonance with the probe laser and conse- 
quently only a small number of photons are present and 
undergo a corresponding red-shift. This argument also 
reveals that the delicate correlation between the atomic 
motion and the photon number variation, underlying the 
cooling effect, imposes an upper bound on the atomic ve- 
locity, kv < K, which sets the velocity capture range of 
cavity cooling. 
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FIG. 3 Experimental scheme used for the observation of cav- 
ity cooling. Single atoms are captured by an optical dipole 
trap formed by far red-detuned light in a longitudinal cavity 
mode which is different from the one used for cavity cooling. 
The characteristic parameters of the interaction between the 
weak probe and the atoms are {g,K,y) = 2tt x (16,1.4,3) 
MHz. From Maunz et al. (2004). 

In the experiment, single atoms injected into the cavity 
are trapped at the field antinodes of a strong intracavity 
dipole trap. To be detectable in cavity transmission of 
the weak probe beam, the atoms simultaneously have to 
be close to an antinode of the probe field mode. As the 
probe field induces cavity cooling, the resulting stronger 
confinement can be directly read out of the transmitted 
signal, as shown in Fig. 4. Time-resolved detection of 
the cavity transmission allowed to extract a cooling rate 
of /3/m = 21kHz, which is large compared to the es- 
timated coohng rate of 4 kHz expected for blue-detuned 
Sisyphus cooling of a two-level atom in free space, or with 
the Doppler cooling rate of 1.5 kHz at equivalent atomic 
saturation. 



2. Cavity cooling and trapping with far red-detuned light 

Far-off-resonance dipole traps are commonly used for 
long-time capturing and localization of neutral atoms 
(Grimm et al., 2000). The suppression of spontaneous 
emission results in an almost conservative trapping po- 
tential. However, with the elimination of spontaneous 
emission (|Ayi| 3> 7), any free-space cooling mechanism 
also disappears. The far-off-resonant trapping scheme 
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FIG. 4 Demonstration of cavity cooling. Time-resolved re- 
duction of the averaged cavity transmission of a weak resonant 
probe beam indicating improved localization of the atom at 
the antinodes of a far red-detuned trapping field. The closer 
the atom resides at an antinode, the larger the detuning of the 
cavity resonance with respect to the probe frequency. From 
Maunz et al. (2004). 



has been revisited for a strongly coupled atom-cavity sys- 
tem where the cavity mode provides a new dissipation 
channel. 

Surprisingly, cavity cooling can remain very efficient in 
the limit of large atomic detuning |A^| — >■ 00. For opti- 
mal cooling the driving frequency has to be set slightly 
below the cavity resonance frequency ~ —n -f Uq. 
Underlying the cooling mechanism is a polariton reso- 
nance of the strongly coupled atom-cavity system (corre- 
sponding to the dressed state |— ) in the weak excitation 
limit where only the lowest excitations manifold of the 
Jaynes-Cummings spectrum matters, see Fig. 2). Even if 
uja and loc are very different, the bare cavity resonance 
is slightly modified because the photon excitation mixes 
with a small amount of the atomic excitation. In an inho- 
mogeneous system, the mixing leads to a dependence of 
the polariton resonance on the atomic position (see Fig. 
2). Although the modulation is tiny in amplitude, the 
resonance is comparably narrow, having a width in the 
range of k for the cavity-type polariton. Thus the sys- 
tem can be very sensitive to the atomic motion and even 
slow atom velocities induce large non-adiabatic modula- 
tions of the steady-state field amplitude (Domokos et al., 
2004; Murr, 2006). 

For demonstration, we consider the simplest case of an 
atom moving along the cavity axis in the field generated 
by an external driving laser. It has been shown that the 
standing- wave cavity field simultaneously traps and cools 
the atom (Vukics and Domokos, 2005). For a standing- 
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wave mode, f{x) = cos{kx), the cooling rate is given by 
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4sin^ {kx)x 



2g'^{Ac - Uq cos^ {kx)){K + Tq cos^ (kx)) 
({Ac - Uo cos2 {kx)f + (k + To cos2 {kx)f^ ' 



(14) 



where Pe denotes the mean population in the excited 
state |e). Choosing the cavity detuning as Ac ~ —k + Uq 
leads to the optimum friction coefficient which, spatially 
averaged, reads 



/3 



7 \K 



(15) 



On the left-hand side the cooling rate is normalized to the 
rate of spontaneous photon scattering. This expression 
shows that the friction coefficient at a fixed saturation Pg 
is independent of the atomic detuning, which gives rise 
to the perspective of cooling molecules or other objects 
without closed cycling transition (see Sec. II. C. 3). 
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FIG. 5 Transverse pump scheme of cavity QED. Atoms are 
transported into the cavity using an optical conveyor belt. 
Instead of driving the cavity directly, the atoms are trans- 
versely laser-driven giving rise to photon scattering into the 
cavity mode. The standing-wave dipole trap yields a large 
differential AC-Stark shift, i.e. , a modulation of the atomic 
detuning AA(r), see Eq. (1). In this geometry the cavity 
vacuum field, the weak driving laser tuned according to the 
optimum choice in Eq. (15), and the trap laser together form 
a very efficient three-dimensional cooling scheme. From Nuss- 
mann et al. (2005). 



contribute to the cooling efficiency. The resulting cool- 
ing effect has been demonstrated experimentally by Nuss- 
mann et al. (2005) , making use of an orthogonal arrange- 
ment of a cooling laser, a trapping laser and a cavity 
vacuum mode (see Fig. 5). This combination gives rise 
to friction forces along all three spatial directions. The 
achieved cooling efficiency led to microkelvin tempera- 
tures and to an average single-atom trapping time in the 
high-finesse cavity as long as 17 seconds, during which 
the strongly coupled atom could be observed continu- 
ously, see Fig. 6. 
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FIG. 6 Demonstration of cavity cooling and long-time trap- 
ping of a controlled number of atoms inside the cavity, (a) 
Shown is a single trace of the recorded photon-count rate in- 
dicating the capture of an atom 75 ms after switching on the 
pump laser (see Fig. 5). Within 100 ^s, the scattering rate 
reaches a steady-state value, (b) The recorded photon-count 
rate allows for determining the atom number and the trap- 
ping time, (c) The analysis of 50 traces, each of 6 s duration 
and starting with one atom, yields an average lifetime r of 
17s (upper curve), whereas single atoms that are not exposed 
to the pump laser only reside for 2.7 s in the cavity volume 
(lower curve). From Nussmann et al. (2005). 



Further studies revealed that the previous setting of 
detunings and intensities can be extended to more gen- 
eral geometries, including the motion perpendicular to 
the cavity axis, or the external driving of the atom in- 
stead of the cavity. What is required for cooling is an in- 
homogeneity in the system on the wavelength scale which 
leads to a position-dependent steady-state of the coupled 
atom-cavity system. This inhomogeneity can arise from 
the cavity mode function, as for the result in Eq. (14), 
but also from a standing-wave pump field, or from the 
spatially modulated AC-Stark shift in a strong standing- 
wave laser field (Murr et al., 2006b). All these sources 



In the experiment (Fig. 5), a far-detuned standing- 
wave dipole trap which is oriented perpendicular to the 
cavity axis is used to transport atoms into the cavity 
(Dotsenko et al, 2005; Kuhr et al, 2003). The combi- 
nation of controlled insertion of single atoms into and 
retrieval out of a high-finesse optical resonator with cav- 
ity cooling led to a deterministic strategy for assembling 
a permanently bound and strongly coupled atom-cavity 
system. Long storage times well above 10 seconds and the 
controlled positioning of single or a given small number 
of atoms on the submicrometer scale are simultaneously 
available (Khudaverdyan et al, 2008; Nussmann et al. 
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2005). 



The exploration of cavity cooling was a stimulating 
and essential step for the experimental achievement of 
strongly-coupled cavity QED systems combined with the 
control over the atomic motion. With the implementa- 
tion of free-space laser cooling and trapping techniques 
in cavity experiments sufficiently long atom-cavity inter- 
action times were demonstrated. These achievements led 
to remarkable experimental breakthroughs and applica- 
tions in single-atom cavity QED, recently. For exam- 
ple, high-precision measurements demonstrated the ba- 
sic cavity QED model in the optical domain, i.e. , by re- 
solving the doublet of the lowest-lying excitations of the 
atom-cavity system Boca et al. (2004) and Maunz et al. 
(2005), as well as the quantum anharmonic domain of 
the Jaynes-Cummings spectrum (Kubanek et al., 2008; 
Schuster et al., 2008) where squeezed light can be readily 
generated (Ourjounitsev et al., 2011). Furthermore, the 
achieved trapping times permitted the development of 
a deterministic single-photon source (Kuhn et al., 2002; 
McKecvcr et al., 2004), for having full polarization con- 
trol (Wilk et al., 2007b), and to realize the long time 
sought atom-photon quantum interface (Boozer et al., 
2007; Wilk et al., 2007a) and single-atom quantum mem- 
ory (Spccht et al., 2011). Many prosperous directions can 
grow out from the realization of the elementary case of 
electromagnetically induced transparency with a single 
atom (Figueroa et al., 2011; Kampschulte et al., 2010), 
such as for example all-optical switching with single pho- 
tons. 

3. Temperature limit 

The extra friction term induced by the cavity is closely 
connected to the modification of the zero-point field fluc- 
tuations. Indeed, even for large atom-field detuning the 
diffusion within a cavity-sustained optical dipole trap can 
be an order of magnitude larger than for a free-space 
field (van Enk et al, 2001; Murr et al, 2006a; Puppe 
et al, 2007b). The heating rate due to fluctuations can 
be explicitly calculated in a semiclassical approach (see 
Sec. II. A. 3) which allows to estimate the stationary tem- 
perature attained by the atoms. Under optimal condi- 
tions one finds the intuitive result 

ksT « hn , (16) 

which is independent of the atomic parameters. This re- 
sult was confirmed by numerical simulations (Domokos 
et al, 2001) and fits very well to experimental obser- 
vations. Interestingly, the result remains largely valid 
in the limit where the temperature reaches the recoil 
limit ksT K hujn which is not governed by the semi- 
classical description anymore. For a particle trapped in 



a harmonic potential with vibrational frequency v > k 
(resolved-sideband regime) efficient ground state cool- 
ing was proposed by Zippilli and Morigi (2005). Quan- 
tum interference effects in the spontaneous emission of a 
trapped particle in a cavity allow for ground-state cooling 
even in the bad-cavity regime where v < k (Cirac et al, 
1995). This prediction does not contradict Eq. (16), since 
it was made for a strongly localized trapped atom (Lamb- 
Dicke regime) at a precisely given position, whereas the 
temperature limit above assumes spatial averaging over 
the cavity wavelength. 

According to Eq. (16), there seems to be no lower 
bound on the temperature as long as the cavity finesse 
can be increased. However, with decreasing loss rate, 
the capture range of the cavity cooling mechanism also 
shrinks. This relation is thoroughly discussed by Murr 
(2006), based on an explicit expression for the friction 
force obtained for arbitrary velocity. When applying very 
strong cavity fields, a 3> 1, in close analogy to the mean 
field treatment of optomechanical models (Genes et al, 
2008), it is possible to effectively enhance the weak atom- 
field interaction appearing at very large detunings to an 
effective strong coupling gcs = goct at the expense of 
introducing extra fluctuation terms (Nimmrichter et al, 
2010). This setting can considerably speed up the cooling 
process and enhance the capture range, while still leading 
to a similar final temperature as given by Eq. (16). 

4. Cooling in multimode cavities 

The atom-field dynamics qualitatively changes when 
invoking several cavity modes to participate as dynam- 
ical degrees of freedom. In simple terms, not only the 
magnitude but also the spatial shape of the optical poten- 
tial and the associated light forces become a dynamical 
quantity. 

This can be easily demonstrated at the generic ex- 
ample of a ring cavity geometry (Gangl and Ritsch, 
2000). In the regime of dispersive atom-field coupling, 
the atom not only modifies the resonance frequencies of 
the two counter-propagating cavity modes, thereby tun- 
ing their field amplitudes, but also gives rise to phase 
locking by coherent photon redistribution between the 
cavity modes. This determines the position of the nodes 
and antinodes of the emergent standing- wave interference 
pattern of the cavity radiation field. For a red-detuned 
pump field, Aa < 0, the particle is drawn to an antinode 
of the field which, at the same time, gets dragged along 
with the slowly moving atom (assuming kv < ^,7). Due 
to the delayed response of the intracavity field, however, 
the particle is permanently running uphill and thus ex- 
periences a friction force. The two-mode geometry of a 
ring cavity has been shown to result in faster cooling and 
larger velocity capture range as compared to a single- 
mode standing-wave cavity (Gangl et al, 2000; Schulze 
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et al., 2010). Moreover, the laser pump configuration 
used for polarization gradient cooling or velocity-selective 
coherent population trapping can be envisaged within a 
ring cavity, for which case very efficient cavity cooling is 
predicted without fundamental lower limit on the tem- 
perature (Gangl and Ritsch, 2001) 

The more modes in a cavity are available in the vicin- 
ity of the pump frequency, the smaller is the trans- 
verse length scale on which the field shape gets mod- 
ulated in the presence of an atom. On the one hand, 
this leads to stronger three-dimensional localization of 
atoms around their self-generated intensity maximum 
(Salzburger et ai, 2002). On the other hand, the cool- 
ing time reduces in a highly-degenerate confocal cavity 
more or less quadratically, whereas the diffusion increases 
only linearly with the effective number of modes involved 
(Domokos et al., 2002; Nimmrichter et al, 2010). 

The scope of cavity-mediated optical manipulation of 
atoms significantly enlarges also in the case of many-atom 
systems, which we will review in Sect. IV. D. 3. 

C. Extensions of cavity cooling 

The general principle of cavity cooling is expected to 
be applicable in a broad range of other systems with dif- 
ferent radiation field geometries or other material com- 
ponents. 



1. Cooling trapped atoms and ions 

There are several experimental systems in which 
trapped atoms are strongly coupled to a high-finesse cav- 
ity. Ion trap setups have been combined with high- 
finesse cavities in the moderate coupling regime (Her- 
skind et al, 2009; Keher et al., 2004). There are all- 
optical schemes, too, where different longitudinal modes 
of a standing-wave cavity are used to separate the op- 
tical trap modes from the cooling ones (Maunz et al., 
2004; Schleier-Smith et ai, 2011). State-insensitive cool- 
ing and trapping of single atoms employing light field 
at magic wavelengths, which induces an almost identi- 
cal AC-Stark shift of the two relevant electronic states, 
has been demonstrated (McKeevcr et al, 2003). Further, 
trapping of atoms in low field regions of a blue-detuned 
intracavity dipole potential has been investigated exper- 
imentally (Puppe et al., 2007a). In a similar intracavity 
dipole trap, the axial atomic motion was cooled down to 
the ground state by way of coherent Raman transitions 
on the red vibrational sideband, meanwhile the atomic 
motion was inferred from the recorded Raman spectrum 
by Boozer et al. (2006). 

The cavity cooling mechanism operates also in the case 
of tightly confined particles. In the Lamb-Dicke regime 
for tightly confined particles, ^Ju^jijv <^ 1 where v de- 



notes the harmonic trap frequency, explicit expressions 
for the cooling and heating rates of the CM motion of an 
atom, trapped in an optical resonator and driven by a 
laser field, have been derived both in the regime of weak 
and strong atom-cavity coupling. In the former, a variant 
of sideband cooling appears (Cirac et at, 1995; Vuletic 
et ai, 2001). Experimentally, the cavity cooling of a sin- 
gle trapped ^^Sr"*" ion in the resolved-sideband regime 
has been demonstrated and quantitatively characterized 
recently (Leibrandt et ai, 2009). The spectrum of cav- 
ity transmission, the heating and cooling rates, and the 
steady-state cooling limit have been measured in perfect 
agreement with a rate equation theory. The final tem- 
perature corresponding to 22.5(3) occupied vibrational 
quanta was limited by the moderate coupling between 
the ion and the cavity. 

The calculations have been extended to the strong- 
coupling regime, where higher-order transitions between 
eigenstates of the coupled system have been identified 
and novel non-trivial parameter regimes leading to cool- 
ing have been revealed (Blake et ai, 2011; Zippilli and 
Morigi, 2005). In the resolved-sideband regime, ^ k, 7, 
the discreteness of the vibrational spectrum, which is the 
same for the electronic ground and excited states, gives 
rise to interference between different transition paths in 
analogy to the cooling of trapped multilevel atoms (Mo- 
rigi et at, 2000). Ground state cooling is achievable ac- 
cording to the theoretical predictions (Zippilli and Mo- 
rigi, 2007). 



2. Cooling nanoparticles and relation to optomechanics 

The fact that cavity cooling requires only linear po- 
larizability suggests that it could be directly applicable 
to large objects, such as nanobeads (Barker and Shnei- 
der, 2010; Chang et ai, 2009), thin reflective membranes 
(Genes et ai, 2009), or even small biological objects such 
as viruses (Romero-Isart et al, 2010). Moreover, since 
membranes, being macroscopic objects, can have large 
static polarizability (refractive index) , the cooling can be 
much more efficient than for single atoms or molecules. 
Indeed, there is a strong connection between cavity cool- 
ing of atoms and dispersive cavity optomechanics (Jayich 
et ai, 2008; Thompson et ai, 2008), which can easily be 
seized in the framework of the scattering models (Xuereb 
et ai, 2009). Cavity cooling of membranes experimen- 
tally shows great success down to the vibrational quan- 
tum ground state (Jayich et al., 2011). 

As the local field strength is strongly enhanced inside 
a resonator, optical dipole traps can be operated at very 
large dctunings, where only the static polarizability of 
the particle is relevant (Deachapunya et ai, 2008; Nimm- 
richter et ai, 2010). In such a setting of coupled optical 
and mechanical systems, the ring cavity with degenerate 
pairs of counterpropagating modes, or other configura- 
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tions where degenerate modes are available, can offer the 
realization of various effective models. 

Consider, for example, a symmetrically pumped ring 
cavity. The field can be written as a superposition of the 
strongly pumped and thus highly excited cosine mode 
and the empty sine mode. The cosine mode fulfills two 
purposes: (i) it generates the trapping potential, and (ii) 
it feeds the sine mode through photon scattering off the 
particle (atom, molecule, membrane). The model Hamil- 
tonian is of the form (Schulze et ai, 2010): 



H 



P 

2m 



at a.. 



'hU{x)-\-ih (rjal — ri*ac 



(17) 

where U (x) is the dispersive interaction potential, and Uc 
(fls) denote the field amplitudes of the cosine (sine) mode, 
respectively. Linearizing the position around the trap 
minimum, we can recover the standard optomechanical 
Hamiltonian, 



H 



2m 



-m2hUoa\ac{kx) 



- ?i(Ac - UQ)a\ac ~ hAcalas 



- hUl^ia, + al)x , (18) 

with quadratic coupling to the cosine trapping mode and 
linear coupling to the sine cooling mode. As the particle 
couples the two modes there appears an energy splitting 
which allows to extract via inelastic scattering kinetic en- 
ergy from the vibrational motion in the optical trap (El- 
sasser et ai, 2003). As for standard cavity cooling the 
final temperature in the classical regime is again limited 
by the cavity linewidth fc^T « Hk. However, in a very 
good cavity, when the pump field is sufRciently strong, 
one can reach the resolved-sideband regime, where the 
trap frequency v exceeds the cavity linewidth, and the 
final temperature would correspond to less than a sin- 
gle excitation fcsT < hv. In this ground-state cooling 
limit one has to resort to a quantum description of mo- 
tion and the optical fields. Interestingly, the sine mode 
automatically acts as a built-in monitoring system which 
continuously observes the vibrational quantum state of 
the particle in the cosine mode. Hence close to T = 
one can observe quantum jumps of the particle via the 
sine mode photon counts (Schulze et ai, 2010). 



3. Cooling molecules 

a. Cooling the translational motion of molecules Molecular 
structure fundamentally alters and complicates the pic- 
ture conceived for laser-cooling two-level atoms. Upon 
excitation from the pump field, the molecule can relax 
by either Rayleigh scattering back to the ground state 
\g) at the rate 7Ry or Raman scattering to metastable 
states at the rate 7Rn- There is a multitude of metastable 



molecular states (spin-orbit, rotational, and vibrational) 
available via inelastic Raman scattering. The generally 
low free-space branching ratio 7Ry/7Rn results in popu- 
lation shelving after only a few photon scattering events, 
thereby prematurely quenching the cooling process. Be- 
cause of the prohibitive expense of building multiple re- 
pumping laser systems, optical cooling of molecules via 
free-space dissipative scattering of photons is thought not 
to be practicable. 

Since cavity-assisted laser-cooling relies on the cavity 
dissipation channel, it has been suggested as a potential 
method to mitigate Raman loss. Spontaneous photon 
scattering, in principle, can be entirely suppressed by us- 
ing large detuning. However, as discussed in Sec. II. B. 2, 
in order to keep the cooling efficiency constant, one needs 
to preserve a given level of excitation in the atom or 
molecule. Therefore, merely the large detuning does not 
solve the branching ratio problem of molecules (Lev et al., 
2008). To overcome this severe problem, the use of an 
optical cavity with cooperativity parameter much larger 
than unity is mandatory, in accordance with Eq. (15). In 
this case the enhanced coherent Rayleigh scattering into 
a decaying cavity mode can ensure a vanishingly small 
probability of the molecule to Raman scatter during the 
cooling time. For CN diatomic molecules, the cavity cool- 
ing process has been calculated numerically (Lu et ai, 
2007). 



b. Cooling the rotation and vibration of molecules While 
theoretical models and experiments have mostly concen- 
trated on the center-of-mass motion of structureless po- 
larizable particles or two-level atoms, the complex rovi- 
brational structure of molecules is one of the central ob- 
stacles preventing efficient laser-cooling of molecules. In 
many common beam sources the initial temperature can 
be designed to be low enough to freeze most vibrations 
and only leave few rotational quanta (Rangwala et ai, 
2003). Nevertheless, the interaction with the cooling 
laser light will in general start to redistribute the popula- 
tion within the rovibrational manifolds strongly altering 
the optical properties of the molecules and hampering 
further cooling. Only a few exceptions of this rule have 
been discovered and investigated lately (Shuman et ai, 
2010). Cavity cooling, however, can in principle be de- 
signed to counteract this heating process and even fur- 
ther cool the rovibrational energy of molecule. As an 
enourmous spread of transition frequencies is required to 
facilitate this, it proves advantageous to simultaneously 
apply a multitude of different longitudinal cavity modes 
(Kowalewski et ai, 2007; Morigi et al, 2007). Simula- 
tions show that the rovibrational cooling can be com- 
bined with motional cooling, e.g., in a trap (Kowalewski 
et ai, 2011), to get a cold molecular gas in all degrees of 
freedom. At this point a practical implementation would 
require precooling by other methods, such as the opto- 
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electrical scheme proposed by Zeppenfeld et al. (2009), to 
achieve sufficient interaction times and densities within 
the cavity mode volume. 



4. Cooling and lasing 

Collective coherent emission of a laser-driven atomic 
ensemble into the field of an optical cavity accompa- 
nied by a very fast and efficient cooling of atomic mo- 
tion was observed in an experiment conducted by Chan 
et al. (2003). Although the effect has not yet been fully 
understood, it is attributed to Raman gain within a Zee- 
man manifold. The combination of cavity cooling with 
intracavity gain is an intriguing prospect. It was initially 
suggested by Vuletic (2001) to transform a bad cavity ef- 
fectively into a good cavity with fast cooling towards an 
even lower temperature. While the principle idea proves 
to be correct, a more realistic and detailed modeling, 
which accounts for fluctuations to consistently treat the 
gain, gives a higher limit of the achievable temperature 
(Salzburger and Ritsch, 2006). This observation was also 
confirmed in the optomechanical regime of cavity cooling, 
where intracavity gain leads to faster cooling but a higher 
final temperature (Genes et al., 2009). 

In a standard setup, the intracavity gain could be gen- 
erated by an additional inverted medium placed within 
the cavity. This would lead to a technically challeng- 
ing setup, if one aims to operate in the strong-coupling 
regime. Interestingly, it turns out that in a conceptually 
much simpler configuration, the gain can also be provided 
by the same atomic medium which is aimed to be cooled 
in the setup. Of course, such a scheme requires a suit- 
able pumping mechanism which transfers atoms from the 
lower to the upper level of the cooling transition, with- 
out introducing too much extra noise. In the ultimate 
limit one can envisage a single atom, which is exter- 
nally pumped within a high-finesse cavity. Stimulated 
emission into the cavity mode provides gain to create 
a trapping potential for the atom. For a blue-detuned 
cavity this gain simultaneously extracts motional energy 
from the particle and thus provides cooling (Salzburger 
and Ritsch, 2004). Fortunately, an inverted atom is a 
high-field seeker in the blue-detuned light field, so that 
it will be trapped close to optimal gain. Hence this 
setup provides for lasing, trapping and cooling of a sin- 
gle atom within a resonator forming the most minimal- 
istic implementation of a laser (Salzburger et al., 2005). 
The system can be generalized to several particles, which 
strongly reduces the requirements on the pump mecha- 
nism (Salzburger and Ritsch, 2006). In the limit of ul- 
tracold gases in an optical lattice, stimulated optical gain 
occurs concurrent with Bose enhanced coherent popula- 
tion of the lowest energy band. While for a pulsed setup 
this constitutes in principle a very fast and efficient cool- 
ing method, a CW setup could provide a possible route 



towards the realization of a CW atom laser (Salzburger 
and Ritsch, 2007, 2008). 

5. Monitoring and feedback control 

Starting from the early days of cavity QED, a strongly 
coupled atom-cavity system was considered as a number- 
resolving neutral particle detector (Mabuchi et al., 1996), 
a concept which is still being developed and implemented 
in miniaturized devices (Tcper et al., 2006). Going one 
step further, the high-finesse resonator acts as a micro- 
scope with which the trajectory of individual atoms can 
be reconstructed from the recorded cavity transmission 
with high spatial (< ^m) and temporal (< /is) resolution 
(Hood et al., 2000). The method can be considerably 
improved by the use of multimode cavities. The particle 
does not only modify the phase and intensity of the intra- 
cavity field, but redistributes light between the different 
spatial modes. The output field imaged on a CCD cam- 
era therefore allows to monitor directly and in real time 
the motion of the particle (Horak et al., 2002; Maunz 
et al., 2003). Note that even for incomplete position in- 
formation at any given time, the most likely trajectory 
of single atoms can be reconstructed with the help of 
inversion algorithms based on the coupled equations of 
motion. 

Once the position and motion of the particle are 
known, it is straightforward to apply feedback on the 
motion of a single atom by adjusting the pump lasers to 
steer the particle motion within the cavity and increase 
its trapping time (Fischer et al., 2002). The cavity field 
both provides particle detection and mediates the feed- 
back force. This method has been successfully refined 
by several groups and resulted in an increase of single- 
particle trapping times by several orders of magnitude 
(Kubanek et al., 2009, 2011). By applying controlled 
and delayed feedback forces on the particle, its kinetic 
energy can be reduced as well. This kind of feedback 
cooling resembles stochastic cooling techniques applied in 
high-energy physics. Strongly enhanced cooling has been 
predicted when the feedback scheme, consisting of time- 
dependent switching of the trapping field as a function 
of the intracavity intensity, is operated in the dispersive 
bistability regime (Vilensky et al., 2007). This method 
should also give new prospects to optomechanical setups. 
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III. COLD ATOMIC ENSEMBLES IN A CAVITY 

New research directions opened in cavity QED when 
cold and uhracold atomic ensembles were successfully 
prepared within high-finesse optical resonators. In the 
many-body configuration, the common coupling of atoms 
to the cavity field creates a wealth of new possibilities 
to implement tailored atom-atom interactions over large 
distances, an ingredient which usually is absent in free- 
space cold atom experiments. 

The atom-atom coupling is mediated by the cavity ra- 
diation field between the AC electric dipole moments. 
However, its nature is inherently different from the free- 
space dipole-dipole interaction. In a cavity, the inter- 
action strength does not decay with the interatomic dis- 
tance and depends only on the local coupling of the atoms 
to the cavity field. Fundamentally, the interaction is not 
binary: the ensemble of atoms collectively acts onto the 
state of the radiation field which then reacts back on the 
individual atoms. This scenario is generally referred to as 
global coupling. The range of the interaction is given by 
the size of the cavity mode, which can be macroscopic. In 
cases where single-atom strong coupling is not achieved, 
the collective energy exchange still can be dominated by 
coherent interaction. 

After discussing the nature of the long range atom- 
atom interaction mediated by a cavity field in various 
geometries, we consider first the many-body influence 
on the cavity cooling scheme. Then we address the 
most spectacular collective effects realized by cold atoms 
within linear and ring cavities. Critical phenomena, in- 
stability thresholds, and scaling laws will be discussed 
by means of various mean field theories in the end of this 
section. 



A. Collective coupling to the cavity mode 

Resonant coherent coupling between an ensemble of N 
two-level atoms and a single standing-wave cavity mode 
is described by the many-body generalization of Eq. (1) 

H/h = -Aca^a - ^ AA(rj)cr]crj- 
j 

+ J2igfir,)ia]a~a^a,) (19) 



where j = 1 . . . N labels the atoms, and the mode func- 
tion /(r), for simplicity, is real. The atomic ensemble can 
be represented by a single collective dipole with effective 
coupling strength only if (i) the atomic motion can be av- 
eraged out (ii) only the cavity mode is laser-driven, and 
(iii) the atoms are in the low saturation regime. In this 
case the atoms collectively couple to the cavity mode with 
an effective strength of ^off = 9\jY^j Pi'^j) (the sum- 
mation index runs from 1 to N). Correspondingly, an 



A^-fold enhancement appears for the many-atom system 
in terms of the single-atom cooperativity, C = g'^ /{2k-j), 
which measures the ratio of light scattering into cavity 
mode versus surrounding vacuum modes (Tuchman et al., 
2006). For example, the strong distortions of the single- 
atom normal mode splitting in the cavity transmission 
spectrum induced by a thermal beam of atoms crossing 
the cavity could be interpreted by such a collective mode 
picture (Raizen et al., 1989). Employing an optical con- 
veyor belt, an adjustable number, N = 1 . . . 100, of cold 
atoms has been transported into a microcavity, and large 
nonlinearities have been achieved as evidenced by the ob- 
servation of absorptive optical bistability in the real-time 
transmission spectrum (Saucr et al., 2004). 

In general, however, one has to consider the many- 
body system composed of a large number of internal and 
motional degrees of freedom. We will exhibit this in the 
following at the simplest nontrivial case of two atoms in 
the same mode. 



1. Cavity-mediated atom-atom interaction 

Let us discuss the character of the cavity-mediated 
atom-atom interaction in two different pump geometries, 
namely, pumping the cavity field either directly or indi- 
rectly via light scattering off the laser-driven atoms. 



a. Cavity pumping Consider iV atoms moving in the field 
of a laser-driven optical cavity. The detuning between the 
driving laser and the dispersively shifted cavity resonance 
depends on the position of all atoms, which in turn expe- 
rience the optical dipole force of the intracavity field. For 
small atomic velocities and in the low saturation limit, an 
adiabatic potential can be deduced (Fischer et al, 2001) 



T^(ri,...r^) 



-atan-^'^'^-^^^+^'ff 



Aak + Ac7 



Aak + Ac7 



(20) 

which is analogous to the Born-Oppenheimer approxima- 
tion used for describing the motion of nuclei in molecules 
in the averaged electronic potential. The potential V de- 
pends on the atomic positions solely via the collective 
coupling strength ges, and thus is valid for any num- 
ber of atoms. This is not surprising as the adiabatic 
force is calculated by freezing the atomic motion. The 
resulting cavity-mediated long-range atom-atom inter- 
action gives rise to an asymmetric deformation of the 
normal-mode splitting as was observed experimentally by 
Miinstermann et al. (2000). 

In the case of two atoms with positions xi and X2 the 
interaction potential landscape V{xi, X2) along the cavity 
axis is shown in Fig. 7 for two different parameter set- 
tings. The upper graph corresponds to the experimental 
parameters used in the Garching group (Miinstermann 
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FIG. 7 The adiabatic cavity potential as a function of the 
atomic positions xi and X2- The cavity is quasi-resonantly 
excited by a pump laser with detuning Ac = — k + Uo. The 
detuning from the atomic resonance is set to A a ~ —507 
to ensure the suppression of spontaneous photon scattering. 
In the upper graph, typical experimental cavity parameters 
(ft — 7/2, g = 57) (Miinstcrmarm et al, 2000) have been 
used, whereas in the lower graph g was increased fourfold. In 
the first case, the potential is well approximated by a sum of 
two single-particle potentials. In the second case, either both 
atoms are trapped or free. From Asboth et al. (2004). 



b. Atom pumping The situation drastically changes if 
the atoms are laser-driven from a direction perpendic- 
ular to the cavity axis. Intracavity photons are then cre- 
ated by Rayleigh scattering of laser photons into the cav- 
ity mode. Due to hght interference, the scattered intra- 
cavity field exhibits a very sensitive dependence on the 
interatomic distance. For two atoms separated by odd 
integer multiples of the half- wavelength, the correspond- 
ing scattering amplitudes into the mode have the same 
magnitude but opposite sign, resulting in destructive in- 
terference and a vanishing cavity field amplitude. On the 
other hand, for atoms separated by even integer multiples 
of the half-wavelength, the field components scattered off 
the two atoms interfere constructively. Compared to the 
field intensity created by a single scatterer, the latter case 
yields a fourfold enhancement of the intensity, referred to 
as superradiance (DeVoe and Brewer, 1996; Dicke, 1954). 

In the case of pumping directly the atoms the force 
along the cavity axis acting on the individual atoms due 
to light scattering cannot be expressed as a gradient of 
a collective potential, at variance to Eq. (20). One can 
admit this by checking that V^Fj 7^ VjF^, where V; is 
the gradient with respect to the coordinate r^, and Fj 
is the force acting on atom j. If there was a potential 



V such that F, 



-V,y, the two sides should be equal 



as they are the second derivatives of the potential and 
the order of taking the derivatives is irrelevant according 
to Young's theorem. The fact that the force can not be 
derived from a potential is not so surprising, in hindsight, 
as we are dealing with an open system with continuous 
energy exchange with the environment and an unlimited 
energy resource in the form of the pump laser. Actually, 
the existence of a potential Eq. (20) for the cavity-driving 
geometry is the exceptional case. 

Approximately, in the limit of Uq^Tq 0, more pre- 
cisely N'^Uo {k, I Ac I), the motion of the atoms is gov- 
erned by the collective potential 



et at, 2000). Although the single-atom light shift is com- 
parable with the cavity linewidth |J7o| ~ k, the effec- 
tive interaction between the atoms is relatively weak and 
the potential resembles the familiar "egg-carton" surface 
proportional to sin^(A:a;i) -I- sin^(A;x2). For 'artificially' 
enlarged atom-field coupling g = 2O7, shown in Fig. 7 
(lower graph), the atom- atom interaction strongly affects 
the potential landscape felt by the second atom, depend- 
ing on the position of the first atom and vice versa. For 
this parameter setting the single-atom light shift is suf- 
ficiently large, Uq 3> k, that removing one atom from 
the cavity antinode makes the potential experienced by 
the other atom vanish. Note that the trap is deeper for 
the smaller coupling of the upper graph. Interestingly, 
as was shown by Asboth et al. (2004), the motion of the 
two atoms gets correlated even for the parameter setting 
of the upper graph as a consequence of additional non- 
conservative forces (see Sec. III. A. 2). 



V{ri, . . .Tat) = ?i 



cos{kxj) cos(fczj) 



(21) 

where cosine mode functions were assumed for the cav- 
ity and the pump laser field. The interference effect is 
manifest: when scanning the atom-atom distance over a 
wavelength, the contrast of the interference in the cavity 
field intensity is unity regardless the atom-cavity cou- 
pling constant g. This is not the case for cavity pump- 
ing, where, in the limit of small coupling constant g, the 
atoms cause only a small modulation of the cavity in- 
tensity. Therefore, the atom pumping geometry lends 
itself to observe spectacular many-body effects even in 
the weak-coupling regime. 

The superradiant light scattering into the cavity is the 
basis of various collective dynamical effects, which have 
been more profoundly studied theoretically. Although we 



17 



photon number atomic populations 




X2 / A X2 / A 

FIG. 8 Collective scattering of two atoms into the cavity 
mode including internal atomic excitation, (a) Mean number 
of photons scattered into the cavity mode at pump-cavity res- 
onance, Ac ~ 0, with coupling strength gcos{kx). The atoms 
are placed at a;i = and X2, respectively, (b) Excited state 
populations Hi and 112 of the two atoms (dashed and solid 
lines). The parameters are k = O.27, 77cff = 'y,g ~ IO7, Aa = 
IOO7. From Zippilh et al. (2004b). 

mostly neglect atomic saturation effects in this review, it 
is important to reveal modifications of the interference 
effect in the collective scattering when a small but finite 
atomic saturation is taken into account. Since the satura- 
tion also depends on the relative distance of the particles, 
new types of nonlinear behavior take place. For exam- 
ple, as shown in Fig. 8, the destructive interference for a 
separation of half-wavelength between the atoms is not 
perfect any more and the photon scattering generates a 
nonclassical cavity field with zero amplitude but finite 
photon number (Vidal et a/., 2007; Zippilh et al, 2004a). 

2. Collective cooling, scaling laws 

As was discussed in Sec. II.B.I, the cavity cooling force 
on single atoms stems from a delicate correlation between 
the atomic motion and the retarded dynamics of the cav- 
ity field. In a many-atom system it is at first unclear 
what happens to these correlations in the presence of 
other moving atoms. Furthermore, the cavity-mediated 
crosstalk between atoms has a component sensitive to the 
atomic velocities (Domokos and Ritsch, 2003), i.e., atom 
1 moving at velocity vi induces a linear friction force on 
atom 2, which might yield correlations in velocity space. 
To answer this question, one can straightforwardly gen- 
eralize the semiclassical model, presented in Sec. II. A. 3, 
for many atoms. In general, however, this leads to an 
analytically intractable problem. The dynamics of the 
many-atom system cannot be reduced to that of an ef- 
fective mode, as was the case for the adiabatic poten- 
tial Eq. (20) for atoms at rest. The two-atom case has 
been discussed in detail by Asboth et al. (2004), who 
found, using the parameter regime of Fig. 7a, a buildup 
of strong correlations in the motion of two atoms due to 
the velocity-dependent cavity forces. 

The scaling of the cavity cooling efficiency with the 



number of particles has been studied by means of numer- 
ical simulations for = 1 . . . 100 in the limit of a weakly 
driven single-mode field, where the optical dipole poten- 
tial negligibly perturbs the free motion of atoms along 
the cavity axis (Horak and Ritsch, 2001c). If the param- 
eter Uq is chosen sufficiently small so that the collective 
light shift is still below the cavity linewidth, NUo < k, 
the rate of kinetic energy dissipation is independent of 
the number of atoms. This suggests that the individual 
atoms in the cloud are cooled independently from each 
other, although they are all coupled to the same cavity 
mode. This holds only in the weak-coupling limit which 
is not practical for cooling since the cooling time is long. 

When the collective coupling to the cavity mode is sig- 
nificant with respect to the linewidth n, the scaling be- 
havior of cooling with the number of atoms has been 
studied by keeping NUq and ri/\/N constant while vary- 
ing the atom number A^. The former ensures an identi- 
cal maximum collective light shift induced by the atoms, 
the latter amounts to a nearly constant optical potential 
depth (proportional to Uo'tj^/k'^)- With this rescaling 
of the parameters, the effect of individual atoms on the 
cavity field diminishes as the number of atoms increases. 
The final temperature was found invariant, however, the 
cooling time increases linearly with the atom number N. 
As long as the driving rj is weak enough to result in a 
shallow optical potential depth, in which the atoms move 
almost freely, all the motional degrees of freedom along 
the cavity axis are cooled. 

In the limit of tightly confined atoms, both theoreti- 
cal calculations (Asboth et al., 2004; Nagy et al., 2006a) 
and experiments (Schleier-Smith et al., 2011) proved that 
only the center-of-mass motion is damped by the cavity- 
induced friction force (Gangl and Ritsch, 1999, 2000). An 
efficient sideband cooling scheme has been proposed by 
Elsasscr et al. (2003) for particles confined in the optical 
lattice potential generated by two counter-propagating 
degenerate modes of a ring cavity. The scheme relies on 
the collective atom-field coupling which lifts the degen- 
eracy and creates two standing-wave modes phase-locked 
by the back-scattering of light. The lower-lying mode sus- 
tains the optical lattice with an intensity adjusted such 
that the upper-lying mode becomes resonant with the 
vibrational anti-Stokes Raman transition. The sideband 
cooling allows to reach the vibrational ground state. 

A collective enhancement of friction on the center-of- 
mass motion has been demonstrated experimentally in 
the transverse pump configuration, as an accompanying 
effect of the self-organization into a Bragg-scattering lat- 
tice (see Sec. III.B.I). Peak decelerations of — lO'^m/s^ 
have been observed, and the damping effect has been 
demonstrated with light-atom detunings up to Aa/'^tt — 
-6 GHz. For smaller detunings AA/27r « -160MHz, 
similarly large velocity-dependent friction forces (up to 
decelerations of —1500 m/s^ and temperatures as low as 
7/iK) have been observed in another set of experiments 
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FIG. 9 Observation of collective friction force on the center- 
of-mass motion. The cavity output power (thin line, left 
scale) is shown during illumination of a freely falling atomic 
cloud (initial velocity 15 cm/s) with a transverse laser beam. 
The initial increase signifies the self-ordering into a Bragg- 
scattering lattice (see Sec. III.B.l). Deceleration of the center- 
of-mass motion is recorded via the beat signal recorded over 
300/is. The modulation stems from the spatial variation of the 
atom-cavity coupling with period A/2, whereby atoms can- 
not scatter at a node of the intracavity standing wave. The 
changing modulation frequency indicates the atomic deceler- 
ation (thick line, right scale). The inset shows the density 
profile (a.u.) of the atomic cloud after free expansion with- 
out (gray line) and with (black line) a 400 fis exposure to the 
pump beam. A fraction of about one third of the atoms is 
delayed significantly in accordance with the measured decel- 
eration. Here the pump I/Is = 420 (Is = l.lmW/cm^ is the 
D line saturation intensity of Cs), AA/2n — —1.58 GHz, and 
Ac/2tv = -10 MHz. The atom number is iV = 2.6 x lO'^. 
From Black et al. (2003). 

(Chan et al., 2003). Both the large friction and the low 
temperature cannot be explained in terms of the interac- 
tion between single atoms and the cavity field. In this de- 
tuning regime a theoretical description is more involved 
because the entire hyperfine manifold has to be taken 
into account and the interaction can lead to Raman las- 
ing between different magnetic sublevels. 



3. Back-action, nonlinear dynamics 

In general, the interplay between the mechanical ef- 
fect of light on the atomic motion and light scattering 
inside the cavity off the spatial atomic density distribu- 
tion can lead to highly non-linear dynamics. Having a 
large atomic ensemble organized in a lattice structure, 
for example, the collective Bragg scattering much more 
efficiently redistributes the light between modes than 
Rayleigh scattering from the individual atoms. The en- 
hancement factor, being on the order of the number of 
atoms, can give rise to significant sensitivity of light scat- 
tering to small variations of the spatial distribution. As 
an example, the back scattering between the counter- 
propagating modes (denoted by -|- and — ) of a ring cavity 
was found to depend strongly on the bunching parameter 



of the atomic distribution around the trapping sites of an 
optical lattice. In experiments performed by the Ham- 
burg group (Elsasser et al, 2004; Nagorny et al, 2003), 
the amplitude a+ in one of the modes was actively sta- 
bilized by a feedback loop, d+ = 0, and thus the other 
mode obeyed the nonlinear equation of motion 

d_ = iNUo B Hia- - iNUoa+ B* + ?/_ . (22) 

a+ 

Here, rj^ denotes the driving amplitude of this 
mode and B — (e~^''^^) the bunching parame- 
ter (see also Sec. III.B.2). For a thermal cloud 
the atomic bunching follows approximately B oc 
al/|a_| exp |— const/-\/|a_||. The resulting nonlinear 
dynamics was exemplified by a new kind of optical bista- 
bility in the dispersive atom-field coupling regime, which 
is outside the range of optical bistability effects relying 
on the nonlinearity of the internal atom-field coupling 
(Lugiato, 1984). Subsequent experiments (Klinncr et al, 
2006) also revealed the mechanical effect of light on the 
atomic distribution in the dispersive regime through the 
normal-mode splitting. 

B. Non-equilibrium pliase transitions and collective 
instabilities 

The nonlinear collective dynamics of thermal atoms in 
a high-finesse resonator can give rise to non-equilibrium 
phase transitions and collective instabilities. In the fol- 
lowing we present two, experimentally evidenced exam- 
ples which have been theoretically studied both in the 
thermodynamic limit and by means of microscopic mod- 
els. 

1. Spatial self-organization into a Bragg-crystal 

A thermal cloud of cold atoms interacting with a sin- 
gle mode of a high-finesse Fabry-Perot cavity undergoes 
a phase transition upon tuning the power P of a far- 
detuned laser beam (wavelength A) which illuminates the 
atoms from a direction perpendicular to the cavity axis 
(Asboth et al, 2005; Domokos and Ritsch, 2002). Below 
a threshold power Per, the thermal fluctuations stabi- 
lize the homogeneous density distribution of the atomic 
cloud, and light which is scattered off the atoms into the 
cavity destructively interferes, rendering the mean cavity 
field amplitude to be zero. Above threshold, P > Pa, the 
atoms self-organize into a A-periodic crystalline checker- 
board order which is bound by the interference between 
the pump field and the macroscopic cavity field, resulting 
from Bragg scattering into the cavity mode. 

This self-organization effect can be described in terms 
of a semiclassical model similar to Eq. (13), generalized to 
many atoms. A set of variables Pj and Vj is introduced. 
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the index j = 1, . . . ,N labeling the atoms. For simplic- 
ity, the atomic motion is considered in two dimensions 
spanned by the cavity axis and the pump laser direction, 
with coordinates x and z, respectively. The equation of 
motion for the coherent cavity field amplitude a is given 
by 



Ac—Uo''^cos'^{kxj) a~ k+Fq ^ cos^(fca:j) 
j 3 

~ ^VcS cos(fca;j) cos{kzj) + ^q, , (23a) 



where the effective pumping strength of the cavity mode 
is denoted by r]es = 2^ , see Eq. (10). Due to the in- 
terference term J^j cos{kxj) cos(fczj) light scattering into 
the cavity vanishes for a homogeneous atomic density dis- 
tribution. It can be small even if all the atoms are max- 
imally coupled but the signs of the summands alternate. 
The light forces exerted on the individual atoms along 
the cavity and pump direction are given by 



-hUo\a\ 



._d_ 
dxi 



cos^ {kxj) 



d 



?i?7cff (a + OL*)— — cos(fcz,) cosikxj) + (,xi , (23b) 
OXi ■' 



= -nUQ{n/gf—cos^{kzj) 
ozj 

d 

- hrjcsia + ot*)jr- cos(/cXj) cos(fczj) + , (23c) 

These equations include Langevin noise terms ^q., S^xj^ 
and £^zj, defined by the non- vanishing second-order cor- 
relations, 

N 

iCU =K + J2^ocos^{kx,) , (24a) 
i=i 

(CnCa) = ihTodnSirj) cos{kxj) , (24b) 

dn£*{r,) d^£{r,) + a„£(r,) *(r,)] , (24c) 

with indices n,m = Xj,Zj. The noise terms associated 
with different atoms are not correlated. The complex 
dimensionless electric field is given by = cos(fca;)Q; -I- 
cos{kz)V,/g. 

In Fig. 10 a numerical simulation of the trajectories 
of 40 atoms during the first 50 fis of illumination are 
shown. The initial configuration is given by an ensemble 
of thermal atoms with random positions from a uniform, 
and velocities from a thermal distribution. The cavity 
mode initially is in the vacuum state (a = 0). With the 
right choice of parameters the emergence of a periodic 
pattern in the spatial density distribution of the atoms 
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FIG. 10 (Color online) Self-organization of laser-driven atoms 
in a cavity. Numerically simulated two-dimensional trajecto- 
ries during the first 50 /isec of transverse illumination. A 
checkerboard pattern of trapped atoms emerges, in which the 
occupied trapping positions are separated by even multiples 
of A/2 (= 1 edge). The grid lines indicate points of maximum 
coupling to the standing-wave cavity or pump field. There is 
a possible complementary configuration with atoms occupy- 
ing the other set of intersections. Parameters: 7 = 20///sec, 
{g, hi) = (2.5, 0.5)7, atomic detuning Aa = — 5OO7, cavity de- 
tuning Ac = —n + NUo, and the pumping strength f2 = 5O7. 
From Asboth et al. (2005). 



is observed, accompanied by the buildup of a coherent 
cavity field amplitude (Fig. 11). In the emerging config- 
uration, assuming a red-detuned pump laser, the trapped 
atoms are oscillating about intensity maxima of the in- 
terfering pump-cavity field. Along the cavity and pump 
direction these are separated by even multiples of the op- 
tical wavelength. Since only the black or white fields of 
the underlying checkerboard lattice pattern are occupied, 
constructive interference leads to efficient Bragg scatter- 
ing of pump photons into the cavity. 

As will be shown in more detail later, the process of 
self-organization relies on the right choice of the detun- 
ing 5c = Ac — NUo/2 between the pump laser and 
the dispersively shifted cavity resonance Ac. For the 
case Sc < 0, the potential term cos(A:a;j) cos(fc2;j) in 
Eqs. (23b, c) attracts atoms towards the "majority" sites 
and repels them from the "minority" sites, providing pos- 
itive feedback. Initiated by density fluctuations, one of 
the two possible Bragg lattices is then formed in a run- 
away process. For the case Sc > 0, the scattered cavity 
field creates potential maxima (minima) at the positions 
of the majority (minority) sites, counteracting the ampli- 
fication of density fluctuations and preventing a dynam- 
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FIG. 11 (Color online) Cavity cooling in the self-organized 
phase. The time evolution of the photon number in the cav- 
ity (left) and the phase space density of the atoms (right) on 
a long time scale, for A*' = 40 and A'^ — 160 atoms (differ- 
ent vertical scalings are used for illustrating the superradi- 
ance |ap oc A''^). In the right panel, the green curve fluc- 
tuating around a constant value corresponds to a uniformly 
distributed A'^ = 40 atoms driven below the self-organization 
threshold. The parameters are the same as in Fig. 10. From 
Asboth et al. (2005). 



ical instability. Furthermore, in this regime the delayed 
cavity response causes cavity heating of the atomic mo- 
tion which obscures an equilibrium situation in the lack 
of other dissipative process. 

For Ac — NUq < 0, the initial fast buildup of a coher- 
ent cavity field continues over a longer timescale. The 
kinetic energy of the oscillating and the untrapped atoms 
dissipates owing to the cavity cooling mechanism, which 
leads to an increase of the number of trapped atoms and 
a stronger localization in the potential wells. This further 
improves coherent scattering into the cavity, as indicated 
by the slow increase in the cavity field intensity shown 
in Fig. 11. Comparing the time evolution of the intra- 
cavity photon number for the self-organization process 
of 40 and 160 atoms (rescaled in Fig. 11 by a factor of 
16) demonstrates the superradiance effect, i.e. the field 
intensity scales cooperatively as the square of the par- 
ticle number. The right panel of the figure shows that 
the cooling rate, described by the decrease of the phase 
space density of the atoms, is also similar for A^ = 40 and 
N = 160 and that the self-organization leads to smaller 



phase space densities than the homogeneous distribution 
below threshold. 

Self-organization of laser-cooled atoms has been ob- 
served in experiments of the MIT group with N « 10^ 
Cs atoms prepared at a temperature of 6 /^K in a nearly 
confocal Fabry-Perot cavity (Black et al., 2003). Above 
a threshold intensity of the transverse pump beam, col- 
lective emission of light into the cavity was observed at 
a rate which exceeds the free-space single-atom Rayleigh 
scattering rate by a factor of up to 10^. This experiment 
clearly demonstrated the process of spontaneous symme- 
try breaking by measuring 7r-jumps in the phase of the 
emitted cavity field relative to the transverse pump field, 
corresponding to self-organization into the black or white 
lattice sites of a checkerboard pattern (Fig. 12). Retar- 
dation between the cavity field and the atomic motion re- 
sulted in a collective friction force on the center-of-mass 
degree of freedom. A deceleration of up to 1000 m/s^ 
has been achieved with atom-cavity detunings as large 
as Aa = -2tt X 1.58 GHz. 

For finite atom number N and finite measurement 
time, an interesting hysteresis effect accompanies self- 
organization, as shown in Fig. 13. The thermodynamic 
limit iV ^ cx) is approached by simulations of Eq. (23) 
with the atomic density N/V cx Ng^ and the cavity loss 
rate kept constant. The percentage of defect atoms af- 
ter 4 ms of simulation time as a function of the pump- 
ing laser strength clearly shows the transition. However, 
the transition point is dependent on N and on whether 
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FIG. 12 Observation of spontaneous symmetry breaking in 
the self-organization phase transition. Simultaneous time 
traces of the (a) intracavity intensity (in arbitrary units) 
and (b) the relative pump-cavity phase. Drops in the in- 
tensity correspond to time intervals during which the beams 
of the magneto-optical trap (MOT) are switched on forcing 
the atomic density distribution to randomize. After switch- 
ing off the MOT beams, the atoms self-organize again into 
one of the two possible checkerboard patterns, as indicated 
by the relative phase signal. Experimental parameters were 
N = 8.2x 10^ Aa = -27r x 1.59 GHz, Ac = -2n x 20 MHz, 
and ///sat = 440. From Black et al. (2003). 
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FIG. 13 (Color online) Hysteresis effect for finite measure- 
ment time. Ratio of atoms in the "defect" positions against 
pumping strength rj, 4ms after the loading of the trap with 
a uniform ("up") or organized ("down") gas of atoms. The 
different curves show the approach towards the thermody- 
namic limit. The parameters are k = 7/2, Aa ~ — 5OO7, 
Ng'^ = 2007^ Ac = ~K- Nq'^/IAa], ksT = hn. From 
Asboth et al. (2005). 

the initial positions were uniformly distributed ("up") 
or at "odd" points of maximal coupling ("down"). The 
breadth of the hysteresis increases with the atom num- 
ber, but decreases with the measurement time. This 
behavior can be explained by taking into account sta- 
tistical fluctuations arising from the finite atom number 
N. Assuming that the self-organization from a uniform 
distribution ("up" curves) is triggered when the fluctu- 
ating energy difference between the even and odd sites 
momentarily exceeds the mean kinetic energy, a scaling 
Ng'^ of the threshold with the pump intensity was found 
in accordance with the numerical results of the plot. The 
disappearance of the lattice pattern for decreasing pump 
power ("down" curves) when the system is started from 
the ordered phase occurs at the half of the mean-field 
threshold (see Eq. (33) in the next subsection), indepen- 
dently of the atom number N. 

Self-organization of laser-driven atoms occurs also in 
a transversally driven ring-resonator geometry support- 
ing two running- wave modes (Nagy et ai, 2006b). In 
contrast to the linear single-mode cavity case, here the 
transition from the homogeneous to the organized density 
distribution involves spontaneous breaking of a continu- 
ous (rather than a discrete) translational symmetry (see 
also IV.D.3). 



2. Collective atomic recoil laser 

Collective atomic recoil lasing (CARL) is the promi- 
nent many-body instability effect in a ring cavity, origi- 
nally predicted by Bonifacio et al. (1994). An ensemble of 
cold atoms couples to two counter-propagating modes of 
a unidirectionally pumped high-flnesse ring cavity. Light 
scattering off the atomic ensemble between these cavity 



modes leads to a collective instability corresponding to 
an exponential gain for the back-propagating field mode 
amplitude in conjunction with an atomic bunching at 
the antinodes of a self-organized optical lattice. In the 
presence of dissipation of the atomic kinetic energy, a 
steady-state operation of CARL can be achieved with a 
self-determined atomic drift velocity and back-reflected 
light frequency. 

The CARL scheme involves an interplay between the 
influence of the atomic motion on the Rayleigh scatter- 
ing of light and, reversely, the mechanical effect of light 
upon the atomic motion. The former effect has been 
previously seen, for example, as the so-called recoil in- 
duced resonance (RIR) in the transmission spectrum of a 
probe beam making a small angle with a one-dimensional 
lin_Llin optical molasses (Courtois et ai, 1994), and also 
in a optical dipole trap formed by counter-propagating 
modes of a ring cavity (Krusc et ai, 2003b). This nar- 
row, dispersion-like resonance around the pump field fre- 
quency originates from a two-photon Raman transition 
between different momentum states of the atoms. Differ- 
ent populations of the corresponding momentum states 
lead to gain or attenuation of a probe beam. For a 
thermal velocity distribution, a probe frequency tuned 
slightly below the pump frequency gives rise to gain, 
whereas for negative detuning there is an attenuation of 
the probe. The probe transmission spectrum measure- 
ment provides information about the temperature, and 
even more about the velocity distribution (Brzozowska 
et al, 2006). 

It has been predicted that, based on the RIR gain ef- 
fect, a weak probe field injected in the direction oppo- 
site to the strong pump field would exponentially am- 
plify due to the self-bunching of a part of the atoms 
into a lattice which Bragg reflects the strong pump beam 
more efficiently than single atoms (Bonifacio et at, 1994). 
However, one needs long enough interaction time so that 
back-action of the light scattering on the velocity dis- 
tribution can have a significant effect (Bcrman, 1999). 
This can be accomplished, for example, by confining 
the light modes into a cavity. Lasing mediated by the 
collective atomic recoil between the counter-propagating 
modes of a unidirectionally ring cavity has been observed 
(Kruse et ai, 2003a). Back-action on the atomic mo- 
tion has been demonstrated by detecting the displace- 
ment of the atoms accelerated by the momentum trans- 
fer process. The self-consistent solution is an accelerat- 
ing Bragg-lattice of atoms co-moving with the standing 
wave formed by the pump and the back-reflected com- 
ponent having a Doppler-shifted frequency Aw = 2kv 
with respect to the pump. The phase dynamics of the 
counter-propagating modes can be monitored as a beat 
signal between the outcoupled beams, which reveals the 
acceleration by an increasingly red-detuned probe as a 
function of time (Fig. 14). 

The runaway process can be counteracted by introduc- 
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FIG. 14 Observation of lasing mediated by collective atomic 
recoil, (a) Recorded time evolution of the observed beat sig- 
nal between the a+ and a_ cavity modes of a ring cavity. 
Initially both modes are pumped to form an optical lattice of 
— 10^ atoms. The initial drop is due to the decay of the 
unpumped mode after switching off the a_ pump a.t t — 0. 
Well beyond the ring-down time of about 10 /xs, the persist- 
ing oscillations demonstrate the coherent backscattering of 
the pumped mode, (b) Numerical simulation with the tem- 
perature adjusted to 200 /iK. (c) The symbols (x) trace the 
evolution of the beat frequency after switch-off (dotted line 
is from numerical simulation). The increase of the beat fre- 
quency corresponds to the acceleration of the Bragg lattice of 
back reflecting atoms, which lasts until the Doppler-shifted 
frequency drops out of cavity resonance. Absorption images 
of a cloud of 6 X lO'^ atoms recorded at (d) ms and (e) 6 ms 
after switching off the probe beam pumping. All images are 
taken after a 1 ms free expansion time, (f ) This image is ob- 
tained by subtracting from image (e) an absorption image 
taken 6 ms after switch-off with low cavity finesse for which 
no collective recoil is expected. The intracavity power has 
been adjusted to the same value as in the high finesse case. 
From Kruse et al. (2003a). 



ing some external friction force on the atom motion. The 
dissipation gives rise to a steady-state solution which in- 
volves a constant drift of the entire atomic cloud at a 
speed V. Accordingly, we transform the atomic position 
variables such + vt, and the coherent field am- 

plitude of the running-wave field mode propagating op- 



posite to the pumped field mode as q:_ — a^e 



2ikvt 



The 



drift velocity v is to be determined in a self-consistent 
manner. The semiclassical equations describing atomic 
motion are given by 



Pj 



Pj + KU{) 2ik [a'^c 



2ikz 



(25a) 

where /3 represents the linear friction arising, e.g., from 
collisions with a buffer gas (Bonifacio, 1996; Pcrrin et ai, 
2001) or laser-cooling in an optical molasses (Kruse et ai, 



2003a). The cavity field amplitudes evolve as 

= {iSc - k) a+ - iNUoBa^ + rj , (25b) 
(5_ = [i {Sc - 2kv) - k] 5_ - iNUoB*a+ , (25c) 

where Sc — Ac — NUq is the effective detuning of the 
pump frequency from the atoms-shifted cavity resonance. 
The atomic positions enter through the bunching param- 
eter 



1 ^ 
N ^ 



(25d) 



A closed set of equations can be formed in which the 
atomic cloud is characterized by the three real param- 
eters w, 6, and 0. A trivial solution of these equations 
corresponds to the case where the atoms are uniformly 
distributed in space (6 = 0, w = 0) and the counter- 
propagating field mode amplitude vanishes (a_ = 0). A 
non-trivial steady-state solution can be obtained numer- 
ically from the coupled algebraic equations (time deriva- 
tives set to zero). This solution can be approximated 
analytically by assuming perfect bunching, h — \ and 
ip — The only remaining free parameter is then given 
by the steady-state drift velocity v which obeys the alge- 
braic equation 



Ikv = 8 



D 



(26) 



where D = {i5c - n)[i{5c - 2fcw) - k] + N^U^b^. In 
order to gain insight into the solution, the following sim- 
plifications can be made: (i) neglect in D the last Uq 
term originating from (second-order) scattering of the 
Q!_ mode back into the pumped a+ one, and (ii) con- 
sider resonance 6c = 0. Then one can identify the typ- 
ical solutions of two different regimes. In the limit of 
a large Doppler shift, kv ^ k, the drift velocity and 
the back reflected power scale with the atom number as 
V oc N^^^ and oc A^**/^, respectively. This is re- 

ferred to as the CARL limit in the literature. In the op- 
posite limit of a small Doppler shift, kv <^ k, the velocity 
obeys v (x N and the intensity exhibits superradiant be- 
havior, |a_P oc N'^, corresponding in this geometry to 
Bragg retro-reflection. The relation between these two 
kinds of superradiant instabilities in the collective inter- 
action of light with an atomic gas has been established 
experimentally by Slama et al. (2007a,b). While superra- 
diant Rayleigh scattering from atomic clouds is normally 
observed only at very low temperatures, i.e., well below 
1 /iK, (Inouye et al., 1999), the presence of the ring cav- 
ity enhances cooperativity and allows for superradiance 
with thermal clouds as hot as several 10 /iK. 

In the experiments of the Tiibingen group (Cube et al., 
2004; Kruse et al., 2003a), an optical molasses has been 
used to impose a motional damping force on the atoms. 
In such a viscous CARL system, there is a steady-state 
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FIG. 15 Pump power threshold for the collective atomic recoil 
lasing. (a) Sectionwise Fourier-transform of the interference 
signal Pbcat with the pump laser power being ramped down 
and up (the dotted line is proportional to the pump laser 
power). At t = the system is in the ordered CARL phase, 
by gradually decreasing the power the drift velocity decreases 
(the peaks in the Fourier spectra shift downward) until the 
back reflection ceases at a threshold pump power. Ramping 
up the pump power from about i = 5 ms, the appearance of a 
peak is delayed and occurs at about the same threshold value 
of the pump power, (b) Dependence of the CARL frequency 
on the intracavity pump power, (c) Dependence of the probe 
field intensity on the pump power. The CARL laser threshold 
is around P+ = 4 W intracavity power. The fitted curves are 
based on a Fokker-Planck theory outlined in Sec. III.C.I. (d) 
Calculated bunching parameter. The parameters are /3 = 4k, 
= lO'' atoms, Aa = -27r x 1.7 THz and T = 200^tK. From 
Cube et al. (2004). 



operation with a self-consistent drift velocity, according 
to Eq. (26), at which the friction compensates the accel- 
eration due to back scattering of photons. The Fourier 
spectrum of the beat signal between the pumped mode 
and the counter-propagating one in steady-state is shown 
in Fig. 15a for various pump strengths. The drift veloc- 
ity can be deduced from the Doppler-shift of the back- 
reflected field as shown in Fig. 15b, the corresponding 
range of velocities is from 7 to 13cm/s. The beat fre- 
quency as well as the drift frequency vary as a function 
of the pump power. However, the dramatic feature of 
Fig. 15 is the clear appearance of a threshold: the self- 
bunching and backscattering starts only above a well- 
defined threshold pump power. This measurement pro- 
vides thus an experimental evidence of a phase transition 
to a state of synchronized atomic motion (Cube et at, 
2004). Underlying the critical behavior is the diffusion 
accompanying friction: the spontaneous photon scatter- 
ing of the molasses laser beams leads to a random heating 
force which stabilizes the homogeneous distribution and 
thereby can prevent the formation of a Bragg-lattice for 
weak pump power (Robb et ai, 2004). Above thresh- 



old, the dissipation and the fiuctuations together lead to 
a position distribution which exhibits a finite bunching 
parameter. These effects are discussed in the framework 
of various mean-field theories for the atomic position dis- 
tribution (see Sect. IILC). 

We note that the CARL system and its phase transi- 
tion has been studied also in the case where the atomic 
transition frequency is close-to-resonance with the pump 
laser field (Perrin et ai, 2001). Then the atomic polar- 
ization plays a dynamical role and the transition does not 
require spatial bunching but the emergence of coherent 
polarization grating (Perrin et ai, 2002). 

The viscous CARL transition exhibits analogy with 
that of the generic Kuramoto model which describes the 
self-synchronization of coupled oscillators with different 
frequencies (Kuramoto, 1975; Strogatz, 2000). To re- 
veal the analogy, one can transform the CARL equa- 
tions using the following assumptions: (i) the motion is 
overdamped {pj = 0); (ii) the pumped field amplitude 
a+ is a constant of time (a+ « ??/«:); (iii) the counter- 
propagating mode amplitude is stationary, oscillating at 
frequency loq, i.e. , effectively Q!_ = —iuQci- with the 
frequency wq being determined by the constant drift ve- 
locity; and (iv) k luq. With these assumptions and by 
using the notation 6j = 2kzj, Eqs. (25) simplify to 



2k 

m 13 



^ + Kb sin((p - Oj) , 



(27) 



which is formally equivalent with the Kuramoto model. 
The Langevin-type random noise ^, associated with the 
friction term ~f3p in Eq. (25a), introduces the random 
frequencies present in the Kuramoto-type systems. The 
coupling strength is i^T = 2lujiNUq \a+\'^/{ujQ (3). The 
mean-field character is obvious: each oscillator couples 
only to the mean-field quantities b and ip. The phase 9j 
is pulled toward the mean-field phase ip, which increases 
the order parameter b. The coupling is proportional to b, 
which sets a positive feedback loop. With the increasing 
coherence b, even more oscillators can be recruited to 
the synchronized pack (those being within the bandwidth 
Kb) , further increasing b. Such a runaway process starts 
only above a critical coupling K. 



C. Phase-space and mean-field descriptions for large 
particle numbers 

Due to the nonlinearity of the coupled atom-field dy- 
namics, exact analytic results for atomic ensembles cou- 
pled to optical cavities are rather sparse and the com- 
putational demand often hinders simulations at realistic 
particle numbers (Deachapunya et aZ., 2008; Salzburger 
and Ritsch, 2009). In the following, we will present mean- 
field methods which allow to predict instability thresh- 
olds and to model the effective dynamics in the thermo- 
dynamic limit. Starting with the assumption of a quasi- 
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thermal distribution, the critical behavior can be revealed 
and the threshold for criticality can be approximated. We 
will continue by adopting phase-space methods that can 
account for general position and velocity distributions. A 
Vlasov-type equation allows for a more accurate estima- 
tion of the threshold, and also for performing stability 
analysis and setting up a phase diagram. This analy- 
sis leads to such intriguing predictions that the cavity- 
mediated interaction combined with the cavity cooling 
effect can be the basis of a generally applicable sym- 
pathetic cooling scheme. Finally, a Fokker-Planck-type 
equation can be constructed in order to determine the 
steady-state of the non-equilibrium systems in the ther- 
modynamic limit. The methods show similarities with 
plasma physics where equally complex, coupled dynam- 
ics of particles and fields occurs (Montgomery, 1971). 



1. Critical point 

The simplest mean-field model is based on the assump- 
tion that the atomic motion is overdamped and the distri- 
bution function of the atomic positions p{x^ t) is a ther- 
mal distribution. The critical point of the CARL in- 
stability (Sec. III.B.2) and that of the self-organization 
(Sec. III.B.l) can be calculated with this approach. 



a. Dynamical equations The motional damping is char- 
acterized by a linear friction coefficient ji (half of the 
kinetic energy damping rate is /3/m with m being the 
atomic mass) and a temperature T. Then the mean atom 
density distribution obeys the Smoluchowski-equation 



dp{.x, t) 
dt 



ld_ 



F{x)p{x,t)~kBT 



dp{x,t) 
dx 



(28) 



In the ring cavity geometry of CARL (Kruse et ai, 
2003a), for example, the force F{x) is given by the last 
term on the right-hand side of Eq. (25a). It contains the 
field mode amplitudes which couple back to the atomic 
density distribution via the bunching parameter B, 



with B= f dxp{x,t)e^''''' . 
Jo 



(29a) 
(29b) 



For simplicity, the center-of-mass velocity and the de- 
tuning Sc were set to zero. Linear perturbation calculus 
leads to the instability threshold of the homogeneous so- 
lution (Robb et ai, 2004). Since the spatial coupling 
functions are sinusoidal, only a few Fourier components 
are involved in the initial dynamics. In particular, in 
order to determine the instability of the homogeneous 
distribution, only the single mode function e"'*^^ of the 
counter-propagating cavity mode needs to be taken into 
account. In the so-called CARL limit (see Sec. III.B.2), 



one obtains the following threshold condition for the cav- 
ity pump amplitude: 

3/2 1 ^ ^5/2 



V 



> 



(29c) 



b. Canonical distribution A further simplification can be 
made in the mean-field approach if k is the largest rate in 
the dynamics. Adiabatic elimination of the cavity field 
dynamics then results in a self-consistent optical poten- 
tial V{x) in which the spatial density p{x) of the atoms 
is determined by a canonical distribution, 

p{x) = ^exp(~V{x)/{kBT)), (30) 

with the partition function Z = J c:iq){—V{x)/{kBT))dx 
ensuring normalization of p{x) to unity. The tempera- 
ture could be identified with the one which is achieved in 
cavity cooling, kgT w hn, but in general it can be set by 
other means, e.g., by laser-cooling in an external optical 
molasses. The nonlinearity enters the equations through 
the dependence of the potential V{x) on the atomic den- 
sity p{x) itself, i.e. , V = V{x,p{x)). 

In principle, the optical force acting on the atoms in 
the cavity does not derive from a potential when the 
back-action of the atomic motion on the radiation field 
amplitude is significant. As was noted by Asboth and 
Domokos (2007), dynamical equations based on forces, 
such as Eq. (28) and Eq. (29a) have to be used. How- 
ever, in the spirit of the mean-field approach, the effect 
of an individual atom on the field amplitude is negligi- 
ble with respect to the summed effect of all the others. 
In the limit of many atoms with small single-atom cou- 
pling, the motion of a single atom is very well described 
by an effective potential determined by the many-body 
ensemble. 

For the example of self-organization in a standing- wave 
cavity, see Sec. III.B.l, the light potential along the cav- 
ity axis is given by 



V{x) = U2 cos'^{kx) + Ui cos{kx) , 



(31) 



which is composed of the sum of a A/2-periodic poten- 
tial stemming from the cavity field and a A-periodic one 
arising from the interference between cavity and pump 
fields. The depths of these potentials are 



U2 = N^{cos{kx)yhIoUo 
Ui = 2N{cosikx)) hlo {Ac 



(32a) 

NUo{cos^{kx))) , (32b) 



where Iq is a dimensionless single-atom scattering pa- 
rameter, Iq (x rf' . Equation (30) has to be solved in 
a self-consistent manner by iteration. In Asboth et al. 
(2005), the threshold has been analytically determined 
to be 



vis, 



kBT{K^ + 6l) 



N\6, 



c\ 



(33) 
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2. Stability analysis and phase diagram 

In the following we will account explicitly for the ef- 
fect of the velocity distribution on the dynamics and 
on the instability threshold. A mean-field model based 
on the Vlasov-equation for the phase space distribution 
f{x,v,t) has been derived from a microscopic theory for 
the infinite system size by Griefier et al. (2010). For one- 
dimensional motion along the cavity axis, the dynamical 
equation reads 



dl 
dt 



dx 



(34) 



where (j){x,a) is the potential corresponding to a mo- 
mentary field amplitude a. For the generic example of a 
laser-driven cold atomic cloud in a single-mode standing- 
wave resonator with mode function cos(fca;), the potential 
is 



{x, a) — 



2h /[/, 



|ap cos(2fca:) -t- ?7cffRe (a) cos{kx) 



(35) 

Similarly to the adiabatic potential, Eq. (31), there is 
a A/2 and a A-periodic term, this latter originates from 
interference between the transverse pump laser and the 
intracavity field. However, in this more general approach 
the cavity field amplitude is kept dynamical obeying the 
self-consistent equation 

a = (— K + iSc)a + rj 

— ia — - — / dv cos{2kx)f{x,v,t)dx 

2 J-oo Jo 

/oo /"A 
dv / cos{kx)f{x,v,t)dx , (36) 
-oo Jo 

where Sc — Ac — NUq/2. This approach based on the 
Vlasov-equation is well suited to study the mean-field 
dynamics at short times in order to test the stability of 
stationary states. Over longer time scales, statistical fiuc- 
tuations have to be taken into account in the framework 
of a Fokker-Planck equation for the velocity distribution, 
which is presented in Sec. III.C.3. 



a. Nonlinear response of a cold atomic cloud in a driven 
Fabry-Perot cavity For cavity pumping only {rjeS = 
0,77 ^ 0), the self-consistent steady-state solution ex- 
hibits a strong nonlinear optical response. Underlying 
the nonlinearity, the particle distribution and thus the 
effective refractive index of the cloud depends on the cav- 
ity pump intensity. Above a sufficient pump strength, 
multiple stationary solutions appear, reminiscent of op- 
tical bistability.One can perform a systematic stability 
analysis of these solutions by studying the dynamics of 
small fluctuations of the field and the particle distribu- 
tion (Griefier and Ritsch, 2011). As shown in Fig. 16, the 
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FIG. 16 (Color online) Normalized solutions for the steady- 
state photon number (red) 7o = |ap versus effective cavity 
detuning 5 = Sc^ NUq for a thermal gas in a driven standing- 
wave cavity. The driving strength is i] — 13k (a) and rj = 18k 
(b). Those parts of the response curve that lie inside the 
instability region (shaded area) are depicted in dashed and 
correspond to linearly unstable steady states. The intervals 
designated A correspond to bistability, the interval designated 
B supports no stable steady state at all. The parameters are 
N = 10^ Uo = 0.04:K, rj = 18k, k = 2000aJii and fcsT = fin. 
From Griefier and Ritsch (2011). 



stability analysis reveals regions of bistability as well as 
parameter ranges where no stable solutions exist. 

In the parameter regions of instability, the numerical 
solution of the dynamical Vlasov equation reveals a limit 
cycle behavior with subsequent appearance of higher fre- 
quencies than the fundamental cycle (Griefier and Ritsch, 
2011). In the quantum regime essentially the same be- 
havior is retrieved with the recoil frequency determining 
the oscillation frequencies v « 4a; (Ritter et al, 2009) 
(see Sec. IV.C). 



b. Self-organization of a laser-driven cloud of atoms For 
a purely transverse pump geometry (r] — 0,7]^^ 7^ 0), 
one can systematically recalculate the critical pump am- 
plitude ?7off,c that marks the transition from the stable 
regime to the unstable one, where small fluctuations are 
amplified and grow exponentially. 

The Vlasov equation Eq. (34) together with the equa- 
tion for the coherent cavity field amplitude a possesses 
an infinite number of stationary solutions with a spa- 
tially homogeneous density distribution and zero cavity 
field but different velocity distribution, which, however, 
are not necessarily stable against fiuctuations. Indeed, 
any symmetric velocity distribution g{v/vT) — Lv'Yf{v) 
for < is stable only if 



vp 



-2^ 



h\Sc 



(37) 



where vp denotes the Cauchy principal value. Here we 
have defined the thermal velocity — 2k-QT/m; L de- 
notes the cavity length. For a Gaussian distribution the 
integral evaluates to one, and the condition is equivalent 
to Eq. (33). 
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FIG. 17 (Color online) Laser illuminated cold gas in a ring- 
cavity: time evolution of the intra-cavity field intensity (left), 
and the instantaneous spatial n{x) (upper right) and veloc- 
ity F{v) (lower right) distributions of the particles along the 
cavity axis at times t — and = 28/k for 6c = — k (upper 
row) and 5c = n (bottom row). For 5c < 0, after a transient 
exponential growth, the field intensity saturates accompanied 
by the trapping of atoms in the self-organized pattern. By 
contrast, for 5c > the transient exponential growth is fol- 
lowed by oscillations. Parameters are A'^ = 10'', Uo = ~k/N, 
77cff = 0.05k, kvT = K, and vr = vt/5. From Griefier et al. 
(2010). 
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FIG. 18 (Color online) Simultaneous self-organization of two 
species. The system is started from a perturbed uniform state 
above the instability threshold, Eq. (39), in such a way that 
species one (two) itself would be pumped six times above (far 
below) the critical point. Figure (a) shows the position dis- 
tributions in the final state, (c) and (d) the momentum dis- 
tributions initially (dashed lines) and after self-organization 
(solid line). Solid lines depict the results of stochastic tra- 
jectory simulations for ensembles of particles as in Eq.(23b), 
while open circles show the predictions of the corresponding 
Vlasov model. Figure (b) shows the time evolution of the two 
order parameters 6i{2) approaching theoretical steady-state 
values. Parameters are A^i = 10*, N2 = 500, m2 = lOmi, 
ksTi = 10**^K, kBT2 = 2.5 x lO^hn, 7?i = 2.4k, rj2 = 27.4k 
and iUR = 10"^K. From Griefier et al. (2011). 



Figure 17 shows the results of a numerical simulation 
of Eqs. (34), (35) and (36), initialized with a perturbed 
gaussian distribution 

f{x,v,0) = , ^ e-'''/^^(l-ecos(fca;)) , (38) 

with e <C 1, for the case of a transversely pumped ring 
cavity, where the light can be scattered into a superposi- 
tion of two resonant cavity modes. Apart from possess- 
ing continuous translational symmetry, the dynamics is 
qualitatively very similar to the single mode case (Griefier 
et al., 2010) shown in Fig. 11 of the previous section. One 
clearly recognizes the striking difference in the dynamical 
behavior for positive and negative values of 6c- While we 
have an instability in both cases, self-organization is only 
found for Sc < 0- 

c. Sympathetic self-ordering and cooling Self-organization 
and collective coherent light scattering into a high-finesse 
cavity in principle allows for trapping and cooling of any 
kind of polarizable particles. In practice, however, the 
required phase-space densities and laser intensities to ini- 
tiate the ordering process are hard to achieve for atomic 



species, molecules or nanoparticles which cannot be ef- 
ficiently optically precooled (Deachapunya et al., 2008; 
Lev et at, 2008). As an alternative approach, the self- 
organization threshold can be achieved by inserting dif- 
ferent species simultaneously into an optical resonator. 

The Vlasov-type model can be generalized to a dilute 
gas of various kinds of Ng polarizable point particles of 
mass nis illuminated by a single transverse standing- wave 
laser field. For the threshold value, a condition analogous 
to Eq. (37) can be obtained (Griefier et al, 2011). The 
homogeneous distribution is unstable if and only if 

where kBTs — msv1/2. Note that the right-hand side 
of Eq. (39) depends only on cavity parameters, and all 
terms in the sum on the left-hand side are positive and 
proportional to the pump intensity. This guarantees that 
inserting any additional species into the cavity always in- 
creases the total light scattering rate and thus lowers the 
minimum power needed to start the self-organization pro- 
cess, regardless of temperature and polarizability or den- 
sity of the additional particles. Moreover, the different 
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FIG. 19 (Color online) Sympathetic cavity cooling. Time 
evolution of the kinetic temperatures of a heavy and a light 
species. The blue dashed line represents the heavy particle 
alone and the blue solid line the enhanced cooling in the 
presence of a lighter species below self-organization thresh- 
old. Parameters are m2 = 200mi, iVi = 200, N2 = 200, 
y/Nlrji = HAujR, y/N2rj2 = ISAljr, k = 200wi{ and dc ~ — k. 
From Griefier et al. (2011). 



species can be located at different regions within the cav- 
ity. Assisted self-organization of a species which, alone, 
would be pumped below threshold is shown in Fig. 18. 

Below the self-organization threshold cooling occurs 
thereby equalizing the stationary momentum distribu- 
tions for all species. Fig. 19 exhibits the enhanced decay 
of the kinetic energy of the heavy particles in presence 
of cavity field and a cold species. Although the distri- 
butions get independent in stationary equilibrium, the 
cooling process itself involves energy exchange between 
different species. Thus if any of the species is cold or can 
be cooled by different means, the other components are 
sympathetically cooled in parallel. 



3. Non-equilibrium steady-state distributions 

Over longer time scales, diffusion has to be accounted 
for in terms of a nonlinear Fokker-Planck equation for the 
statistically averaged velocity distribution. This allows 
for the calculations of cooling time scales and the unique 
steady-state distribution. 



a. Transverse pump configuration below threshold Below 
the instability threshold, Eq. (37), the mean spatial dis- 
tribution / is homogeneous, i.e., independent of x. The 
statistical fluctuations of the potential and the actual 
atomic distribution gets important. For the spatially 
averaged distribution, a lengthy calculation (Nicdenzu 
et al., 2011) leads to a nonlinear Fokker-Planck equation 
for the velocity distribution F{v,t) = f{x,v,t), 



d 



d_ 
dv 



[a[F]i 



F = 



d_ 
dv 



(40) 



with coefficients 

B\F] = 



|D(ifcu)|2 



2m2 



\D{ikv)\^' 



(41a) 
(41b) 



These functionals depend on (F) via the dispersion rela- 
tion 



D{s) - (5 
— ihkSc 



n? 



2m 
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dv 



F'{v) 
s + ikv 



F'{v) 
s — ikv 



(42) 



which encodes all cavity-mediated long-range particle in- 
teractions. Far below threshold the dispersion relation 
reduces to D{ikv) ~ {ikv + k)^ + 6^, which corresponds 
to the case of independent particles. 

Steady-state solutions of Eq. (40) exist only for neg- 
ative detuning Sc < 0, where light scattering is accom- 
panied by kinetic energy extraction from atomic motion. 
Below threshold one obtains non-thermal q-Gaussian ve- 
locity distribution functions (de Souza and Tsallis, 1997): 



F(t.)(X l-(l-g); 



2fcBT 

with q = 1 + iOf(/\Sc\ and the effective temperature 



(43) 



kvtT — h- 



4|<5c 



(44) 



The minimum temperature is reached for dc = —k. The 
magnitude of the detuning \Sc\/i^r determines the shape 
of the distribution. For \Sc\ = lor, it is a Lorentzian 
distribution, whereas for \5c\/'-0r — > 00, i.e. g — )■ 1, it 
converges to a Gaussian distribution with kinetic tem- 
perature /ceTkin = m (w^) . 

Inserting the steady-state g-Gaussian (43) distribution 
into the threshold condition (37) gives a self-consistent 
stability criterion. As a result, the homogeneous distri- 
bution is stable only if 



(45) 



q 



where the equality is reached for Sc = —k,, i.e., for op- 
timum cavity cooling given by Eq. (14). The stability 
criterion can be rewritten in the intuitive form 



N\Uo\Vp < 



(46) 



where Vp = ft^/AA is the optical potential depth created 
by the pump laser, and NUq is the total dispersive shift 
of the cavity resonance. Note that even if the initial 
temperature is too high for the homogeneous distribution 
to be unstable, cavity cooling induced self-organization 
is possible. 
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FIG. 20 (Color online) Schematic view of the phase di- 
agram in the weak-coupling limit (Af|C/o| «^ «) for k = 
lOOojij. Equilibrium solutions exist only for 5c < —uir/2, 
a Lorentzian steady-state velocity distribution is realized for 
the case \5c\ = ijJR- For large negative values of the detun- 
ing 5c, strongly organized equilibrium solutions exist already 
for pump strengths slightly above the critical value. Adapted 
from Niedenzu et al. (2011) 



b. Transverse pump conGguration above threshold Above 
the self-organization threshold, the inhomogeneous spa- 
tial distribution can still be derived from a Fokker-Planck 
equation similar to Eq. (40) by using action-angle vari- 
ables (Chavanis, 2007; Luciani and Pellat, 1987). In the 
limit of deep trapping and a harmonic approximation for 
the potential, the steady state is a thermal distribution 
with a temperature depending both on the effective trap 
frequency ujq and on the cavity linewidth k, 



5^ + 4uj^ Sc^'^ 



4|(5c 



hujt 



'0 ■ 



(47) 



IV. QUANTUM GASES IN OPTICAL CAVITIES 

Quantum gases are considered as ideal model sys- 
tems to study quantum many-body phenomena under 
well-controlled experimental conditions. The possibilities 
which arise from loading ultracold atomic ensembles of 
different particle statistics into various optical potential 
landscapes and to tune the strength of the contact atom- 
atom interaction, make these system well suited for quan- 
tum information and simulation research (Bloch et al., 
2008) . The merger of the field of ultracold gases with that 
of cavity QED provides a set of additional possibilities. 
Cavity-mediated atom-atom interactions can be tailored 
by choosing different resonator and pump geometries and 
give rise to novel quantum phases. Closely related, the 
atomic back-action upon the cavity-generated lattice po- 
tentials can be significant, which paves the way to study 
phonon or soft-condensed matter physics with ultracold 
gases (Lewenstein et al., 2007). Further, coherent scat- 
tering into the cavity field can be used for non-destructive 
and real-time probing of different many-body phases. 

The coupling between a Bose-Einstein condensate and 
an optical cavity is conceptually fundamental since a sin- 
gle mode of a matter wave field interacts with a single 
mode of the light field: as all atoms occupy the same 
motional quantum state they couple identically to the 
optical cavity field. This situation can substantially re- 
duce the number of degrees of freedom necessary to de- 
scribe the system. Therefore the experimental situation 
can often be almost perfectly described by fundamen- 
tal Hamiltonians of matter-light interaction. These in- 
clude the Tavis-Cummings or Dicke model, as well as the 
generic model for cavity optomechanics. 



The effective trap frequency loq can be approximated by A. Experimental realizations 




(48) 



which is valid in the regime |i5c| ^ ojr with r]eg,c given 
by the equality of Eq. (45). As the temperature depends 
explicitly on the pump strength, a stronger pump laser 
beam results in more confined particles with increased 
kinetic energy. The system has the interesting property 
that the more particles we add, the deeper the optical 
potential gets, which shows analogy to self-gravitating 
systems (Posch and Thirring, 2005). 

Finally, putting all this together one obtains the self- 
consistent phase diagram for self-organization which ac- 
counts for cavity cooling, see Fig. 20. 



Experimentally, there have been different approaches 
to realize and study Bose-Einstein condensates or bosonic 
atomic ensembles close to quantum degeneracy in opti- 
cal high- finesse cavities (Brennccke et al. , 2007; Colombe 
et al, 2007; Gupta et al, 2007; Purdy et al, 2010; Slama 
et al., 2007a). So far, all groups used *^Rb atoms. In the 
Tiibingen group (Slama et al., 2007a) a Bose-Einstein 
condensate was loaded for the first time into a ring cav- 
ity with large mode volume using magnetic trapping and 
transport. This experiment extended prior work on the 
collective atomic recoil laser (see Sect. III.B.2) with laser- 
cooled atoms (Cube et al., 2004; Kruse et al., 2003a) into 
the ultracold regime. 

Loading ultracold quantum gases or Bose-Einstein con- 
densates into ultra-high finesse optical cavities of small 
mode volume, which operate in the single-atom strong 
coupling regime of cavity QED, has been achieved by 
applying different concepts. The Berkeley group (Gupta 
et al., 2007) prepared an ultracold gas of up to 10^ atoms 
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FIG. 21 (Color online) Different experimental schemes for 
preparing ultracold atoms and Bose-Einstein condensates in- 
side high-finesse optical Fabry-Perot resonators, (a) Ultracold 
atoms are prepared in a magnetic trap, formed using electro- 
magnets coaxial with the vertically oriented high-finesse cav- 
ity (length = 194 pm) and delivered along the a;-axis towards 
the cavity center. Once overlapping with the cavity mode, 
the atoms are loaded into a deep intracavity lattice poten- 
tial provided by a far-detuned cavity pump field. Adapted 
from Murch et al. (2008). (b) Ultracold atoms, prepared 
in a magnetic trap placed above the optical resonator, are 
loaded into a vertically oriented optical lattice potential and 
transported into the cavity by controlled frequency chirp- 
ing the counter-propagating laser beams. Once in the cav- 
ity (length = 176 /xm), the atoms are loaded into a crossed- 
beam harmonic dipole trap where Bose-Einstein condensa- 
tion is achieved. From Brcnnecke et al. (2007). (c) A Bose- 
Einstein condensate is prepared in an atom-chip-based mag- 
netic trap and positioned afterwards with subwavelength pre- 
cision in the mode of a fibre-based Farby-Perot cavity which 
has a length of 39 ^m. Figure courtesy of J. Reichel. (d) 
A magnetic QUIC trap is used to prepare and to transfer a 
Bose-Einstein condensate into the field of a vertically oriented 
Fabry-Perot resonator with a length of 5 cm. Figure courtesy 
of A. Hemmerich. 



within an ultra-high finesse Fabry-Perot resonator by 
loading it into a vertically oriented, deep intracavity op- 
tical lattice potential, see Fig. 21a. A cavity with similar 
parameters was used in the approach of the Ziirich group 
(Brcnnecke et al, 2007). Here, Bose-Einstein conden- 
sates of typically 2 x 10^ atoms were transported into 
the mode volume of the optical cavity using an optical 
elevator formed by two counter-propagating laser beams 
with controlled frequency difference, see Fig. 21b. The 
Paris group (Colombe et al., 2007) used an atomic chip to 
produce Bose-Einstein condensates of up to 3000 atoms 
and control its position on a sub-wavelength scales within 
a novel type of fibre-based Fabry-Perot cavity with high 
mirror curvature and reduced mode volume, see Fig. 21c. 
Also the Berkeley group (Purdy et al., 2010) achieved 
sub-wavelength positioning of Bose-Einstein condensates 
of a few thousands of atoms inside a conventional small- 
volume high-finesse optical cavity using an atomic chip. 

A novel BEC-cavity system operating in an interesting 
and so far unexplored parameter regime was presented re- 
cently by the Hamburg group (Wolko et at, 2012). Here, 
Bose-Einstein condensates of typical 2 x 10^ atoms are 
prepared and transferred magnetically into the field of a 
5 cm long near-concentric Fabry-Perot resonator result- 
ing in a large single-atom cooperative C ^ 1 and a very 
narrow cavity bandwidth on the order of the recoil fre- 
quency lor, see Fig. 21d. 

The lowest electronic excitation spectrum of degen- 
erate and non-degenerate atomic samples strongly cou- 
pled to the cavity field has been studied by Brcnnecke 
et al. (2007) and Colombo et al. (2007). The presence 
of N atoms which collectively couple to the cavity field 
results in an enhanced collective coupling which scales 
as ^/N. a correspondingly large vacuum Rabi-splitting 
was measured in the experiments (Brcnnecke et al. , 2007; 
Colombe et al., 2007). The energy spectrum of a coupled 
BEC-cavity system together with the square-root depen- 
dence of the energy splitting on the atom number are 
shown in Fig. 22. 

The electronic excitation spectrum is sensitive to the 
effective number of atoms coupled maximally to the cav- 
ity mode, i.e. it depends on the density distribution of 
the atoms integrated over the cavity mode profile. How- 
ever, the above described measurements of the electronic 
excitation spectrum have an energy resolution given by 
the excited state and cavity lifetimes, which is too large 
to probe the low-energy excitations of the external degree 
of freedom of a Bose-Einstein condensate. Probing quan- 
tum statistics and quantum correlations in atomic many- 
body states using the dispersive coupling to far-detuned 
laser and cavity fields is discussed in Sect. IV. E. 2. 
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FIG. 22 (Color online) Collective vacuum Rabi-splitting of a 
coupled BEC-cavity system. The displayed data was obtained 
by cavity transmission spectroscopy using a weak probe laser 
beam (Brennecke et ai, 2007). The detuning of the probe 
beam with respect to the frequency uja of the |-F = 1) — >■ 
\F' = 1) transition of the D2 line of *^Rb is denoted by Ap. 
Two orthogonal circular polarizations of the transmitted light 
were recorded and are displayed as red circles {a'^) and black 
triangles (o"~). (a) Position of the probed resonances as a 
function of the detuning Ac between the empty cavity reso- 
nance and the atomic transition frequency uja for 2.2 x 10"" 
atoms. Bare atomic resonances are shown as dotted lines, 
whereas the empty cavity resonance of the TEMqo is plot- 
ted as a dashed-dotted line. The solid lines are the result 
of a theoretical model including the influence of higher-order 
cavity modes, (b) Shift of the lower resonance of the coupled 
BEC-cavity system from the bare atomic resonance as a func- 
tion of atom number for Ac = 0. The solid lines are fits of 
the square-root dependence on the atom number A'^. Adapted 
from Brennecke et al. (2007). 



B. Theoretical description 

This section provides the theoretical basis for a quan- 
tum many-body description of a coupled and laser-driven 
BEC-cavity system at zero temperature. For simphcity, 
we consider two different pump laser fields with equal 
frequency uj which propagate along and transverse to the 
axis of a Fabry-Perot resonator. A similar many-body de- 
scription for the case of a BEC in a driven ring-cavity has 
been presented e.g. in Moore et al. (1999). The Hamilto- 
nian is composed of an atomic, a cavity and an atom-field 
interaction part 



H = Ha + Hc + Hac ■ 



(49) 



In the following, we will assume a sufficiently large de- 
tuning of the cavity frequency ojc and the pump laser 
frequency w from the atomic transition frequency uja, so 
that the atom-field interaction is of purely dispersive na- 



ture, see Sec. II. A. 2. In this case all excited states can 
be adiabatically eliminated and the atom resides most of 
the time in its electronic ground state. Correspondingly, 
the motional degree of freedom is captured by a scalar 
matter-wave field operator ^'(r). 

The atomic many-body Hamiltonian is given by 



Ha = j d^v^ir) \n^^^ + ^*^(r)^'(r)] ^'(r) , 



(50) 



where u — ATrh^as/m denotes the strength of the 
short-range s-wave collisions with scattering length 
(Pitacvskii and Stringari, 2003). Here, the single-atom 
Hamiltonian 



2m 



+ Ki(r) 



(51) 



includes an external trapping potential 14i(r) which also 
incorporates the potential caused by the transverse pump 
laser field. 

The dynamics of a single, coherently laser-driven cav- 
ity mode with mode function cos(fca;) and resonance fre- 
quency ujc is described in a frame rotating at the pump 
laser frequency lj by the Hamiltonian 



He ~ —hAc a' a + ifirj (a^ — a) . 



(52) 



As before, the detuning between the pump laser fre- 
quency and the cavity resonance frequency is denoted by 
Ac — UJ — uJc- The generalization to multimode cavities 
will be discussed in Sec. IV. D. 3. 

The dispersive interaction between the pump and cav- 
ity radiation fields and the atoms reads (in the frame 
rotating at uj) 

Hac = j d^r*^(r) Wq cos^ {kx) a) a 

+ hrjcS cos(fca:) cos(fcz) (a^ + a) 4'(r) . (53) 

The first term arises from the absorption and stimulated 

2 

emission of cavity photons, with Uq = denoting the 
maximum atomic light-shift for a single intracavity pho- 
ton. As before, g denotes the maximum atom-cavity 
coupling strength and Aa — uj — uja the pump-atom 
detuning. The second term corresponds to the coher- 
ent redistribution of photons between the standing- wave 
transverse pump laser (with mode function cos(fc2:)) and 
the cavity field. The maximum scattering rate for a sin- 
gle atom is given by the two-photon (vacuum) Rabi fre- 
quency 77efi = , where fl is the Rabi frequency of the 
transverse pump laser. Both interaction terms in Eq. (53) 
can be viewed as a four- wave mixing of light and matter 
wave fields (Rolston and Phillips, 2002). 

The system is subject to dissipation due to photon 
leakage through the cavity mirrors. The corresponding 
irreversible evolution can be modeled by the Liouville 
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terms in the master equation, Eq. (5a), or, equivalently, 
by a Heisenberg-Langevin equation (Gardiner and Zollcr, 
2004) for the cavity field operator a. 



d 

—I 
dt 



'i[a,H] 



(54) 



with cavity field decay rate k} The Gaussian noise oper- 
ator ^ maintains the commutation relation for the photon 
operators in the presence of cavity decay. In the optical 
domain, the temperature of the bath of electromagnetic 
field modes can be set to zero. Correspondingly, ^ has 
zero mean value and the only non-vanishing correlation 
function reads 



(55) 



according to the fluctuation-dissipation theorem. Fur- 
ther possible dissipation channels can, for example, act 
directly on the atomic cloud. 

As the cavity field mediates a global coupling among 
all atoms, a mean-field approach is well suited to solve 
the above set of equations (Horak et ai, 2000; Horak and 
Ritsch, 2001a; Nagy et ai, 2008). The mean-field descrip- 
tion assumes the presence of a macroscopically populated 
matter wave field — (vE'(r)) (condensate wavefunc- 
tion) and a coherent cavity field with amplitude a = (a) 
which can be separated from the quantum fiuctuations 
according to 



a — > a + (5a , 
*(r) ^/N~cip{r) + (5*(r) . 



(56a) 
(56b) 



Here, Nc denotes the number of condensate atoms with 
(/3(r) being normalized to 1. Quantum fluctuations are 
assumed to be small and their mean values vanish by 
definition, i.e. {5a) = and ((5^'(r)) = 0. The dynamical 
equations of motion resulting from Eqs. (56) contain a 
hierarchy of terms according to the different powers of 
the fluctuation operators. To zeroth order in the fluctu- 
ations, one obtains a Gross-Pitaevskii-type equation for 
the condensate wavefunction iy3(r,t) coupled to an ordi- 
nary differential equation for a{t): 



.9 , , 



' 2m ^ 
nUQ\a{t)\'^ cos^{kx) 



yo(r) + Afc«|^(r,t)p 
P cos^(fca;) 

2 ?i?7cfrRe{(3!(i)} cos(fca;) cos(fcz) if{Y,t) . (57a) 



^ It is important to note that the original problem is intrinsically 
time-dependent because of the external laser driving, although 
this time dependence has been formally eliminated by going into 
a rotating frame. Nevertheless, the coupling to the reservoir 
occurs at the high, optical frequency range and thus the simple 
form of the loss description can be used irrespective of the low- 
frequency dynamics imposed by the effective Hamiltonian. 



i — a(t) = [-Ac + NJJQ{cos^{kx)) - in] a{t) -f ijy 

+ Ncrjcs{cos{kx) cos(fcz)) , (57b) 

where we used the notation (/(r)) — J d^rf{r)\ip{r,t)\'^. 
The Gross-Pitaevskii equation contains potential-like 
terms which depend on the amplitude a and intensity 
jofp of the cavity field, and express the mechanical effect 
of the cavity light upon the atoms. The dynamics of the 
cavity field involves spatial averages over the condensate 
density distribution. Because of cavity decay, the time 
evolution leads to a self-consistent stationary solution for 
the mean fields, which usually is obtained only numeri- 
cally (Nagy et ai, 2008). 

For a given condensate wave function and coherent 
cavity field amplitude in steady state, the quantum fluc- 
tuations to leading order form a linear system, which 
provides the energy spectrum of excitations. With the 
notation R = [5a, (5a^, 54'(r), (5^'^'(r)], the time evolution 
of the fluctuation operators takes the compact form 



d_ 
dt 



R = MR + H , 



(58) 



where M is the linear stability matrix of the mean field 
solution (Nagy ei a/., 2008), and the term H = 0,0] 
accounts for the quantum input noise of the cavity field. 
In general, the matrix M is non-normal, i.e. it does not 
commute with its hermitian adjoint. Therefore it has 
different left and right eigenvectors, denoted by l^'^^ and 
r^'^^ that form a bi-orthogonal system with scalar prod- 
uct (l(fe)^ !■(')) = 5k,i. The decoupled quasi-normal exci- 
tation modes defined by pk = (I'-'^^R) are mixed excita- 
tions of the photon and the matter wave fields. 

The spectrum of excitations was analyzed first from a 
cavity cooling point of view in the cavity pumping geom- 
etry (77 7^ 0, rjcff — 0). The imaginary part of the spec- 
trum revealed that excitations of the ultracold atomic gas 
can be damped through the cavity loss channel (Gardiner 
et at, 2001; Horak and Ritsch, 2001a), provided the decay 
rate n is on the order of the recoil frequency w/j. The ex- 
citation spectrum was used in further studies to describe 
critical phenomena, such as the dispersive optical bista- 
bility in the cavity pumping geometry (Szirmai et al., 
2010) and the self-organization phase transition (Konya 
et ai, 2011; Nagy et ai, 2008) in the atom pumping ge- 
ometry {r]cs 7^ 0, T] = 0), see Sec. IV. C and Sec. IV. D. 

In the stable regime, second-order correlation func- 
tions can be derived from Eq. (58). Importantly, there 
can be a non-vanishing population of the atomic excited 
modes, ((5\E'^(r)(5^(r)) ^ 0, even at zero temperature. 
This quantum depletion of the condensate is indepen- 
dent of collisional interactions, which are known to cause 
a finite population of the Bogoliubov modes (Pitaevskii 
and Stringari, 2003). At variance, here the condensate 
depletion arises from the cavity-mediated atom-atom in- 
teractions as well as the dissipative process associated 
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with cavity decay. The quantum noise accompanying the 
photon loss process couples into the atomic system and 
excites atoms out of the condensate mode. Formally, it 
stems from the term containing the photon creation oper- 
ator a) in the equations of motion of ^'(r). This noise am- 
plification mechanism is analogous to the Petermann ex- 
cess noise factor in lasers with unstable cavities (Grang- 
ier and Poizat, 1998). It was shown by Szirmai et al. 
(2009) that a depletion on the order of -I- n"^ /ujp; 

can be expected rather independently of the atom-field 
coupling, even for Uq ^ Q and rjes = 0. It is a signa- 
ture of the global coupling that the quantum depletion 
is independent of the total atom number N . In most 
of the experiments with linear cavities, the ratio h/lou 
is large (^ 10'^). Since cavity decay can also be inter- 
preted as a continuous weak measurement of the cavity 
photon number, the depletion can be attributed to quan- 
tum back-action upon the atomic many-body state, as 
discussed in Sec. IV. C. 3. It is also interesting to note 
that the second-order correlation functions reveal an en- 
tanglement between the matter-wave and the cavity field 
modes (Szirmai et al., 2010). 



C. Cavity opto-mechanics with ultracold atomic ensembles 

In this section we focus on the dispersive interaction 
between the collective motion of a quantum gas and 
a single-mode Fabry-Perot cavity, which is coherently 
driven with amplitude 77 by a laser field at frequency uj. 
In this case the dispersive matter- light interaction, (53), 
is given by 



Hac = j Sv^^y) HUq cos^ (kx) a ^{r) 



(59) 



On the one hand, the atomic medium experiences a pe- 
riodic potential, whose depth is proportional to the in- 
tracavity photon number a^a. The potential depth for a 
single cavity photon is L^o = g'^ /Aa and can be tuned in 
the experiment via the detuning between the cavity 
pump frequency uj and the atomic transition frequency. 
On the other hand, the atom-light interaction causes a 
dispersive shift of the empty cavity frequency, which is 
determined by the spatial overlap between the atomic 
density \E''''(r)^'(r) and the cavity intensity mode func- 
tion cos^(A;a;). A change in the atomic density distribu- 
tion caused by the intracavity dipole force can therefore 
dynamically act back on the intracavity field intensity by 
shifting the cavity resonance with respect to the driving 
field. 

In general, the interplay of these two effects results in 
a highly nonlinear evolution of the coupled atoms-cavity 
system. For certain limiting situations, however, the sys- 
tem can effectively be described in the framework of cav- 
ity optomechanics (Kippcnbcrg and Vahala, 2008), which 



studies the radiation-pressure interaction between a har- 
monically suspended mechanical element and the field 
inside an electromagnetic resonator. In a frame rotating 
at Lu this is described by the generic cavity optomechanics 
Hamiltonian 



H, 



OM 



hojmC^ c — h{Sc — GX)a^a + ihr]{a^ — a) (60) 



where and c denote creation and annihilation opera- 
tors of the mechanical oscillator at frequency uJm- The 
mechanical element couples via its position quadrature 
X = {c + c^)/\/2 with coupling strength G to the in- 
tracavity photon number a^a. The detuning between 
the driving laser and the cavity resonance frequency for 
zero displacement X is denoted by Sc- Configurations 
in which the harmonic oscillator couples quadratically in 
X to the cavity field have recently been realized (Purdy 
et al, 2010; Thompson et al, 2008). This offers the pos- 
sibility to detect phonon Fock states of the mechanical 
element and to prepare squeezed states of the mechanical 
oscillator or the optical output field. 



1. Experimental realizations 

Particular experimental situations allow to realize the 
optomechanics Hamiltonian Eq. (60) with an atomic en- 
semble dispersively coupled to the field inside an opti- 
cal cavity. This relies on the fact, that the cavity field 
affects and senses predominantly a single collective mo- 
tional mode, which matches the spatial cavity mode pro- 
file and plays the role of the harmonically suspended me- 
chanical element. Two different approaches for realizing 
cavity optomechanics with ultracold atoms have so far 
been realized experimentally. 



a. Collective center-of-mass motion in the Lamb-Dicke regime 
In experiments performed by the Berkeley group (Gupta 
et at, 2007; Murcli et at, 2008), ultracold atoms are 
loaded into the lowest band of a far-detuned intracavity 
lattice potential, forming a stack of hundreds of tightly 
confined atom clouds (see Fig. 23). Each atom cloud 
is harmonically suspended with oscillation frequency 
and extends along the cavity axis by only a fraction of 
the optical wavelength, thus realizing the Lamb-Dicke 
regime. A cavity mode, whose periodicity differs from 
that of the trapping lattice potential, couples strongly 
to a single collective center-of-mass mode of the atomic 
stack. All remaining collective modes decouple from the 
cavity field and can be considered as a heat bath to which 
the distinguished collective mode is only weakly coupled 
via e.g. collisional atom-atom interactions. The system 
realizes the linear optomechanics Hamiltonian Eq. (60), 
with the optomechanical coupling strength being given 
by G = ^/NcskUoXho- Here, k is the cavity wave vec- 
tor, A"ho = \/^/ '2muJm denotes the harmonic oscillator 
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FIG. 23 (Color online) Scheme for cavity optomechanics with 
ultracold atoms confined in the Lamb-Dicke regime. A high- 
finesse cavity supports two longitudinal modes: one with 
wavelength of about 780 nm that is near the D2 line of **^Rb, 
and another with wavelength of about 850 nm. The latter pro- 
duces a one-dimensional optical lattice potential, with trap 
minima indicated in orange, in which ultracold *^Rb atoms 
are confined within the lowest vibrational band. The atomic 
clouds induce, depending on their trapping position Zi, disper- 
sive frequency shifts on the 780 nm cavity resonance, as shown 
in the bottom line. In turn, the cavity field exerts a position 
dependent force /, as indicated by the red arrows. In the 
Lamb-Dicke regime, the collective atoms-cavity interaction 
reduces to the generic optomechanics Hamiltonian wherein 
a single collective mode of harmonic motion linearly couples 
to the cavity field. From Better et al. (2009). 



length with atomic mass m, and iVoff' ~ N/2 with total 
atom number N. 

The quadratic coupling regime of optomechanics with 
ultracold atoms was realized in an atom-chip-based setup 
(Purdy et ai, 2010) which allows for subwavelength po- 
sitioning of a tightly confined ultracold ensemble. By 
preparing as low as two atomic clouds, tightly confined 
at adjacent lattice sites of a far-detuned intracavity lat- 
tice potential, and controlling their center-of-mass po- 
sition along the cavity axis both linear and quadratic 
optomechanical coupling can be realized, providing an 
atoms-based realization of the " membran-in-the-middle" 
approach (Thompson et ai, 2008). 



b. Collective density oscillations in a Bose-Einstein conden- 
sate A different route to realize linear cavity optome- 
chanics with an ultracold atomic ensemble was experi- 
mentally explored in the Ziirich group (Brcnncckc et ai, 
2008). A Bose-Einstein condensate of typically 10^ atoms 
is prepared in an external harmonic trapping potential, 
extending over several periods of the cavity standing- 
wave mode structure (see Fig. 24). In contrast to the 
Lamb-Dicke regime considered before, here a momen- 
tum picture is more appropriate. Initially, all condensate 
atoms are prepared - relative to the recoil momentum 
hk — in the zero-momentum state \p = 0). The dis- 
persive interaction with the cavity field diffracts atoms 
into the symmetric superposition of momentum states 
I ± 2hk) along the cavity axis. Matter-wave interference 
between the macroscopically occupied zero-momentum 



component and the recoiling component results in a 
spatial modulation of the condensate density with pe- 
riodicity A/2, which oscillates in time at the frequency 
4wfl = 2hk^/m. As long as diffraction into higher-order 
momentum modes can be neglected, the dynamics of the 
coupled system is again captured by the simple optome- 
chanics Hamiltonian Eq. (60). Here, collective excita- 
tions of the recoiling momentum mode play the role of 
phonon excitations of a mechanical mode with oscillation 
frequency = 4:L0r. The coupling rate, G = VNUq/2, 
again scales with the square root of the atoms number, 
indicating the collective nature of the atom-light interac- 
tion. The realization of an optomechanical system em- 
ploying the collective motion of a Bose-Einstein conden- 
sate triggered subsequent theoretical studies along this 
direction (Chen et ai, 2011, 2010; De Chiara et ai, 2011; 
Zhang et ai, 2009). 
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FIG. 24 (Color online) Cavity-optomechanics with a weakly 
confined Bose-Einstein condensate dispersively coupled to the 
field of an optical high-finesse resonator, (a) A collective 
density excitation of the condensate (blue) with periodicity 
A/2 = Tr/k acts as a mechanical oscillator with oscillation 
frequency Alur. Optomechanical coupling to the cavity field 
is provided by the dependence of the optical path length on 
the atomic density distribution within the spatially periodic 
cavity mode structure, (b) Condensate atoms initially pre- 
pared close to zero momentum, p = 0, are scattered ofi' the 
intracavity lattice potential into the symmetric superposition 
of states with momentum p = ±2hk. Matter-wave interfer- 
ence with the macroscopic zero-momentum component results 
in a harmonic density oscillation evolving at frequency Aujji. 
Adapted from Brennecke et al. (2008). 



Further realizations of cavity optomechanics with 
atomic ensembles have been proposed theoretically by 
dispersively coupling a quantum-degenerate Fermi gas 
(Kanamoto and Meystrc, 2010) or the internal spin- 
degrees of freedom of a quantum gas (Brahms and Kurn, 
2010; Jing et ai, 2011) to the field of an optical cavity. 
The latter system has been shown to exhibit a formal 
analogy with a torsional oscillator coupled quadratically 
to the cavity mode. It provides an ideal nondestructive 
tool for the control of quantum spin dynamics, and was 
proposed to resolve the quantum regime of an antiferro- 
magnetic spin-1 condensate. 
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The realization of an optomechanical system using ul- 
tracold atoms offers direct access to the quantum regime 
of cavity optomechanics. In contrast to solid-state real- 
izations of optomechanics, evaporative cooling techniques 
available for atomic gases allow for a natural and very 
pure preparation of the mechanical oscillator mode in its 
quantum ground state. Correspondingly, these systems 
pave the way to directly study quantum effects of the 
optomechanical interaction (Brahms et al., 2012; Brooks 
et al, 2012; Murch et al, 2008). 

The easy tunability of system parameters like e.g. the 
mechanical oscillator frequency (via the external con- 
fining potential), the optomechanical coupling strength 
G (via the atom number or the pump-atom detuning) or 
the initial temperature of the mechanical oscillator allows 
to explore the transition between different regimes of op- 
tomechanics. Most important, the coupling strengths G 
achievable with atomic systems open access to the 'granu- 
lar' regime of optomechanics (Ludwig et al., 2008; Murch 
et al, 2008), where single excitations in either of the two 
subsystems have a non-negligible effect upon the dynam- 
ics of the other. This can be measured by the granularity 
(or quantum) parameter which is defined as C = G/k. 
For C = already a single excitation of the mechanical 
mode shifts the cavity resonance by half its linewidth, 
and already a single photon entering the cavity imparts 
one excitation quantum in the mechanical oscillator. In 
future research, this might allow the generation and de- 
tection of quantum correlations between the mechanical 
and light degrees of freedom. Further research possibili- 
ties based on atoms-based realizations of optomechanics, 
are given by possible implementations of precision sen- 
sors of forces, or devices to manipulate light fields on a 
quantum level. 



2. Nonlinear dynamics and bistability for low photon number 

The optomechanical interaction, Eq. (60), being intrin- 
sically nonlinear gives rise to dispersive optical bistability 
and nonlinear dynamics of the coupled system. Optical 
bistability (Lugiato, 1984), a well studied phenomenon in 
nonlinear optics, refers to the co-existence of two stable 
steady-state solutions when e.g. driving an optical cavity 
filled with a medium whose refractive index depends on 
the light intensity. In typical nonlinear Kerr media and 
solid-state realizations of optomechanics, the occurrence 
of bistability typically requires large intracavity power in 
order to significantly alter the system's optical proper- 
ties. The large coupling strength achieved in the atomic- 
ensemble realizations of optomechanics induces optical 
bistability at a mean-intracavity photon level below one 
(Gupta et al, 2007; Ritter et al, 2009). This achieve- 
ment is desirable for applications ranging from optical 
communication to quantum computation (Cirac et al., 
1997; Imamoglu et al, 1997). 



The occurrence of bistability in optomechanical sys- 
tems can be understood from the corresponding semiclas- 
sical equations of motion for the oscillator displacement 
X and the coherent intracavity field amplitude a derived 
from Hamiltonian Eq. (60) 



X + UJmX = -UJ„iG\c 



(61) 



Q: = {i{Sc ^ GX) — k) a + 1] . 



In the bad cavity regime, k 3> lor, the atoms move on a 
timescale which is large compared to the lifetime {2k)~^ 
of intracavity photons. Correspondingly, the cavity field 
adiabatically follows the atomic dynamics according to 



(62) 



+ {5c - GXy ■ 



Retardation effects resulting in dynamical back-action 
cooling or heating of the mechanical element are ne- 
glected in this approximation. Inserting this expres- 
sion into the equation of motion for X, Eq. (61), yields 
X = —-^VoMiX). The optomechanical potential given 
by 



I^om(^) 



1 



-hU-mX^ 



^arctanf^V (63) 

K \ K J 



captures the combination of the harmonic confinement 
and the cavity dipole forces. Here, A{X) — 5c — GX 
denotes the detuning between the driving laser and the 
atoms-shifted cavity resonance. 
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FIG. 25 (Color online) Optomechanical potential and bista- 
bility. (a-d) Optomechanical potential landscape VbM(^) for 
different pump-cavity detunings 5c, indicated by the dashed 
lines in (e). The shaded regions show the resonance profile 
of the cavity, (e) Mean intracavity photon number lap in 
steady state. Open and close circles correspond to the situa- 
tion shown in (b). Parameters are G = 0.42k and rj — -v/S'ycr. 



The optomechanical potential provides an intuitive 
picture to understand the steady-state as well as the 
dynamical behavior of the coupled system, see Fig. 25. 
Above a critical cavity pump strength r^cri determined 
by rjl^ — ^^^ '^g^ , the optomechanical potential ex- 
hibits (within a certain detuning range) two local min- 
ima, which correspond to different intracavity light inten- 
sities as shown in the bistable resonance curve. Fig. 25e. 
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Depending on the direction into which 5c is adiabati- 
cally tuned, the system remains in either of the two local 
minimum of VoM, following the upper or lower bistable 
resonance branch. When reaching the critical detuning, 
where one of the local minima turns into a saddle point, 
the system starts to perform transient oscillations in the 
remaining potential minimum, which translate into a pe- 
riodically modulated cavity light intensity. Due to damp- 
ing of the collective atomic motion, the system finally re- 
laxes to the steady-state in the remaining potential min- 
imum. 

Optical bistability induced by collective atomic motion 
was observed at low intracavity photon number both in 
the Berkeley group (Gupta et ai, 2007) and the Ziirich 
group (Ritter et ai, 2009). The lower and upper bistabil- 
ity branches were observed in single experimental runs by 
slowly sweeping the frequency of the driving laser twice 
across resonance, first with increasing detuning and then 
with decreasing detuning, see Fig. 26. Upon increasing 
the probe strength, the cavity transmission profile be- 
comes more and more asymmetric and exhibits hystere- 
sis. 




Apc/27t (MHz) Apc/2jt (MHz) 

FIG. 26 (Color online) Dispersive optical bistability with col- 
lective atomic motion, (a) Observed cavity line shapes (red) 
for increasing cavity input power at low intracavity photon 
number n and model line shapes (black), based on the Voigt- 
profile of the bare cavity line (inset) . Ape denotes the detun- 
ing between the probe laser frequency and the empty cavity 
frequency, (b) Lower (blue) and upper (red) branches of op- 
tical bistability observed in single sweeps across resonance. 
From Gupta et al. (2007). 



Dispersive optical bistability with collective atomic 
motion was also studied in the regime of quadratic op- 
tomechanical coupling (Purdy et ai, 2010). Here, instead 
of displacing the center-of-mass motion of the mechani- 
cal element, the intracavity dipole forces increase or de- 
crease the rms width of the compressible atomic ensem- 
ble, depending on whether the atoms are confined at a 
maximum or a minimum of the intracavity probe lattice 
potential. The corresponding change of the dispersive 
cavity shift leads again to bistable resonance curves as 
observed in the experiment (Purdy et al, 2010). 

Dynamical optomechanical effects arise in small- 
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FIG. 27 Nonlinear dynamics of a driven BEC-cavity system, 
(a) Sweeping the driving laser frequency across the bistable 
resonance curve (indicated in dashed, scaled by a factor of 
4) excites large-amplitude density oscillations in the conden- 
sate, (b) The magnified cavity transmission signal indicates 
how the density oscillations tune the cavity frequency peri- 
odically in an out of resonance with the driving laser. The 
mean intracavity photon number on resonance was 7.3 corre- 
sponding to a photon count rate of 5.8 MHz. Adapted from 
Brennecke et al. (2008). 



amplitude oscillations of the system around steady state. 
As a result of the optomechanical interaction the fre- 
quency of such oscillations is shifted with respect to the 
bare oscillation frequency in the literature often re- 
ferred to as the "optical spring effect". In the case of 
linear optomechanical coupling this can be inferred from 
a quadratic expansion around the steady state minima of 
the optomechanical potential Vqmj Fig- 25. Experimen- 
tally, the optomechanical frequency shift for the collective 
atomic motion was observed and quantified in agreement 
with theory both in the linear and the quadratic coupling 
regime (Purdy et ai, 2010). Highly nonlinear oscillations 
in the optomechanical potential with relatively large am- 
plitude have been excited and observed in cavity trans- 
mission either by a sudden displacement of the optome- 
chanical potential (Gupta et ai, 2007) or by crossing the 
instability point of the bistable curves (Brennecke et ai, 
2008), see Fig. 27. 



3. Quantum-measurement back-action upon collective atomic 
motion 

The accuracy of any position measurement of a me- 
chanical element is limited by quantum mechanics. Re- 
ferred to as the standard quantum limit, this has been ex- 
tensively studied in connection with the development of 
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gravitational- wave detectors (Caves, 1980). In a generic 
optomechanical setup, which allows for high-precision 
measurements of the position of a mechanical element, 
the standard quantum limit arises from the balance be- 
tween two noise terms: (i) detection shot noise, given 
by the random arrivals of photons on the detector and 
(ii) radiation-pressure induced displacement noise caused 
by the quantum fluctuations of the intracavity photon 
number. Whereas detection noise can be decreased by 
increasing the light power, this comes at the expense 
of increased radiation pressure force fluctuations. The 
optimal sensitivity is achieved if these two noise contri- 
bution are balanced. A direct experimental observation 
of the intracavity photon number fluctuations, requires 
large optomechanical coupling strengths between the in- 
tracavity field and the mechanical element in combina- 
tion with the suppression of thermal or technical noise 
sources which perturb the mechanical motion. 

The utilization of collective atomic motion of an ul- 
tracold gas strongly coupled to the field inside a Fabry- 
Perot resonator, allowed for the first observation of 
measurement-induced back-action upon a macroscopic 
mechanical element formed of 10^ atoms, caused by intra- 
cavity quantum force fluctuations (Murch et at, 2008). 
In the non-granular regime ^ = G/k <^ 1, the spectral 
density of intracavity photon number fluctuations (Mar- 
quardt et ai, 2007; Nagy et ai, 2009) agrees with that 
in an empty driven cavity, and reads 



2nK, 



k2 + (A(X)+w)2 



(64) 



Here, n = jap denotes the steady-state mean-intracavity 
photon number given in Eq. (62). These photon num- 
ber fluctuations are transmitted into the momentum of 
the mechanical element via the optomechanical interac- 
tion, giving rise to a diffusion-like increase of the phonon 
number 



at 



(65) 



as can be derived e.g. from an effective master equation 
for the mechanical oscillator (Nagy et ai, 2009). 

Murch et al. (2008) measured the corresponding heat- 
ing rate of the atomic ensemble in a bolometric way by 
quantifying the evaporative atom loss, see Fig. 28. After 
preparing the mechanical oscillator close to its ground 
state, the cavity transmission of a weak probe beam at 
fixed frequency is monitored on a single-photon count- 
ing module. Continuous background atom loss tunes 
the atoms-shifted cavity frequency in resonance with the 
driving laser. The atom loss rate is deduced from the 
comparison between the recorded transmission curve and 
the empty-cavity resonance curve. The corresponding 
single-atom heating rate is found to exceed the free-space 
spontaneous heating rate, which was deduced by measur- 
ing the atom loss rate far from the cavity resonance, by 
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FIG. 28 Observing measurement-induced back-action upon 
the collective motion of an ultracold atomic ensemble, (a) 
Mean-intracavity photon number n (points), monitored as the 
system is brought across the cavity resonance due to evapora- 
tive atom loss. The expected photon number including (solid 
line) and excluding (dashed line) measurement back-action 
is shown, (b) Total atom number A'^ as a function of time 
as inferred from data shown in (a), using the empty cavity 
line shape and the linear relation between atom number and 
dispersive cavity shift (inset). From Murch et al. (2008). 



a factor of 40, in agreement with the theoretical expec- 
tation. As cavity-mediated coherent amplification and 
damping of the mechanical oscillator is negligible in the 
experiment, the observation of back-action heating can 
be interpreted as a direct measurement of photon num- 
ber fluctuations in a coherently driven cavity. 

Another direct signature of quantum back-action of 
light upon collective atomic motion was obtained by mon- 
itoring the Stokes and anti-Stokes sidebands of the cavity 
transmission subsequent to the preparation of the collec- 
tive motional degree of freedom close to its ground state 
(Brahms et a/., 2012), see Fig. 29. The observed side- 
band asymmetry provides a direct measurement of the 
quantized collective motion and serves as a record of the 
energy exchanged between motion and the light in agree- 
ment with a continuous back-action limited quantum po- 
sition measurement. 

The disturbance of collective atomic motion via the 
intracavity quantum force fluctuations acts back again 
onto the intracavity light field. In particular, the result- 
ing motional-induced modulation of the cavity field can 
interfere with the coherent or vacuum cavity input field 
giving rise to nonlinear optical parametric amplification 
and - for negligible technical or thermal fluctuations - 
to ponderomotive squeezing (Fabre et ai, 1994; Mancini 
and Tombesi, 1994). Only recently, these effects where 
observed for the first time in the Berkeley group utilizing 
the optomechanical coupling between collective atomic 
motion and an optical cavity field (Brooks et ai, 2012). 
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these processes which is achieved by driving the cavity 
near the anti-Stokes sideband results in a steady-state 
phonon occupation number of (5^)^ ^ 1 (Marquardt 
et ai, 2007). Optomechanical cooling in the resolved 
sideband regime is equivalent to optical Raman sideband 
cooling of tightly confined atoms or ions. 

Cavity cooling in an optomechanical-type BEC-cavity 
system which ranges in the good cavity regime, k < 
ujm = 'iujR, was demonstrated recently by Wolke et al. 
(2012). By driving the cavity field selectively close to 
the Stokes or anti-Stokes sidebands atoms where trans- 
ferred via cavity-stimulating backward scattering from 
the macroscopically populated zero-momentum state into 
a superposition of momentum states | ± 2hk) and back, 
see Fig. 30. This experiment paves the way towards 
the achievement of quantum degeneracy starting from a 
thermal gas without relying on evaporative cooling tech- 
niques. 



FIG. 29 Optical detection of quantization of collective atomic 
motion in the cavity output spectrum. Shown are mea- 
sured Stokes sidebands (left panels) and anti-Stokes side- 
bands (right panels) for increasing intracavity photon num- 
ber (bottom to top) together with the theoretical prediction 
(solid lines). The observed Stokes asymmetry provides a 
calibration-free measure for the mean occupation number of 
the mechanical oscillator, which was deduce for the lowest 
graph to be 0.49. The mechanical oscillation frequency was 
aj„ = 27r X 110 kHz. From Brahms et al. (2012). 



<3 



FIG. 30 Observation of sub-recoil cavity cooling with a *^Rb 
BEG in a narrow-bandwidth Fabry-Perot resonator. Shown 
are atomic momentum distributions after driving the cavity 
field with a far-detuned laser field at 803 nm. First, a 400 /is 
long pulse, blue-detuned from the cavity resonance, transfers 
atoms into the momentum states | ±2Tik) (left), subsequently 
a 200 /^s long red-detuned pulse transfers the atoms back into 
the zero-momentum state. Binary atomic collisions result in 
a substantial depletion of the ±2?ifc-momentum state popu- 
lations visible as a diffusive halo. From Wolke et al. (2012). 



4. Cavity cooling in the resolved sideband regime 

For the small-volume cavities which were employed in 
the experiments performed by the Berkeley and Zurich 
group, the cavity decay rate k exceeds the mechanical 
oscillation frequency uJm of the collective atomic degree 
of freedom by more than one order of magnitude. In this 
non-resolved sideband regime of cavity optomechanics 
cooling of the mechanical oscillator into its ground state 
utilizing cavity dissipation is not possible ^. Rather, the 
minimal steady-state occupation number when driving 
the cavity field with a laser field which is red-detuned by 

a/w^ -I- k2 from the cavity resonance is given by > 1 D. Non-equilibrium phase transitions 

for weak optomechanical coupling strength (Marquardt 
et al, 2007). 

Ground-state cooling becomes possible only in the re- 
solved sideband regime where uj^ 3> n (Kippenberg and 
Vahala, 2008). Here, the cavity is able to resolve the 
Stokes resp. anti-Stokes sidebands which correspond to 
adding resp. removing motional quanta from the mechan- 
ical degree of freedom. The large asymmetry between 



^ Indeed, in those experiments the preparation of the mechanical 
degree of freedom close to its ground state was achieved directly 
by evaporative cooling. 



Self-organization of a laser-driven atomic ensemble in- 
side an optical resonator, as was considered for ther- 
mal atoms in Section III.B.l, was extended both theo- 
retically and experimentally into the ultracold regime, 
where atomic motion becomes quantized. Correspond- 
ingly, the transition point to the self-organized phase is 
not determined anymore by thermal density fluctuations, 
rather by the competition between kinetic energy cost 
and potential energy gain associated with a spatial mod- 
ulation of the atomic matter-wave in the cavity-induced 
lattice potential. In case of a weakly interacting Bose- 
Einstein condensate, the reduced number of momentum 
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states accompanying in the dynamics allows for a sim- 
plified description in terms of a collective spin degree of 
freedom, providing a direct link between self-organization 
and an open-system realization of the Dicke quantum 
phase transition. 



1. Self-organization of a Bose-Einstein condensate 

Self-organization of a dilute Bose-Einstein condensate 
(BEC), which is located in a single- mode optical cav- 
ity and illuminated transversally to the cavity axis by a 
far-detuned laser field, was studied theoretically by Nagy 
et al. (2008). In terms of a mean-field description, the 
steady-state of the system was obtained from the equa- 
tions of motion for the coherent cavity field amplitude a 
and the atomic mean- field </5(r), see Eq. (57b) and (57a), 
setting the on-axis pump strength 77 to zero. For simplic- 
ity only atomic motion along the cavity axis was taken 
into account. The numerical solution for the steady-state 
order parameter Q = {ip\ cos{kx)\(p), obtained by numeri- 
cally propagating the equations of motion into imaginary 
time, is shown in Fig. 31. Above a critical two-photon 
Rabi frequency r^eff, the order parameter takes a non- 
zero value indicating self-organization of the atoms in a 
A-periodic density pattern. A stability analysis of the 
non-organized steady state, = 0, yields the following 
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FIG. 31 Self-organization of a driven Bose-Einstein con- 
densate in a standing-wave cavity. Plotted is the steady- 
state order parameter as a function of the effective cavity 
pump strength r^cff , as obtained from a numerical solution of 
the mean-field equations. Parameters are NUo = — lOOojjj, 
Ac = —SOOujR, K = 200aJi{ and fio = lOhujR in the homoge- 
neous phase. According to Eq. (66) the homogeneous phase is 
stabilized in this parameter regime dominantly by coUisional 
interaction energy. The inset shows the condensate wave func- 
tions (solid lines) for y/N rjcs = IOOivr (broader, brown solid) 
and 300a;_R (narrower, blue solid) and the corresponding opti- 
cal dipole potentials (dashed lines). Adapted from Nagy et al. 
(2008). 
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FIG. 32 Collective excitation spectrum of the transversally 
driven condensate-cavity system. Shown are the eigenfre- 
quencies of the lowest six collective atomic-like and the first 
cavity-like (divided by 5) excited states as a function of 
the transverse pump amplitude. For va nishing pump am- 
plitude the Bogoliubov spectrum Q„ = y^n'^ujii(n'^u!ii + 2/io) 
for a condensate in a ID box potential of size A is retained. 
Self-organization is indicated by the softening of the lowest 
lying collective mode towards the critical pump amplitude 
V^^cff.c ~ 65.5aj_R. Parameters are the same as in Fig. 31. 
From Nagy et al. (2008). 



analytic expression for the critical point 



-25n 



(66) 



where /io denotes the chemical potential of the homoge- 
neous condensate and 5c the detuning of the pump laser 
from the dispersively shifted cavity resonance. In con- 
trast to the thermal case Eq. (33), the critical transverse 
pump power scales in the zero-temperature limit with the 
recoil frequency (and the chemical potential), which re- 
flects the fact that the homogeneous phase is stabilized 
by the kinetic energy (and atom-atom collisions). 

A deeper understanding of the process of self- 
organization is gained from the collective excitation spec- 
trum on top of the steady-state mean-field solution. This 
was calculated in Konya et al. (2011) and Nagy et al. 
(2008) using a Bogoliubov-type approach based on the 
separation ansatz Eq. (56). The eigenvalues of the lin- 
earized equations for condensate and cavity fluctuations, 
Eq. (58), yield the energy spectrum of excitations (polari- 
tons) shown in Fig. 32. For the considered case where the 
pump-cavity detuning 5c is large compared to the recoil 
frequency lor, the excitations separate into two classes, 
according to whether they are dominantly cavity-like or 
atom-like. The occurrence of self-organization is recog- 
nized in a characteristic softening of the atom-like excita- 
tion mode which matches the spatial interference pattern 
between cavity and transverse pump mode (Fig. 32, red 
solid line). 
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FIG. 33 Observation of self-organization with a Bose-Einstein 
condensate. Simultaneous time traces of the mean-intracavity 
photon number (middle panel) and the relative pump-cavity 
phase A</> (lower panel) while ramping the transverse pump 
power twice across the critical point at « 0.35 mW. The 
absorption images (upper panel) show the atomic momen- 
tum distribution for the indicates times. The line of sight 
is perpendicular to the pump-cavity plane. Parameters are 
Ac = -27r X 20 MHz, k = 27r x 1.3 MHz and N = 10^ 
Adapted from Baumann et al. (2010, 2011). 



Self-organization of a Bose-Einstein condensate was 
observed by the Zurich group (Baumann et al., 2010). 
A BEC of about 10^ atoms, harmonically confined inside 
a high-finesse optical Fabry-Perot resonator, was illumi- 
nated by a far red-detuned standing-wave laser beam. 
By gradually increasing the power of the transverse laser 
beam, the transition to the self-organized phase was ob- 
served in a sharp raise of the intracavity light intensity ac- 
companied by the build-up of macroscopic populations in 
the momentum states {px,Pz) — (±?ifc, ihk), see Fig. 33. 
Above the critical pump power, the relative phase A0 
between pump field and cavity field is observed to stay 
constant, which demonstrates that the system reached a 
steady state. By controlling the transverse pump power, 
the system can be transferred repeatedly from the normal 
in the self-organized phase and back (Fig. 33). 

The process of symmetry breaking at the transition 
point was studied in Baumann et al. (2011). In re- 
peated realizations of the self-organized phase, two pos- 
sible values of the relative phase A(j) with a difference 
of TT were observed, according to self-organization into 
either the even (m(x, z) > 1) or odd (w(x, z) < 1) 
sites of the underlying mode interference mode profile 
u{x,z) = cos{kx) cos{kz) (Fig. 34). The finite spatial 
extent of the atomic cloud results in a tiny imbalance 
between the even-odd populations in the non-organized 
phase. This effectively acts as a symmetry breaking field. 




time (ms) 

FIG. 34 Observation of symmetry breaking at the self- 
organization transition with a BEC. Shown is in red the rela- 
tive pump-cavity phase (j> monitored on a heterodyne detector 
while repeatedly entering the self-organized phase by tuning 
the transverse pump power P (dashed). The system orga- 
nizes into one out of two possible checkerboard patterns cor- 
responding to the two observed phase values differing by tt. 
From Baumann et al. (2010). 



which favors the realization of one particular organized 
pattern, as was observed in the experiment. The influ- 
ence of the symmetry breaking field could be overcome 
by increasing the speed at which the transition is crossed, 
in accordance with a simple model description based on 
the adiabaticity condition (Baumann et al., 2011). 

In the limit where the cavity field adiabatically fol- 
lows the atomic motion, the process of self-organization 
can also be understood as a result of the cavity-mediated 
atom-atom interactions, see III.A.l. On a microscopic 
level, these are induced by the virtual exchange of cavity 
photons between different laser-driven atoms, accompa- 
nied by the creation of atom pairs recoiling with momen- 
tum hk along the pump and cavity direction. The result- 
ing A-periodic density correlations in the atomic cloud 
energetically compete with the cost in kinetic energy, 
which gives rise to a characteristic roton-type softening 
in the dispersion relation of the condensate at momenta 
{±.fik,zLhk), see also Fig. 32. Once the softened excita- 
tion energy reaches the ground state energy, the system 
self-organizes by macroscopically occupying those mo- 
mentum states. Such mode softening was observed by the 
Zurich group (Mottl et ai, 2012) using a variant of Bragg 
spectroscopy (Stenger et al., 1999) where the cavity field 
was probed with a weak laser pulse whose frequency was 
detuned by a variable amount from the transverse pump 
laser field. The observed excitation spectrum as a func- 
tion of sign and modulus of the cavity-mediated atom- 
atom interaction strength V is depicted in Fig. 35. The 
vanishing of the excitation gap at the transition point to- 
wards the organized phase is accompanied by a diverging 
susceptibility of the system to A-periodic density pertur- 
bations (Mottl et at, 2012). As was theoretically con- 
sidered by Oztop et al. (2012), the softened excitation 
spectrum can also be probed parametrically via ampli- 
tude modulation of the transverse pump laser. 

Conceptually, the self-organized BEC can be regarded 
as a supersolid (Gopalakrishnan et ai, 2009; Leggett, 
1970), similar to those proposed for two-component sys- 
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FIG. 35 (Color online) Observation of mode softening in- 
duced by cavity-mediated atom-atom interactions in a Bose- 
Einstein condensate. Shown is the motional atomic excitation 
energy at momenta (ikhk, ±hk) along the cavity and pump di- 
rection as a function of the transverse laser power P , which 
sets the modulus \ V\ of the cavity-mediated atom-atom inter- 
action. The sign of V is determined by the sign of 5c- For 
negative interaction strength V, the system organizes at the 
critical pump power Per, while for positive interaction an in- 
creased excitation energy is observed in accordance with the 
absence of a phase transition. From Mottl et al. (2012). 



terns (Biichler and Blatter, 2003). Non-trivial diagonal 
long-range order is induced by the cavity-mediated long- 
range interactions, which restricts the periodic density 
modulation to two possible checkerboard patterns, in 
contrast to traditional optical lattice experiments with 
laser fields propagating in free space. Simultaneously, 
the organized phase exhibits ofF-diagonal long-range or- 
der which is not destroyed while crossing the phase tran- 
sition. Only when deeply entering the organized phase, 
tunneling between different sites of the optical checker- 
board potential gets suppressed and phase coherence is 
lost (Vidal et a/., 2010). 



2. Open-system realization of the Dicke quantum phase 
transition 

Self-organization of a laser-driven BEG in an opti- 
cal resonator can be considered as an open-system re- 
alization of the Dicke quantum phase transition, where 
the quantized atomic motion acts as a macroscopic spin 
which strongly couples to the cavity field. The Dicke 
model goes back to the pioneering work of R. W. Dicke 
(Dicke, 1954) and describes the collective interaction be- 
tween matter and the electromagnetic field. Consider N 
two-level systems with transition frequency wq, forming 
a collective spin variable J, which couple identically to a 
single resonator mode at frequency w. This system can 
be described in terms of the Dicke Hamiltonian (also re- 
ferred to as the Tavis-Cummings model (Tavis and Cum- 



mings, 1968)) 

-ffDickc/?i = uja)a + u}qJz + --j={J+ + J-)(a + a)) (67) 

with the collective coupling strength denoted by A oc 
^/N . The ladder operators J± = i iJy describe cre- 
ation and annihilation of collective atomic excitations. 

According to Dicke (1954), a collectively excited 
medium, which carries correlations among the different 
atomic dipoles, decays within a much shorter time into 
its ground state than a single atom. This phenomenon, 
termed superradiance (or superfluorescence) , originates 
from spontaneous phase-locking of the different radia- 
tors resulting in a short radiation burst whose intensity 
is proportional to the number of atoms squared. Super- 
radiant emission of laser-excited media has been studied 
extensively in the past (Gross and Haroche, 1982). 

In contrast to this transient phenomenon, the Dicke 
Hamiltonian Eq. (67) was shown in 1973 to exhibit also a 
ground-state version of superradiance (Carmichael et al., 
1973; Hcpp and Licb, 1973; Lambert et al, 2004; Wang 
and Hioe, 1973). When the collective coupling strength 
A reaches the critical value Acr = v^wwo/2, the Dicke 
model undergoes a quantum phase transition from a nor- 
mal into a superradiant phase, which is characterized by a 
macroscopic cavity field amplitude (a) and a macroscopic 
polarization (J_) of the atomic medium. Apart from 
its fragility upon the inclusion of the term originat- 
ing from the minimal coupling Hamiltonian (Rzazewski 
et al, 1975), the experimental realization of the super- 
radiant Dicke phase transition with direct dipole transi- 
tions was obscured in the past due to practical limitations 
in the available dipole coupling strengths. 

The proposal by Dimcr et al. (2007) circumvents these 
issues by considering a pair of stable atomic ground states 
which are coupled via two different Raman transitions in- 
volving a single ring cavity mode and external laser fields. 
This scheme realizes the Dicke model through an effec- 
tive Hamiltonian in an open-system dynamics, including 
external driving and cavity loss, where the critical cou- 
pling strength can be reached for realistic experimental 
parameters. 

The transversally driven BEC-cavity system is for- 
mally equivalent to this proposal upon replacing the elec- 
tronic atomic states by a pair of motional atomic states, 
as was shown in Baumann et al. (2010) and Nagy et al. 
(2010). The two motional states are given by the flat 
condensate mode \px,Pz) = |0,0) and the coherent su- 
perposition of the four momentum states | ± hk,±hk), 
where x and z denote the cavity and pump direction, re- 
spectively. Coherent light scattering between the trans- 
verse pump beam and the cavity mode couples these two 
momentum states via two distinguishable Raman chan- 
nels, resulting in a dipole-type interaction between cav- 
ity mode and the corresponding collective spin degree 
of freedom, see Eq. (67). The parameters {luo,ijj,X) of 
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FIG. 36 (Color online) Dicke model phase diagram, (a) Dis- 
played is in color the recorded mean intracavity photon num- 
ber n as a function of the transverse pump power P and the 
pump-cavity detuning Ac. A sharp phase boundary is ob- 
served in agreement with a mean-field description (dashed 
line). The dispersively shifted cavity resonance for the non- 
organized BEC is indicated by a horizontal arrow. (b,c) Time 
traces of n while gradually increasing the pump power to 
1.3 mW for the indicated pump-cavity detuning. From Bau- 
mann et al. (2010). 



the corresponding realization of the Dicke Hamiltonian 
are given by the energy difference 2ujji betvifeen the two 
momentum modes (neglecting atom-atom collisions), the 
effective detuning —Sc between the pump laser frequency 
and the dispersively shifted cavity mode frequency, and 
the collective two-photon Rabi frequency be- 
tween pump laser and cavity mode, respectively. In the 
experiment (Baumann et al., 2010), 5c exceeds the recoil 
frequency by three orders of magnitude, thus realizing 
the dispersive regime of the Dicke model. Higher-order 
momentum modes do not contribute in the phase transi- 
tion dynamics itself and are populated only when deeply 
entering the self-organized phase (Konya et al., 2011). 

From the analogy to the Dicke model the following 
expression for the critical coupling strength is obtained 
upon including cavity decay (Dimcr et al., 2007) 



A. 



(68) 



In the absence of atom-atom collisions, this condition 
agrees with the result obtained from the stability analy- 
sis of the mean- field equations, Eq. (66). Experimentally 
(Baumann et al., 2010), the phase boundary was mapped 
out as a function of pump-cavity detuning in agree- 
ment with the theoretical prediction, see Fig. 36. 

It is instructive to contrast the Dicke quantum phase 
transition, realized with a BEC in a single-mode cavity, 
with the occurrence of free-space superradiant Rayleigh 
scattering off an elongated BEC which is driven by off- 
resonant laser light (Inouye et al., 1999). In this ex- 
periment, a superradiant light pulse was emitted along 



the axial direction of the atomic cloud accompanied by 
the creation of recoiling matter-wave components, once 
the pump intensity exceeded a critical value. This dy- 
namical effect is equivalent to Dicke superradiance of a 
collectively excited medium (Dicke, 1954), where matter- 
wave amplification phase-locks the spontaneous emission 
events into the continuum of optical field modes. The 
minimal pump intensity required for superradiance to oc- 
cur is determined by the balance between loss and gain 
processes. In contrast, light scattering off the BEC into a 
single cavity mode is a reversible process and the critical 
pump strength dominantly results from the finite pump- 
cavity detuning. 

In the self-organized phase the detuning 6c between 
pump laser and the dispersively shifted cavity resonance 
becomes a dynamic quantity. This effect is not taken 
into account by the description in terms of the Dicke 
model Eq. (67) which is a valid approximation as along 
as the maximum dispersive cavity shift UqN is small 
compared to the pump-cavity detuning |Ac|. In the 
case where U^N exceeds |Ac| the system can exhibit 
dynamically frustrated behavior characterized by a pe- 
riodic sign change of the effective pump detuning from 
the dispersively shifted cavity resonance, as was observed 
in Baumann et al. (2010), see Fig. 36c. Theoretically 
the influence of the additional nonlinear dispersive term 
^ UoJza^ci appearing in the Dicke model Eq. (67) upon 
the dynamics of the coupled BEC-cavity system was 
investigated in great detail by Bhascen et al. (2012); 
Keeling et al. (2010b); and Liu et al. (2011). Employ- 
ing a semiclassical description, Bhaseen et al. (2012) re- 
veal a rich phase diagram including distinct superradiant 
fixed points, bistable and multistable coexistence phase 
and regimes of persistent oscillations, and explore the 
timescales for reaching these asymptotic states. It is em- 
phasized that the behavior of the open system is con- 
trolled by the stable attractors, which do not necessarily 
coincide with the points of minimal free energy. As such, 
there is a crucial distinction between the k limit of 
the dynamical system and the equilibrium behavior at 
K = 0. Similar conclusion can be drawn in the quantum 
case, as discussed in the following. 

The coupling of the cavity field to the electromag- 
netic field environment, causing cavity decay, amounts 
to a weak measurement of the coupled BEC-cavity sys- 
tem. The corresponding quantum back-action results in 
a diffusion-like depletion of the ground state of the Dicke 
Hamiltonian, even at zero temperature. The underly- 
ing physics is similar to that described in section IV. C. 3, 
with the important difference that the system gets in- 
creasingly susceptible to quantum back-action when ap- 
proaching the critical point. 

The rate at which the ground state of the Dicke Hamil- 
tonian initially gets depleted due to cavity decay was 
calculated in Nagy et al. (2010) based on the Langevin 
equation approach Eq. (58). In the dispersive regime 
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FIG. 37 (Color online) Criticality in the closed and open- 
system Dicke phase transition. The mean values (thick, right 
axes) and the incoherent excitation numbers (thin, left axes) 
of the photon (a) and atomic (b) fields are plotted as a func- 
tion of relative coupling strength jz/^cr ~ A/Acr. The incoher- 
ent excitation numbers in steady-state for the open system 
(thin solid) diverge at the critical point with exponent —1 
in contrast to the ground-state number of excitations for the 
closed system (thin solid), which diverges with the mean-field 
exponent —1/2. From Nagy et al. (2011). 



\5c\ ^ Wij, the ground state depletion happens mostly 
in the atomic space and the corresponding diffusion rate 
can be approximated below threshold by the expres- 
sion oj /\5c\{^/ y^cT:)"^ ■ Per atom, this corresponds for 

2 

\5c\ ^ K to a heating rate of Note, the formal 

"c 

equivalence of this result with the spontaneous heating 
rate in a far-detuned dipole trap. Importantly, the use 
of a large detuning \5c\ removes the time limitation im- 
posed by measurement-induced back-action. 

Measurement-induced back-action drives the BEC- 
cavity system into a steady state which is a dynamical 
equilibrium between diffusion and damping. It is very 
interesting that this limiting state is not the same as 
the equilibrium state of the system, i.e., for k — Q the 
ground state at T = 0. Namely, the order of the two 
limiting procedures, i — > oo and k — > cannot be in- 
terchanged. The steady-state (t — >■ oo limit) occupation 



of the cavity field and the excited momentum state were 
calculated in Nagy et al. (2011). Flux and second-order 
time correlations of the cavity output signal were investi- 
gated theoretically in Oztop et al. (2012). The mean-field 
obtained in the thermodynamic limit is a smooth func- 
tion of K, and the steady-state solution tends to that of 
the ground state of Eq. (67) for k — >• 0. By contrast, 
the comparison of the quantum fluctuations present in 
the ground state of the Dicke model and in the steady- 
state of the damped-driven system exhibit a significant 
difference. In the ground-state of Eq. (67) the second- 
order correlation functions diverge towards the critical 
point with the exponent —1/2, indicating a mean-field- 
type transition, whereas in the non-equilibrium case the 
quantum fluctuations exhibit a divergence with exponent 
— 1 (see Fig. 37). At the same time, the singularity of 
the ground-state entanglement (Lambert et al., 2004) be- 
tween cavity and atomic subsystem is regularized at the 
critical point by the quantum noise associated with cavity 
decay. The non- vanishing entanglement, however, shows 
that the quantum character of the Dicke quantum phase 
transition (self-organization at zero temperature) is not 
fully destroyed in case of an open-system dynamics, and 
it cannot be exactly mapped to a thermal noise-driven 
phase transition. 



3. Phases in highly degenerate cavities 

Self-organization of polarizable particles into periodic 
structures induced and stabilized by the intracavity light 
field resembles the process of crystallization. In a cavity 
with only a single standing-wave mode tuned into reso- 
nance with the pump field, only the amplitude of the cav- 
ity field is a dynamical quantity. The formation of a pe- 
riodic crystal from a spatially homogeneous distribution 
breaks the discrete symmetry corresponding to the even 
and odd antinodes of the standing-wave mode profile. 
Already in the two-mode setting of a ring cavity sustain- 
ing two degenerate counter-propagating modes, the self- 
organization is accompanied by spontaneous breaking of 
a continuous translational symmetry (Nagy et al., 2006b; 
Nicdenzu et al., 2010), which induces rigidity against lat- 
tice deformations. In the case of highly degenerate multi- 
mode cavities the field has much more freedom to adjust 
locally to the particle distribution. 

The general structure of the resulting complex phase 
diagram was studied in Gopalakrishnan et al. (2009) and 
Ritsch (2009). In general such setups allow to realize 
conceptually novel systems and to explore and discover 
properties of crystalline and liquid-crystalline ordering, 
including intrinsic effects as dislocations, the growth and 
arrangement of crystal grain boundaries (see Fig. 38), 
and the nature of the phonon spectrum. Multimode cavi- 
ties also offer a natural connection to models developed in 
the field of neural networks as the Hopfield model or sim- 
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ilar spin models with infinite range statistical couplings. 
First ideas about this relation have been raised very re- 
cently in Gopalakrishnan et al. (2011a, b); and Strack 
and Sachdcv (2011). Extensions to fermionic atoms in 
multimode cavities have been considered in Miiller et al. 
(2012). 




FIG. 38 (Color online) Self-ordered states in a concentric 
multimode cavity forming two-dimensional patterns. The di- 
agram shows a regime near threshold, with domains locally 
populating distinct TEM^xy cavity modes in the equatorial 
plane. Domains can be punctuated by dislocations (shown 
in the left half of the figure), but might also show textural 
variation in space (right half of the figure). The black lines 
represent nodes of the cavity field, which separate 'even' and 
'odd' antinodes. As the atoms are Bose-condensed, the atomic 
population per site is not fixed. From Gopalakrishnan et al. 
(2009). 

Gopalakrishnan and coworkers generalized and 
adapted a field-theoretical framework, also successfully 
used in solid-state physics, to describe many-body 
systems coupled to a multitude of degenerate modes of a 
high-finesse cavity (Gopalakrishnan et al., 2010; Keeling 
et al., 2010a). For a quasi- two-dimensional cloud of 
atoms confined in the equatorial plane of a concentric 
optical cavity, the transition from the homogeneous 
distribution into a spatially modulated one is of the 
Brazovskii type (Brazovskii, 1975), which describes the 
phase transitions from isotropic to striped structures 
in liquid crystals. The description is based on an 
effective equilibrium theory which is valid when the 
effective cavity loss rate K?7^g/A^ is smaller than the 
recoil frequency w]^^, see IV. D. 2. Here, the dispersive 
cavity shift was assumed to be much smaller than the 
pump-cavity detuning Ap. Unlike the Landau theory 
of crystallization, here the free energy of the system 
does not have a cubic term that breaks the symmetry 
at the phase transition. The transition persists at zero 
temperature, hence it realizes a quantum phase transi- 
tion of an unusual university class. The non-equilibrium 
extension of this theory, which includes the effect of 
photon leakage out of the cavity as a perturbation, 
leads to the conclusion that the photon loss corresponds 
to an effective temperature and quantum correlations 
are washed out by decoherence on timescales longer 
than the cavity decay time. Note that the Bogoliubov- 
type mean-field model description of the open-system 
Dicke-model in a single mode cavity also predicts the 



depletion of the ground state (Nagy et al, 2010) due 
to measurement- induced back-action (Murch et al., 
2008), or, in other words, due to the diffusion induced 
by fiuctuations accompanying cavity photon loss (Nagy 
et al., 2009). However, unlike the Brazovskii transition, 
the Dicke-model system is driven into a steady-state, 
significantly different from the ground state, which has 
a critical point (Nagy et al., 2011). 



single mode case 





multimode case a , 




a . 










y SS ><; 








\ uniform 
normal\ \ superfluid 
solid / \. 


K.E./C 




/ ~ supersolid/ ^^^^ 1 (Braz 


ovsldi, weakly 




/ supersmectic '^"^^^ 


1st order) 




K.E.IC 



FIG. 39 Schematic zero-temperature phase diagram for a 
BEG in a concentric cavity. The control parameters are the 
atomic scattering length a and the inverse effective atom- 
cavity coupling (^~^ with — -q^ff/Ac. For weak, repulsive in- 
teractions and increasing atom-cavity coupling, the superfluid 
first undergoes self-organization via the Brazovskii transition, 
thus forming a supersolid. If the transverse laser intensity is 
increased further, the supersolid undergoes a transition into 
a normal solid (i.e., a Mott insulator). For strong, repulsive 
interactions, the uniform BEG can lose phase coherence con- 
currently with a first-order self-organization transition. This 
situation is to be contrasted with that for the case of a single- 
mode cavity (inset), in which there should always be a su- 
persolid (SS) region separating the uniform fiuid (SF) and 
normal solid (S) regions. First- and second-order transitions 
are marked by 1 and 2, respectively. From Gopalakrishnan 
et al. (2010). 

Similarly to the case of the single-mode experiment 
performed in the Zurich group (Baumann et al, 2010), 
the emergent crystalline state in a transversally driven 
multimode cavity can be considered as a supersolid phase 
where crystalline order and off-diagonal long-range order 
(long-range phase coherence) coexist. The phase dia- 
gram, schematically shown in Fig. 39 for a multimode 
cavity is strikingly different from the single-mode case 
in that a region with direct uniform superfluid-to-normal 
phase transition occurs, whereas in the single-mode cav- 
ity there is always a supersolid state between the uniform 
and the normal solid phases (Vidal et al., 2010). It is also 
observed that, for a strongly layered three-dimensional 
structure, the inter-layer frustration precludes global or- 
dering and the system breaks up into inhomogeneous do- 
mains. 
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E. Extended Hubbard-type models for ultracold atoms in 
cavities 

The theoretical description of the quantum many-body 
dynamics of ultracold atoms confined in optical lattices 
and strongly interacting with a quantized cavity field 
can be based on a sophisticated extension of the Bose- 
Hubbard (BH) model (Fisher et al, 1989). In static op- 
tical lattices, the BH model properly accounts for the 
quantum statistical properties of bosonic atoms at the 
lattice sites, as well as the inter-particle quantum corre- 
lations (Bloch et ai, 2008). The basic assumption, valid 
in the limit of very low temperature, is that the dynam- 
ics can be restricted to the lowest (or lowest few) Bloch 
bands of the periodic optical potential. Correspondingly, 
the many-body wavefunction can be expressed in terms 
of Wannier functions localized at individual lattice sites. 

However, if the optical lattice potential is sustained 
by the mode of a high-finesse cavity, thus becoming a 
dynamical degree of freedom, it gets a highly nonlin- 
ear problem to determine the Wannier functions them- 
selves and thereby the ground state of the many-body 
system. For a laser-driven cavity, e.g., the atoms disper- 
sively shift the cavity resonance, thus affecting the intra- 
cavity field amplitude which itself determines the lattice 
depth. Hence, the optical lattice potential and the state 
of the atoms have to be evaluated in a self-consistent 
way (Larson et ai, 2008a; Maschler et ai, 2008; Maschler 
and Ritsch, 2005; Vidal et al, 2010), as will be presented 
below. We focus on the most studied particular case 
of spinless bosons, and mention only that the cases of 
fermions and spin particles are expected to lead to in- 
teresting novel effects (Larson et al., 2008c; Sun et al., 
2011). 



1. Bose-Hubbard model with cavity-mediated atom-atom 
interactions 

Consider an ensemble of N bosonic particles subject to 
an optical lattice potential which is generated by the field 
of an optical resonator and possibly by an additional, 
far off-resonant standing-wave laser field. The latter is 
represented by the external potential term Vci{r) in the 
single-atom Hamiltonian Eq. (51). Restricting the mo- 
tional dynamics to the lowest energy band (lowest vibra- 
tional state), we expand the atomic field operator 



M 

*(r) = ^6,w(r-r,), 



(69) 



in the Wannier basis of atomic states localized at sites 
i = 1 . . . M, where bi denotes the associated annihilation 
operator. Upon inserting this expansion into Eqs. (49), 



one obtains 
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where several cavity modes with mode functions /m(r) 
and corresponding photon annihilation operators 
considered. The coefficients Eij and Jfj are defined as 
in the standard BH Hamiltonian 

— j™(r-r,), (71a) 

d^rw{r ~ Ti) fci{r)w{r - Tj) , (71b) 



Jf 



where we separated the characteristic magnitude Vd of 
the classical trapping potential (difference between max- 
imum and minimum) from its spatial form /ci(r). The 
last term of Eq. (70) describes the on-site interaction 
with U = {inash^ /m) J c?^r|w(r)|^. Of primary interest 
are the extra couplings generated by the cavity modes 
with matrix elements 



47 = / d'vw{v - r,)/;(r)/„(r)zi;(r - r,) , 



(71c) 



From the cavity field point of view, the diagonal elements, 
I = m, correspond to the atomic state-dependent disper- 
sive shifts of the cavity mode frequency wc.mj whereas 
the off-diagonal elements, I ^ m, describe photon scatter- 
ing between different cavity modes. The Wannier func- 
tions w(y — Yi) appearing in these integrals, in principle 
depend on the dynamic potential terms generated by the 
cavity field. This renders the problem highly non-trivial. 

In the most general case, the Wannier functions have 
to be calculated for each photon number state to define a 
corresponding manifold of parameters in the BH model, 
Eq. (70). In other words, the couphngs J-'J' are replaced 
by operators which can be easily expressed in a Fock 
basis. Such a brute force approach is necessary if the 
effect of the cavity field on the trapping potential is sig- 
nificantly different for adjacent Fock states (Horak and 
Ritsch, 2001b). Typically, numerical simulations have to 
be performed to study, e.g., microscopic processes under- 
lying many-body effects that are understood in the mean- 
field limit (Maschler et al, 2007; Niedenzu et al, 2010; 
Vukics et al, 2007). Obviously, this approach is limited 
to small system sizes of a few particles moving in a few 



45 




L?/ /probe flo 



yprol 

f 



FIG. 40 (Color online) Scheme for quantum non-demolition 
measurement of atomic many-body states in an optical lattice. 

atoms are trapped in a one-dimensional lattice potential 
(green, M lattice sites) which partially overlaps with a cavity 
mode Hi (blue) and a transverse probe mode ao (red). The 
number of illuminated lattice sites is denoted by K. Depend- 
ing on their many-body state the atoms act as a quantum 
refractive index whose statistical distribution with respect to 
the probe and cavity modes can be mapped out via transmis- 
sion or diffraction spectroscopy as a function of the probe- 
cavity detuning or the angles 0o and Q\. 



cavity modes. Most of the works, however, used approx- 
imations to treat the cavity-generated optical potential 
in which the localized Wannier functions are defined in a 
self-consistent manner. 



2. Cavity-enhanced light scattering for quantum measurement 
and preparation 

Before addressing the problem of dynamical cavity- 
induced potentials within the framework of the BH 
model, we note that a lot of applications have been 
developed based on the coupling of quantized cavity 
field modes to trapped, ultracold atomic systems in a 
simple scattering regime, as was exhaustively reviewed 
by Mekhov and Ritsch (2012b). In the scattering sce- 
nario, the external lattice potential Ki(r) is taken strong 
enough to define solely the localized Wannier functions, 
and their modification due to the cavity light forces is 
negligible. The quantized cavity field modes are a pertur- 
bative probe which can yield a mapping between quan- 
tum properties of atomic many-body states and light ob- 
servables. This system gives, for example, a means to 
determine the quantum state of ultracold atoms by light 
scattering (Miyake et ai, 2011). 

A typical quantum measurement scheme involving a 
single cavity mode is depicted in Fig. 40. It was shown 
that various quantum states of ultracold bosons trapped 
in the lowest band of an optical lattice and having equal 
mean densities can be distinguished (Chen et ai, 2007; 
Mekhov et ai, 2007a). As a characteristic example, the 
very different transmission spectra of the Mott insula- 
tor (MI) and the superfluid (SF) states are exhibited in 
Fig. 41. In contrast to standard techniques this measure- 
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FIG. 41 (Color online) Cavity transmission spectra showing 
the atom number distribution of an ultracold gas in the in- 
tracavity part of an optical lattice, see Fig. 40. Shown are 
transmission profiles as a function of the probe-cavity detun- 
ing Ap of a MI (red) and a SF (blue) state for (a) a good 
cavity with k = 0.1?7ii and (b) a bad cavity with k = Uu, 
where Uu = gl/AA- In (b) the satellites are not resolved but 
the spectra for SF and MI states are still different. Parame- 
ters are A'^ = M = 30 and K = 15. (c) Spectra for a SF state 
with N — AI = 70 and different number of illuminated lattice 
sites K. From (Mekhov et ai, 2007c). 



ment is non-destructive, limited only by quantum mea- 
surement back-action. Depending on the chosen geome- 
try, light scattering is sensitive to the global and local 
atom number fluctuations (Bhattacherjee et al., 2010; 
Chen and Meystrc, 2009; Mekhov et ai, 2007c), or to 
long-range correlations between two or more lattice sites 
(e.g. four-point correlations) (Bux et ai, 2011; Mekhov 
et al, 2007b). 

The back-action of repeated quantum non-demolition 
measurements drives the atomic many-body state to- 
wards specific states (Mekhov and Ritsch, 2009). This 
is in full analogy to microwave cavity QED experiments 
(Brunc et ai, 1990; Gueiiin et ai, 2007) where, comple- 
mentarily, the field state is driven into nonclassical Fock 
states by measuring the state of a train of atoms crossing 
the cavity. 



3. Self-consistent Bose-Hubbard models in cavity mean-field 
approximation 

Genuine cavity-induced dynamical effects appear when 
the excited cavity modes significantly modify the trap- 
ping with respect to the external potential Vci(r). Even 
without noticeably reshaping the Wannier-functions at 
the trapping sites, the perturbative light probe can mod- 
ify the tunneling rates as was shown for free-space Bragg 
scattering by Rist et al. (2010). In a cavity, local changes 
of the atomic distribution influence the whole cavity- 
sustained optical lattice potential. Thereby a new type of 
long-range interaction between the particles appears and 
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gives rise to resonant nonlocal co-tunneling or momen- 
tum space pairing (Mekhov et ai, 2007b), effects which 
go far beyond the standard Bose-Hubbard model. 

When the Wannier functions themselves are dynam- 
ically influenced by the cavity, a self-consistent mean- 
field approach, similar to the one in Sec. IV. B, has 
been broadly adopted to describe the nonlinear dynam- 
ics of trapped, ultracold atoms in a cavity (Chen et aZ., 
2009; Larson et ai, 2008b; Larson and Lewenstein, 2009; 
Maschler et ai, 2008; Maschler and Ritsch, 2005; Nimm- 
richter et al, 2010; Vidal et ai, 2010). Splitting the 
cavity field amplitude a = (a) + Sa into its mean value 
and fluctuations, the main assumption is that only the 
highly excited mean field can modify the trapping poten- 
tial, whereas the fluctuations amount to a perturbative 
probe. Since the cavity mean-fleld amplitude (a) depends 
on the momentary atomic quantum state, so does the 
depth and shape of the on-site potential. Therefore, the 
Wannier functions in Eq. (69) and hence the coefficients 
Eq. (71) in the Hubbard-type Hamiltonian Eq. (70), have 
to be determined self-consistently in conjunction with the 
proper cavity mean field (a). It is noteworthy that, often, 
the self-consistent calculation does not lead to a unique 
solution. 



a. Phases in dynamical optical lattices Consider a linear 
cavity with only a single mode being driven and over- 
lapping with a static optical lattice potential Vci{r). As- 
suming only nearest-neighbor hopping to be relevant we 
keep adjacent pairs from the double sum over in- 

dices i and j in Eq. (70). The corresponding many-body 
Hamiltonian is given by 

H = EqN + Ei3+ {hUoa^a + Ki) [JqN + JB 

-nAca'a-inr][a-a^) + ^C. (72) 

The relevant atomic degrees of freedom are the total atom 
number N and the collective nearest-neighbor coherence 
B defined as 



M 



M 



(73) 



respectively, and the operator C = ^ ■ b\bi {h\hi — for 
the two-body on-site interaction. The coefficients Eq, 
E, Jq and J derive from Eq. (71) contracted to a single 
cavity mode and assuming uniform coupling along the 
lattice. 

To exhibit the underlying physics one may neglect the 
photon number dependence of the Wannier functions and 
adiabatically eliminate the cavity field via the Heisenberg 
equation of motion 



\ n —1 


(a) 


\ "o-^ 


A =2k 

c 

TT —-9 w 


















d = |i Ac-Uo (^JqN + JB^ - k| 



a + 1] . 



(74) 



FIG. 42 Phase diagram with overlapping Mott insulator 
states. Boundaries of different Mott lobes (shaded regions) 
as a function of the rescaled chemical potential p. and the 
inverse of the pump strength rj (in units of k) in the ID 
cavity lattice potential of K = 50 sites. Parameters are (a) 
(Ac,r/o) = (2k, -2k) and (b) {Ac,Uo) = (0,2k). The Mott 
lobes are labeled by the number of atoms per site no. The 
dashed Hnes show the boundaries of zones which are hidden. 
Outside the shaded parameter regions, the state of the system 
is superfluid in most cases. From (Larson et al, 2008a). 



This is a good approximation as long as the cavity field 
decays fast compared to the timescale of atomic motion. 
Since tunneling in deep lattice potentials is slow com- 
pared to the recoil frequency, this applies widely in ex- 
perimental setups. 

An effective atomic Bose-Hubbard model, formally 
identical to the usual one but with coefficients Jq and J 
depending on the many-body state, has been considered 
systematically in the thermodynamic limit. By evaluat- 
ing the stability of the Mott insulator states, the phase 
diagram has been constructed. As shown in Fig. 42, the 
model predicts the existence of competing Mott insula- 
tor states (Larson et ai, 2008a,b). The overlapping Mott 
lobes indicate the possibility of bistability in this laser- 
driven, nonlinear system, see IV. C. 2. The state of the 
system can be controlled by fine tuning the pump pa- 
rameters near the shifted cavity resonance. For certain 
parameters a state with two atoms per site can lead to 
a much higher photon number and thus deeper optical 
potential, so that its energy falls below that of the state 
with unity filling. 
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In order to gain insight into the nature of atom-atom 
couphng via the cavity field, a simple effective Hamilto- 
nian can be constructed. The adiabatic field amplitude 
can be expanded to second order in the small tunneling 
matrix element J 



V 

K — i6c 



Up J ^ {Up J? 
;-i6c {k-iSc)' 



, (75) 



where the effective detuning 6c ~ Ac — UqJoN was in- 
troduced, and the atom number was set to N. Inserting 
this solution back into the Hamiltonian Eq. (72) and the 
Liouville operator Eq. (5a), which accounts for cavity 
damping, leads to an effective adiabatic model. It com- 
prises the nonlinear Hamiltonian 
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and Liouville operator 



2BgI3 - I3^g- gI3^ 



(76a) 



(76b) 



which describes decoherence in the basis of the eigen- 
states of the operator B. Note that the adiabatic elimi- 
nation procedure described above is not rigorous mathe- 
matically, since we adiabatically approximated the solu- 
tion of a nonlinear dynamical equation in Eq. (75), which 
appears also as an ordering ambiguity of the involved op- 
erators (Larson et ai, 2008b; Maschler et ai, 2008). 

For a small, numerically tractable system, the lowest 
energy eigenstate of H^d can be calculated. As a key 
example, the SF-to-MI quantum phase transition in an 
optical lattice sustained entirely by a quantized standing- 
wave cavity field was analyzed (Maschler et al., 2008). 
The dynamical response of the photon number to the 
atomic motion is able to strongly modify atomic num- 
ber fluctuations and hence to drive the phase transition. 
Depending on the cavity parameters (e.g. the detuning 
between cavity and external pump laser), the photon fluc- 
tuations can either suppress or enhance atomic fluctua- 
tions and hopping, therefore pushing the system towards 
or outwards the MI or SF states. Accordingly, as depicted 
in Fig. 43, the position of the SF-to-MI phase transition 
in a cavity optical lattice potential can be shifted (keep- 
ing the mean potential depth constant) towards either 
smaller or larger values of the coUisional atom-atom inter- 
action strengths, depending on whether the pump-cavity 
detuning is chosen positive (a) or negative (b). 

Figure 44 demonstrates the importance of the photon 
number fluctuations in a quantum potential by testing 
the stability of the Mott phase in a potential of fixed av- 
erage depth but different mean photon number. While 
in an almost classical field (highly excited coherent state 




FIG. 43 Mott insulator (MI) to superfluid (SF) phase tran- 
sition in a cavity optical lattice. The probabilities pmi and 
PSF to find the atoms in the states I^Pmi) and I'I'sf) as a 
function of the dimensionless ID on-site interaction strength 
giD / {dEiij [d is the lattice constant, En — Tiuir is the recoil 
energy) are compared for two cases: first, for an optical lattice 
sustained by the quantum field of a cavity mode (Vci = 0), 
and second, for a purely classical optical lattice (771 = 0). We 
choose r] such that in each of the two examples (a and b) 
both potentials have equal depth for zero on-site interaction 
giD. The quantum (QM) and classical (class) cases are de- 
picted with solid and dashed lines, respectively. Parameters 
are {Uo,K,ri) = {— 1,1/^/2, y/K5)u!R. The detuning between 
the probe and the dispersively shifted cavity frequencies af- 
fects the position of the phase transition. In (a) resp. (b) 
this detuning is positive, Ac — UoN — k, resp. negative 
Ac — UqN — —K and the transition point is shifted towards 
smaller resp. higher interaction strengths in comparison to 
that in a classical lattice. From Maschler et al. (2008). 



with many photons) the Mott phase is stable, photon 
number fluctuations (uncertainty) inherent in a weak co- 
herent state of few photons enhance tunneling and decay 
of the perfect order. 

This is explicitly shown in Fig. 44 depicting the decay 
of an initially prepared perfectly ordered atomic state. 
We choose different mean intracavity photon numbers 
and keep the average depth of the potential constant by 
readjusting the couphng strength Uq- In the classical 
limit (very large photon number and small atom-cavity 
coupling), the system remains in the initially prepared 
MI state. For smaller photon numbers n ~ 20, the ini- 
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FIG. 44 (Color online) EfTect of photon number granularity 
upon the Mott insulator (MI) state for two atoms in a cavity- 
sustained optical lattice. The time evolution of the occupation 
probability pmi of the MI state is shown for various mean 
intracavity photon numbers n. The atom-cavity coupling g 
is adjusted such that the average potential depth (SEr) is 
identical for all curves. Whereas for a static optical lattice 
potential of equal depth the system remains in the initially 
prepared MI state (solid line), photon number fluctuations for 
low n deplete the MI state. From Maschler et al. (2008). 
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FIG. 45 Self-organization as a quantum seesaw effect, (a) 
Atoms which are trapped in a free-space ID lattice poten- 
tial with two adjacent sites (left and right) tunnel with rate 
J between the corresponding Wannier states \l) and |r). (b) 
Coupling the atoms in addition to the field of a cavity whose 
axis is perpendicular to the optical lattice induces light scat- 
tering between the optical lattice laser and the cavity field 
with opposite phases from the two sites. The modified po- 
tential resulting from the interference of the lattice field and 
the cavity field discriminates the two sites and causes positive 
feedback and atomic ordering into one of them. The process 
starts by spontaneous symmetry breaking, and depends on 
the quantum statistics of the initially prepared many-body 
state. 



tial MI state only slowly degrades in time. However, 
when the photon number fluctuations become compara- 
ble to the mean, i.e. for mean photon numbers as small 
as n ~ 1, the system quickly escapes from the MI state 
via fluctuation induced tunneling. Note that in order to 
keep the average optical potential constant, lower pho- 
ton numbers are connected to a larger potential per pho- 
ton, so that the potential fluctuations are additionally 
enhanced at low photon numbers. The classical limit is 
also approached in the bad cavity limit, k ^ uir, where 
number fluctuations occur so fast, that particles only see 
the average and do not have the time to tunnel during 
an intensity fluctuation. 

The quantum properties of the cavity light become pre- 
dominant if already single intracavity photons create an 
optical potential of considerable depth, capable of trap- 
ping numerous atoms. As quantum mechanics allows for 
the existence of superpositions of photon number states, 
one may obtain superpositions of several potentials with 
different depths (Horak and Ritsch, 2001b). 

The complete phase diagram of ultracold atoms in two- 
band BH models coupled to a cavity light field has been 
calculated by Silver et al. (2010) by means of a vari- 
ational approach and the analogy to the Dicke-model 
superradiant phase transition has been pointed out, see 
Sec. IV.D.2. 



b. Self-organization within the Hubbard-model approach A 
possible influence of quantum statistical properties on the 



spatial self-organization process, described in Sec. III.B.l 
and IV. D, can be studied within the framework of an 
extended BH model. For simplicity, the geometry is 
modified in comparison with the generic case of self- 
organization, as depicted schematically in Fig. 45. Atoms 
are confined in a static optical lattice potential which 
is oriented perpendicular to the cavity axis. As before, 
large atom-laser detuning and negligible atomic satura- 
tion are assumed. The laser fields providing the optical 
lattice potential are considered to be tuned close to res- 
onance with a cavity mode, therefore inducing coherent 
cavity driving via Rayleigh scattering off the atoms. The 
single-atom Hamiltonian corresponding to this geometry 
reads (Maschler et a/., 2007) 

2 

H = ^ + Vci cos^(kx) -h(Ac- Uo) a'' a 
2 m 

-I- \/fiVciUo cos{kx) (a + a^) , (77) 

where Vd denotes the depth of the static lattice poten- 
tial, Ac the detuning between the lattice laser and the 
cavity resonance, and Uq the light shift of the cavity res- 
onance frequency per atom. In a simple and intuitive 
picture the dynamic cavity field plays the role of a "see- 
saw" potential. Interference between the cos{kx) poten- 
tial, generated through photon scattering, and the static 
cos^ {kx) potential determines the overall potential felt by 
the atoms. The spatial symmetry of the system allows 
for the emergence of two possible ordered configurations, 
with all atoms residing at either odd (cos(fcx) = 1) or 
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FIG. 46 Entanglement-assisted self-organization in a quan- 
tum optical lattice potential, see Fig. 45. Shown are atom- 
light entanglement (solid lines), mean cavity photon number 
(dotted-dashed lines), and the two-site atom-atom correla- 
tion function (dotted lines) for two atoms in two adjacent 
lattice wells, denoted by left (l) and right (r). Lines with ex- 
tra circles show the case of exactly one atom in each well at 
start (MI state), while the other lines show the evolution for 
an initially symmetric superposition state for each atom (SF 
state). The parameters are Uo = —2k, Ac = —6k, J = k/100 
and J = 1.6«:. From Maschler et al. (2007). 

even (cos(fca;) = —1) lattice sites. 

The many-body Bose-Hubbard Hamiltonian, Eq. (70), 
adapted to this scheme reads 

^ = Y. J^'^^l^^ -hlAc-UoY. bib, J a^a 

+ {a + a^)J2^J,Abr (78) 

Here, the standard matrix elements for the kinetic and 
potential energy /2m + V'cicos^(fcx) between sites i 
and j are denoted by Jij, whereas J^.j gives the ma- 
trix elements of the interference term •\/?7oKi/?icos(fcx). 
On-site interactions are neglected at this point and, for 
consistency, we require weak couphng per atom, i.e. 

n\Uo\ < |Fci|. 

Essential dynamical properties of this system beyond 
the mean-field approximation become evident already for 
two atoms. Monitoring the microscopic physics of self- 
organization as shown in Fig. 46, the process resembles 
the decay of a homogeneously filled lattice with one par- 
ticle per site on average to the self-ordered state, where 
both particles occupy even or odd sites. Following the 
decay of the probability for the two atoms sitting in dif- 
ferent wells {{nirir)), we first note that the formation of 
the self-organized state is accompanied by a fast growth 
of atom-field entanglement. Most importantly, however, 
one finds a striking dependence of the self-organization 



dynamics on the initial quantum fluctuations. A SF state 
with both atoms prepared in the symmetric superposition 
of the two wells, 1/2 (6| -I- bl)'^\0), self-organizes much 
faster than a perfectly ordered MI state, bjbl\0), with 
exactly one atom per well. In the latter case, the cavity 
field remains in the vacuum state until a tunneling event 
induces atomic coherence between the left and right lat- 
tice site, triggering the decay of the MI state towards the 
self-organized state. 

Under some approximations this model can be extrap- 
olated to the thermodynamic limit where quantum phase 
transitions similar to the one predicted in the mean-field 
approach can be studied. Among various other properties 
this leads to the coexistence of diagonal long-range order 
and long-range coherence (Vidal et at, 2010), indicat- 
ing new phases to appear in the gaps between Mott-like 
states with different integer filling factors. 



c. Ring cavity When several independent cavity modes 
are dynamically interacting with the atoms, not only the 
depth but also the shape and the spatial periodicity of 
the potential can change. In the generic case of a ring 
cavity the depth and the longitudinal position of the lat- 
tice is dynamical. While, already in standard optical lat- 
tices the validity of the lowest-band assumption is often 
doubtful and corrections are necessary, this approxima- 
tion loses its meaning in a ring cavity. Expansions based 
on a single set of Wannier functions cannot be consis- 
tently formulated since the lowest-band Wannier func- 
tions for a given position contain contributions from a 
large number of higher bands for a slightly shifted posi- 
tion. Hence small lattice shifts immediately involve many 
higher-order bands. 

The simplest example of two cavity modes sustained by 
a ring cavity reveals, that a naive crude truncation of the 
Bose-Hubbard model with respect to the pumped cosine 
mode at the lowest band decouples the atoms from the 
associated sine mode. This immediately eliminates cen- 
tral dynamical effects of the system as overall momentum 
conservation and nonlocal correlated hopping (Nicdenzu 
et al, 2010). 

As an example. Fig. 47 shows that a quantum jump 
in the lattice photon number is accompanied by a sud- 
den change of the tunneling rate between adjacent sites. 
After the second jump in the trajectory shown in the 
Figure, the system returns to the original mean values, 
however, both the position and photon number quanti- 
ties exhibit much larger noise, which demonstrates the 
effect of heating stemming from the photon number fluc- 
tuations. 
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FIG. 47 (Color online) Correlated photon jumps and tunnel- 
ing of an atom in a symmetrically driven ring cavity. The sam- 
ple trajectory, showing the position expectation value (kx) 
(left axis) and the mean photon number {a' a) of the un- 
driven sine mode (right axis), presents two quantum jumps 
occurring at URt ~ 70 and URt ~ 100. They lead to a 
simultaneous change of the photon number and the atomic 
band excitation corresponding to different tunnel oscillation 
frequencies. In a higher band much faster oscillations are 
obtained between neighboring sites. This generates a higher 
effective hopping amplitude and heating on average. The pa- 
rameters are Uq — ~2u}ji, — \/6 (field amplitude of the 
driven cosine mode), Ac = f/o — ^ and k — 500aj_R. From 
Niedenzu et al. (2010). 



V. OUTLOOK 

After almost two decades of active research in optical 
cavity QED with cold and ultracold atoms, initially dom- 
inated by theoretical investigations, the field presently 
exhibits fast growth with several experimental groups 
demonstrating spectacular effects. Single atoms are rou- 
tinely cavity-cooled and trapped over seconds within op- 
tical high-finesse resonators, providing a well-controlled 
quantum system for quantum information science. Ul- 
tracold quantum gases prepared in magnetic or opti- 
cal traps are now reliably coupled to high-quality cavity 
modes. Even in the far dispersive regime, these systems 
are governed by strong back-action effects of the collec- 
tive atomic motion on the cavity field degrees of freedom. 
Cavity decay offers a unique channel to monitor the com- 
plex coupled atom-light dynamics non-destructively and 
in real time. Many central atomic variables can be ac- 
cessed by quantum non-demolition measurements which 
minimize quantum back-action. 

With the aim of trapping and cooling ensembles, nano- 
and even microscopic particles or arrays of thin mem- 
branes, the research field of cavity QED, reviewed in 
this paper, overlaps and unifies more and more with the 
rapidly growing field of optomechanics (Stamper-Kurn, 
2012). Trapping and cooling arrays of membranes in 
cavities is only one very striking example (Xuereb et ai, 
2012). Practical applications as ultra-sensitive detectors 



of mass, acceleration or magnetic fields or even tests of 
general relativity seem within range of current technol- 
ogy- 
Initiated by theoretical studies and early experiments 
with cold atomic ensembles, a non-equilibrium quantum 
phase transition between a superfluid and a supersolid 
phase has been investigated experimentally. Recently, 
theoretical investigations opened new directions and pos- 
sibilities towards controlled preparation and investiga- 
tions of the physics of spin glasses, more complex su- 
persolid and superglass phases (Gopalakrishnan et at, 
2011b; Strack and Sachdev, 2011). Further prominent 
solid-state Hamiltonians involving phonons or polarons 
could be studied with unprecedented control and obser- 
vation possibilities (Mekhov and Ritsch, 2012a). 

A recent breakthrough experiment demonstrating sub- 
recoil cavity cooling towards quantum degeneracy (Wolke 
et ai, 2012) opens the prospect of replacing evaporative 
cooling techniques by cavity cooling and direct prepa- 
ration of exotic quantum states from a thermal gas. 
This also paves the way towards implementing a contin- 
uous atom-laser as a new tool in ultracold atom physics 
(Salzburgcr and Ritsch, 2007). 

Still, important challenges as cooling and trapping of 
molecular samples or large suspended objects have not 
been experimentally demonstrated. The prospects of 
multi-species implementations in multimode cavity en- 
vironments still have to be fully evaluated. Apart from 
this point, experiments seem ahead of theoretical and nu- 
merical simulation possibilities, where theory has to be 
improved and better suited models need to be developed. 

In a long term vision, cavity-sustained light fields 
allow to couple hybrid systems of very different physical 
nature like superconducting qubits, cold quantum gases 
and micromechanical oscillators without destroying 
quantum coherence of the systems, brought about 
by any classical coupling of such systems. In this 
way cavity-based setups with ultracold gases could 
develop into an important building block for quantum 
information processing or other quantum-based future 
technologies (Hcnschel et ai, 2010) or a route to an even 
better atomic lattice clock (Nicholson et ai, 2012). 
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